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A  microcontroller-based  system  of  oceanographic  instrumentation  providing  a 
comprehensive  set  of  measurements  relevant  to  sediment  transport  processes  has  been 
developed.  Analysis  of  the  data  provided  by  the  system  yields  time  series  of  vertical 
profiles  of  mean  sediment  size  and  concentration,  horizontal  profiles  of  bedform 
geometry,  and  single  location  measurements  of  flow  velocity,  pressure,  turbidity,  and 
water  temperature.  Details  of  the  system  architecture,  including  capabilities  provided  by 
both  hardware  and  software  contained  within  the  system  are  given. 

An  improved  method  for  the  determination  of  suspended  sediment  size  and 
concentration  from  the  system's  acoustic  backscatter  intensity  measurements  is 
presented.  By  retaining  the  size  dependence  throughout  the  derivation  for  an  explicit 
solution  for  concentration,  a  new  explicit  solution  to  the  acoustic  backscatter  equation 


results.  This  new  concentration  solution  improves  the  technique  for  determining  median 
sediment  size  by  incorporating  sediment  attenuation  in  the  calculation.  Because  this  new 
technique  relies  on  the  minimization  of  the  variance  in  concentration  as  determined  by 
different  frequency  transducers,  the  previous  technique  of  pairing  transducers  of  different 
frequencies  is  replaced  by  a  technique  making  use  of  any  number  of  different  frequency 
transducers.  The  new  size/concentration  inversion  technique  is  tested  using  both 
simulated  and  laboratory  data.  Numerical  precision  is  shown  to  be  the  only  source  of 
error  with  the  use  of  simulated  data.  Laboratory  tests  result  in  less  than  20%  error  in  the 
determination  of  both  concentration  and  size  over  a  range  of  nearly  one  meter. 

Finally,  suspended  sediment  concentration  data  from  the  nearshore  region 
obtained  from  an  experiment  performed  in  Duck,  North  Carolina,  are  examined  to  find 
the  relevant  time  scales  of  sediment  suspension.  In  this  location,  low  frequency  forcing 
mechanisms  are  as  significant  in  suspending  sediment  as  the  incident-band  wave  forces 
typically  used  to  model  suspension.  Like  wave  groups,  this  low  frequency  forcing  results 
from  the  linear  superposition  of  velocity  components  in  a  narrow  band  of  frequencies. 
When  these  frequency  interactions  are  considered,  coherence  greater  than  60%  is  found 
between  the  velocity  squared  and  the  near-bed  concentration  across  most  of  the  spectrum. 


CHAPTER  1 
INTRODUCTION 

With  continuing  coastal  development,  the  need  for  proper  coastal  planning  and 
management  of  its  resources  grows.  Towards  meeting  this  need,  several  models  of  beach 
evolution  are  currently  used  as  tools  for  predicting  shoreline  change  and  for  design  of 
coastal  structures.  Improving  these  models  requires  a  better  understanding  of  the 
underlying  physical  processes  of  sediment  transport.  The  capability  of  a  theory  to 
describe  the  basic  physics  of  a  process  is  dependent  on  the  accuracy  of  measurements  of 
the  process.  Refinement  of  the  theory  of  the  basic  mechanics  of  sediment  transport  thus 
depends  upon  improved  measurements  of  the  processes. 

Developments  in  computers  and  instrumentation  have  not  only  improved  the 
capability  of  measurement  of  sediment  transport  phenomena  in  the  laboratory,  but  have 
made  possible  field  measurements  not  possible  a  decade  ago.  While  the  control  and 
repeatability  of  laboratory  experimentation  is  appealing,  there  are  factors  which  influence 
sediment  transport  in  the  field,  which  are  neglected  in  laboratory  studies.  The  evidence 
of  this  is  the  disagreement  often  encountered  between  results  from  laboratory  and  field 
experiments.  Thus,  it  is  desirable  not  only  to  verify  laboratory  results  in  the  field,  but 
also  to  improve  field  experimentation  apparatus  and  technique  such  that  other  elements 
influencing  the  motion  of  sediment  can  be  found. 

In  this  dissertation,  a  new  system  of  instrumentation  developed  for  field 
measurement  of  small-scale  sediment  transport  and  the  associated  hydrodynamics  is 
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presented.  Chapter  2  begins  the  presentation  with  details  of  the  instrumentation  and 
architecture  of  the  most  recently  developed  system.  This  system  was  developed  after 
several  refinements  of  earlier  systems.  Each  refined  system  was  developed  in  response  to 
newly  discovered  limitations  of  the  prior  system.  Improvements  from  earlier  versions 
include  new  and  better  instrumentation,  faster  communications,  and  less  restrictive  data 
storage  requirements.  In  all,  five  systems  were  developed,  beginning  with  the  system 
used  in  the  Supertank  project  (Hanes  et  al.,  1993;  Thosteson,  1995).  Subsequent  systems 
were  used  for  the  Vilano  Beach  project,  the  Duck94  project,  and  the  SIS95  and  SIS96 
projects.  During  the  writing  of  this  text,  the  latest  system  is  being  used  in  the 
SandyDuck97  project.  In  later  chapters,  the  data  presented  are  not  that  from  the  system 
described  in  Chapter  2,  but  instead  are  from  the  next  most  recent  version  of  the  system. 
In  an  effort  to  present  the  reader  with  the  most  recent  and  useful  information,  Chapter  2 
presents  the  newer  system.  The  overall  architecture  has  changed  little  since  the  version 
used  in  the  Duck94  experiments,  but  the  reader  will  be  alerted  to  differences  in  the 
systems  when  appropriate. 

Perhaps  the  most  powerful  instrument  and  certainly  the  most  demanding 
computationally  on  the  system  is  an  acoustic  backscatter  system,  used  for  obtaining  high 
resolution  measurements,  both  spatially  and  temporally,  of  suspended  sediment 
concentration  and  size.  Recent  improvements  to  the  technique  of  determining 
concentration  from  the  measured  backscattered  intensities  are  further  extended  in  chapter 
3  to  include  determination  of  the  mean  grain  diameter  of  particles  in  suspension.  The 
new  technique  presented  removes  attenuation  assumptions  previously  required  in  order  to 
determine  the  particle  size  and  reduces  the  computational  effort  in  determination  of  both 


concentration  and  size.  Also,  this  technique  allows  more  than  two  unique  frequency 
transducers  to  be  used  in  a  single  computation.    This  reduces  the  possibility  of  obtaining 
multiple  solutions  for  size,  which  can  arise  when  only  two  unique  frequencies  are  used, 
and  again  reduces  computational  effort. 

In  chapter  4,  the  field  experiments  performed  the  fall  of  1996  at  the  SIS96  project 
are  reviewed  and  the  analysis  of  the  resulting  data  is  covered.  Field  application  of  the 
acoustic  analysis  techniques  presented  in  chapter  3  for  determination  of  suspended 
sediment  size  will  be  investigated.  The  remainder  of  the  analysis  is  centered  on 
determining  the  frequency  ranges  in  which  sediment  suspension  is  most  prevalent.  It  is 
seen  that  often  the  most  dominant  suspension  events  occur  with  a  frequency  much  lower 
than  the  frequencies  of  the  incident  waves,  though  very  little  energy  is  found  in  the 
spectrum  within  the  infragravity  band.  A  coherency  function  analysis  with  the  square  of 
the  bottom  velocity  shows  that  the  low  frequency  forcing  responsible  for  the  sediment 
suspension  results  from  the  interaction  between  frequencies  in  the  incident  band  of  the 
spectrum.  Interaction  between  frequency  components  is  reviewed,  showing  that  the 
resulting  low  frequency  forcing  is  caused  by  wave  groups.  In  chapter  5,  conclusions 
from  this  research  and  suggestions  for  future  research  will  be  presented. 


CHAPTER  2 
SYSTEM  DEVELOPMENT 

With  more  sophisticated  oceanographic  instrumentation  comes  the  need  for  more 
advanced  data  acquisition  and  instrument  control  systems.  Evolving  from  acquisition 
systems  used  in  several  earlier  experiments,  a  new  system  has  been  developed  for  use  in 
the  SandyDuck97  project,  held  at  the  Field  Research  Facility  (FRF)  in  Duck,  North 
Carolina  in  late  1997  (Hanes  et  al.,  1993).    This  new  system  delivers  highly  accurate 
measurements  of  both  the  hydrodynamics  and  the  resulting  sediment  processes  with  high 
temporal  and  spatial  resolution.  The  focus  during  the  development  of  this  system  has 
been  on  the  ability  to  collect  an  unbroken  record  of  data,  consisting  of  highly 
synchronous  data  from  all  instruments  attached  to  the  system,  with  duration  of  several 
days.  Such  a  record  allows  the  investigators  to  utilize  standard  time  series  analysis 
techniques  to  examine  sediment  processes  on  time  scales  ranging  from  seconds  to  days 
without  making  assumptions  typically  required  to  account  for  gaps  between  data  records. 
In  its  current  form,  outlined  in  figure  2-1,  the  system  utilizes  a  three  frequency  acoustic 
backscatter  system  (ABS),  two  acoustic  Doppler  velocimeters  (ADVs),  three  multiple- 
transducer  arrays  (MTAs),  an  optical  backscatterance  sensor  (OBS),  a  pressure  sensor,  a 
compass,  a  tilt  meter,  an  external  temperature  sensor,  and  an  underwater  video  camera. 
Internal  monitoring  includes  leak  detection,  internal  temperature,  unregulated  and 
regulated  shore  supplied  voltages,  and  the  interna]  battery  voltage.  This  chapter 


Figure  2-1.  External  components  to  acquisition  system. 

describes  the  system  architecture,  essential  algorithms  responsible  for  data  collection  and 
processing,  and  the  system's  capabilities. 


2.1.  Instruments 
The  choice  of  instrumentation  used  by  the  current  system  was  made  to  most 
accurately  record  as  many  parameters  relevant  to  sediment  suspension  processes  as 
possible  with  minimal  flow  obstruction.  Waves  and  currents  are  measured  using  a 
Transmetrics  pressure  sensor  and  Sontek  three-axis  acoustic  Doppler  velocimeters 
(ADVs).  ADVs  were  chosen  due  to  their  ability  to  provide  high-resolution  velocity 
measurements  for  long  periods  of  time  without  the  zero-drift  difficulty  associated  with 
electromagnetic  devices  -  an  important  consideration  with  the  long  records  which  will  be 
collected  with  this  system.  Tilt  and  compass  transducers  will  ensure  accurate  positioning 
of  the  ADVs  and  monitor  any  frame  motion  or  failure.  For  a  combined  concentration  and 
turbidity  measurement,  a  single  D&A  Instruments  OBS  is  utilized.  Three  Seatek  MTAs, 


two  operating  at  2.00  MHz  and  the  third  at  5.00  MHz,  are  used  for  bottom  bedform  and 
slope  measurements.  MTAs  are  linear  arrays  of  acoustic  transducers  which  provide  high 
resolution  measurements  of  bedforms  and  the  local  slope  of  the  seabed.  The  5.00  MHz 
unit  contains  32  transducers  spaced  1.5  cm  apart,  making  it  ideal  for  measurement  of 
small-scale  bedforms  (wavelengths  of  order  5  to  50  cm).  Each  2.00  MHz  unit  contains 
16  transducers  spaced  6  cm  apart,  providing  excellent  measurement  of  larger  scale 
bedforms  (wavelengths  of  order  20  to  200  cm).  To  this  point  in  time,  all  three  MTAs 
have  been  used  together  to  form  a  single  2.5  m  array,  giving  bedform  measurements 
ranging  from  small  to  large  scale  (Jette  and  Hanes,  1997).  Since  orientation  of  bedforms 
can  not  be  determined  when  the  MTAs  are  used  in  this  manner,  orientation  is  determined 
from  video  observations  with  a  DeepSea  Micro-SeaCam  1050  when  turbidity  is  low  and 
from  an  independent  Simrad  Mesotech  model  900  rotating  scanning  sonar  system,  not 
described  in  this  text.  Profiling  for  suspended  sediment  size  and  concentration  is 
accomplished  using  an  acoustic  backscatter  system  (ABS)  made  by  the  Centre  for 
Environment,  Fisheries  and  Aquaculture  Science  (CEFAS,  formerly  known  as  The 
Fisheries  Laboratory),  a  description  of  which  follows.  Finally,  to  determine  speed  of 
sound  and  acoustic  attenuation  parameters  for  use  in  evaluating  data  from  the  acoustic 
instruments,  water  temperature  is  measured.  Further  information  on  the  instrumentation 
vendors  is  found  in  appendix  B. 

Simultaneous  operation  of  multiple  acoustic  concentration  profilers  at  different 
frequencies  has  successfully  provided  measurements  of  profiles  of  both  sediment 
concentration  and  mean  sediment  size  (Hay  and  Sheng,  1992).  The  present  system 
includes  an  ABS  utilizing  1 .0,  2.2  and  5.0  MHz  transducers.  Sync  and  trigger  lines 


provide  precise  control  over  selection  and  firing  of  the  transducers.  A  single  analog 
signal  line  returns  the  envelopes  of  the  backscattered  analog  voltages.  Upon  receipt  of  a 
pulse  on  the  sync  line,  the  ABS  begins  a  cycle  such  that  the  first  pulse  on  the  trigger  line 
fires  the  first  transducer  and  switches  the  signal  line  to  transmit  data  from  the  first 
transducer  to  the  acquisition  system.  The  next  pulse  on  the  trigger  line  causes  the  second 
transducer  to  be  fired  and  switches  the  signal  line  to  transmit  the  second  transducer's 
voltage  envelope.  The  cycle  continues  with  the  third  pulse  selecting  the  third  transducer. 
Additional  pulses  on  the  trigger  line  restart  the  cycle,  beginning  again  with  the  first 
transducer. 

Of  the  instruments  listed  above,  the  system  is  responsible  for  analog  to  digital 
conversion  and  subsequent  storage  of  the  data  from  all  but  the  MTAs  and  the  video 
camera.  Responsibility  of  the  system  to  the  MTAs  is  limited  to  setting  acquisition 
parameters  and  collecting  data  through  a  digital  serial  connection  shared  between  all 
three  MTAs  and  the  main  system.  The  video  camera  is  switched  on  and  off  and  provided 
power  by  the  system.  The  video  signal  simply  enters  the  underwater  package  and  exits 
immediately  through  the  shore  cable.  Of  the  remaining  instruments,  all  but  the  ABS  are 
instruments  that  collect  data  from  a  single  point  in  space  only.  For  this  reason,  they  will 
henceforth  collectively  be  referred  to  as  the  single  point  instruments  (SPIs).  The  system 
has  been  designed  to  provide  data  at  a  maximum  rate  of  four  final  measurements  per 
second.  Signals  from  each  of  the  SPIs  are  passed  through  two-pole,  linear-phase,  anti- 
aliasing filters  with  a  cut-off  near  10  Hz.  This  filter,  combined  with  over-sampling  of  the 
SPIs  at  100  Hz  and  subsequent  averaging  to  the  desired  final  measurement  rate, 
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eliminates  the  possibility  of  aliasing  and  minimizes  signal  degradation  due  to  the 
frequency  response  of  the  filter. 

2,2  System  Demands 

Fundamental  to  the  design  of  this  system  is  the  desire  to  minimize  the  travel 
distance  of  analog  signals  from  the  instruments.  The  purpose  of  this  is  to  minimize 
filtering  effects,  both  in  amplitude  and  phase,  and  reduce  the  introduction  of  additional 
signal  degradation  caused  by  the  otherwise  necessary  modulation  and  demodulation  of 
the  raw  high-frequency  signals.  We  therefore  place  the  system  offshore,  which  results  in 
new  design  demands  related  to  power  consumption,  size,  storage  requirements,  and 
system  monitoring  functions.  We  also  desire  that  the  system  be  fairly  self-contained, 
requiring  only  minimal  external  equipment  support  and  user-intervention.  User- 
intervention  consists  of  configuring  sampling  schemes  before  data  collection  and 
offloading  of  any  data  residing  on  the  system  after  collection  is  complete.  Finally,  given 
the  specific  modes  of  operation  of  the  instruments  to  be  used  in  the  study  as  well  as  any 
other  demands  imposed  by  them,  synchronous  measurements  from  all  instruments  is  of 
the  utmost  importance. 

Ideally,  each  of  the  three  ABS  transducers  should  measure  the  identical  profile  at 
the  same  instant  in  time,  giving  three  coincident  and  collocated  measurements  of  acoustic 
backscatter.  In  practice,  this  is  difficult  for  several  reasons.  First,  physical  limitations 
due  to  the  size  of  the  transducers  and  the  acoustic  beams  from  the  transducers  prevent 
achieving  measurements  from  perfectly  over-lapping  profiles.  Second,  to  statistically 
reduce  the  random  fluctuations  in  the  backscattered  signal  due  to  coherence  of  the 
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returned  signal,  for  each  transducer,  single  profiles  of  the  mean  power  are  constructed 
from  several  consecutive  profiles.  It  is  important  to  realize  that  the  mean  power  must  be 
used  rather  than  the  mean  voltage,  as  the  concentration  of  suspended  sediment  is 
proportional  to  the  power  of  the  backscattered  signal.  This  is  shown  by  the  acoustic 
backscatter  equation  written  in  terms  of  the  mass  concentration,  CfThorne  et  al.,  1991). 


c{z)  JHz)mzwYe^^.) 


(2-1) 


In  this  equation,  P  is  the  backscattered  pressure  (proportional  to  the  voltage  measured  at 
the  crystal),  z  is  the  distance  from  the  transducer,  tyMs  a  function  describing  near-field 
characteristics  of  the  transducer,  k(l  is  the  system  constant,  and  a„.  and  a,  are  water  and 
sediment  attenuation  parameters  respectively.  The  system  constant  characterizes  specific 
physical  and  electrical  properties  of  a  particular  transducer.  These  parameters  and  the  use 
of  the  acoustic  backscatter  equation  are  examined  in  more  detail  in  chapter  3  and 
appendix  A.  Since  the  sediment  size  and  concentration  profiles  are  in  general  functions 
of  time,  additional  fluctuations  in  the  backscattered  signal  due  to  random  reconfiguration 
of  the  particles  in  the  ensonified  volume  and  due  to  changing  sizes  and  concentrations 
certainly  result.  To  account  for  these  limitations,  several  conditions  shall  be  imposed  on 
the  acquisition  system.  Profiles  must  be  acquired  at  a  high  enough  rate  such  that  the 
changes  with  time  in  size  and  concentration  profiles  are  small.  Within  this  period,  a 
sufficient  number  of  profiles  need  to  be  collected  such  that  a  representative  average  will 
statistically  minimize  the  fluctuations.  Finally,  the  ensemble  of  samples  from  any  one 
transducer  should  be  well  distributed  over  this  entire  period. 
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Acoustic  concentration  profilers  very  quickly  produce  large  volumes  of  data.  To 
meet  the  three  conditions  imposed  above,  the  current  system  collects  100  profiles  per 
second  from  each  of  the  three  transducers  on  the  ABS.  Each  individual  profile  consists 
of  up  to  120  12-bit  samples  from  the  analog  signal,  sampled  at  100  kHz.  For  ease  of 
processing  and  storage,  each  12-bit  sample,  hereafter  referred  to  as  a  bin,  is  stored  in  a 
zero-padded  16-bit  word.  As  each  profile  is  collected,  the  squared  value  of  each  bin  is 
calculated  using  a  highly  efficient  table-lookup  algorithm.  This  squared  value  is  then 
added  to  a  running  sum  of  squares,  such  that  for  each  bin  in  the  profile,  the  root  mean 
squared  (RMS)  value  can  be  calculated  after  a  specified  number  of  profiles  have  been 
collected.  Again  in  the  interest  of  minimizing  computational  time,  RMS  values  are 
calculated  from  the  sums  through  the  use  of  a  binary  search  method  through  a  previously 
calculated  table  of  sums.  Figure  2-2  outlines  the  entire  reduction  procedure  for  the  case 
of  computing  RMS  profiles  of  120  points  from  23  consecutive  profiles  -  a  common 
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Figure  2-2.  Real-time  processing  of  collected  back  scattered  signal  profiles. 
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configuration  of  this  system.  Algorithms  for  computing  the  square  and  the  RMS  will  be 
detailed  later  in  this  text.  The  final  outcome  of  this  on-the-fly  processing  of  the 
backscattered  signal  profiles  is  a  set  of  statistically  meaningful  profiles  requiring 
significantly  less  storage  space  and  significantly  less  communications  bandwidth  to 
transmit  real-time  digital  data  from  the  system. 

2.3    System  Architecture 
To  meet  the  control,  communications,  synchronization,  processing,  and  data 
collection  requirements  given,  a  low-power  Onset  Tattletale  model  7  data  logging  engine 
is  employed.  The  model  7  provides  512  kilobytes  of  flash  EEPROM  storage  and  256 
kilobytes  of  static  RAM  for  system  and  application  software,  as  well  as  an  additional  2 
megabytes  of  pseudo-static  RAM  for  data  storage.  In  addition,  it  includes  28  digital  I/O 
lines,  a  4  channel  analog  to  digital  converter,  a  parallel  I/O  port,  RS-232 
communications,  a  real-time  clock,  and  500  megabytes  of  hard  drive  storage.  Central  to 
the  model  7  is  the  Motorola  68332  microcontroller.  The  68332  incorporates  the  CPU32 
central  processing  unit  executing  a  superset  of  MC68000  instructions.  For  use  in  the 
present  system,  the  microcontroller  operates  at  16  MHz.  Also,  the  68332  incorporates  a 
time  processing  unit  (TPU),  which  is  essentially  a  special  purpose  slave  processor  that 
controls  two  timers  and  sixteen  I/O  lines.  In  the  present  system,  this  TPU  is 
indispensable  for  handling  synchronization  as  well  as  hardware  and  software  triggering. 
Figure  2-3  outlines  the  additional  interfacing  made  to  the  model  7  in  the  system. 
Additional  information  and  descriptions  follow. 


12 


Further  extending  the  features  of  the  model  7  is  a  Daedulus  Research  MUX32 
board.  This  board  further  multiplexes  each  of  the  four  analog  to  digital  channels  by 
eight,  allowing  the  model  7  to  sample  up  to  32  analog  signal  lines.  These  32  channels 
will  be  referenced  by  a  port,  numbered  0  to  3,  corresponding  to  the  original  channel 
number  on  the  analog  to  digital  converter,  and  a  channel  number  numbered  0  to  7,  now 
referring  to  one  of  the  newly  multiplexed  lines.  Besides  increasing  the  number  of  analog 
lines,  the  MUX32  also  buffers  the  16  TPU  digital  I/O  lines  and  replaces  the  model  7's 
analog  to  digital  reference  with  a  voltage  reference  of  higher  precision. 


com.  to  MTA's 


power  from  shore 


com.  to  shore 


Figure  2-3.  Block  diagram  of  internal  system  components. 
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Interfaced  to  the  model  7  is  a  filter  board,  which  provides  the  anti-aliasing  for  the 
analog  signals,  as  previously  mentioned,  a  power  board,  and  several  communications 
components.  On  the  shore  end,  the  DC  voltage  supplied  to  the  package  is  such  that  30  V 
is  input  to  the  system  package  offshore  with  all  instruments  turned  on.  For  example,  if 
the  resistance  in  500  meters  of  cable  is  such  that  15  V  is  lost  due  to  this  resistance,  45  V 
would  be  supplied  from  shore.  Given  30  V  input,  the  power  card  regulates  power  for  all 
components  of  the  system  in  addition  to  providing  power  switching  for  each  of  the 
instruments  independently.  Also,  the  power  card  keeps  24  V  of  NiCd  batteries  charged, 
for  use  in  supplying  short  bursts  of  high  current  to  power  the  model  7's  hard  disk  drive. 
These  batteries  also  can  power  the  system  for  a  short  time  in  the  event  the  power 
connection  from  shore  is  compromised.  Communications  from  the  package  to  shore  is  at 
57.6  kilobaud  through  two  50  ohm  coaxial  cables  using  RS-422  transceivers  on  both  ends 
of  the  connection.  Such  a  connection  has  been  shown  in  previous  experiments  to  be 
reliable  through  500  meters  of  coaxial  cable.  For  reasons  discussed  shortly,  RS-232 
communications  from  the  model  7  go  through  an  Onset  Tattletale  model  8.  Between  the 
model  7  and  model  8,  serial  communications  are  at  38.4  kBd.  The  serial  communications 
path  then  goes  from  the  model  8's  57.6  kBd  RS-232  connection  to  the  RS-422  transceiver 
on  the  communications  card. 

By  default,  the  main  clock,  TCR 1 ,  on  the  TPU  of  the  68332  runs  at  one  fourth  of 
the  CPUs  clock  frequency.  Instead,  with  a  CPU  clock  frequency  of  16  MHz,  TCR1  is  set 
during  system  initialization  to  operate  at  1  MHz  due  to  the  range  of  clock  frequencies  that 
must  be  generated  by  the  TPU.  Without  making  this  change,  the  risk  of  losing 
synchronization  between  channels  because  of  counter  overruns  within  the  TPU  increases. 
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TPU  channels  are  used  in  the  present  system  for  hardware  synching  and  triggering  of  the 
ABS.  hardware  triggering  of  the  analog  to  digital  converter,  and  for  generation  of  the 
interrupts  necessary  for  the  triggering  of  software  exceptions  used  for  data  acquisition. 
By  using  a  mode  of  operation  on  the  TPU  known  as  synchronized  pulse  width 
modulation,  several  TPU  lines  are  set  to  run  continuously  as  perfectly  synchronized 
clocks.  Figure  2-4  illustrates  the  clock  rates,  pulse  widths,  and  timing  between  channels. 

Once  these  clocks  have  been  started,  they  will  run  continuously  and  perform  their 
associated  tasks  without  any  intervention  from  the  CPU.  Analog  signals  are  continuously 
sampled,  so  the  program  simply  decides  which  samples  to  keep.  Data  acquisition  is 
performed  entirely  within  exceptions,  described  in  more  detail  later,  so  it  can  be  started 
by  simply  enabling  the  interrupts  generated  by  TPU  channels  9  through  1 1 . 

Since  both  ACP  profiles  and  SPI  data  are  collected  at  100  Hz  and  subsequently 
processed  down  to  the  desired  final  data  acquisition  rate,  at  some  instant  immediately 
after  processing,  a  large  burst  of  data  is  both  saved  to  the  data  memory  within  the  model 
7  and  output  from  the  model  7.  This  burst  of  data,  up  to  approximately  1000  bytes  of 
data,  must  be  transmitted  from  the  model  7  within  roughly  10  ms  to  prevent  interfering 
with  the  timing  of  the  data  acquisition.  Rather  than  complicating  the  communications 
requirements  by  use  of  a  MBd  serial  line,  these  data  are  simply  output  to  the  existing 
parallel  port  on  the  model  7  into  a  FIFO  buffer  and  then  into  the  model  8  mentioned 
above.  By  continuously  monitoring  the  parallel  and  serial  connections  to  the  model  7,  the 
model  8  buffers  all  communications  and  allows  the  use  of  a  single,  57.6  kBd,  serial 
connection  to  the  package.  This  57.6  kBd  line  is  sufficient  to  allow  transmission  four 
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times  per  second  of  120  point  profiles  from  each  of  the  three  ACPs  and  the  additional  14 
channels  of  interest. 

2,4  Software  and  Algorithms 

A  program  called  Monlog  1.0  controls  the  system.  Monlog  is  written  mainly  in 
C,  with  additional  assembly  coding  of  the  actual  data  acquisition  and  data  processing 
routines  for  speed.  Most  of  the  code  handles  user-interface  functions,  such  as  creating 
menus,  handling  user  I/O,  setting  the  time,  setting  record  numbers,  defining  the  sampling 
scheme,  and  offloading  stored  data  records.  The  algorithms  necessary  to  handle  these 
functions  are  trivial  and  will  not  be  covered  in  this  text,  although  a  full  listing  of  the 
program  is  given  in  appendix  C. 

Typical  operation  of  the  system  begins  at  power-up.  The  system  software  handles 
initialization  of  the  core  components  of  the  Model  7.  To  customize  specific  portions  of 
the  initialization,  such  as  the  necessary  modification  of  the  TPU  configuration  register  to 
establish  the  TPU's  1  MHz  clock  rate,  the  contents  of  the  Model  7's  serial  EEPROM  are 
modified.  Immediately  following  the  system  initialization,  control  is  given  to  Monlog. 
Monlog  starts  with  system  specific  initialization,  enabling  charging  of  the  system 
batteries  and  ensuring  that  all  interfaced  instruments  are  started  in  a  known  state. 
Communications  baud  rates  are  established  and  system  variables  initialized.  The  TPU 
channels  are  set  up  and  started  as  was  shown  in  figure  2-4,  with  the  interrupts  they 
generate  initially  disabled,  and  appropriate  default  values  are  set  to  govern  data 
acquisition.  Finally,  the  internal  monitoring  channels  are  sampled,  and  the  results,  along 
with  the  main  system  menu  are  displayed  to  the  user. 
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Once  the  main  menu  is  displayed,  several  options  are  available  to  the  user.  A 
standard  input  routine  is  used  throughout  the  program,  such  that  a  timeout  feature  can  be 
utilized  at  each  opportunity  for  input.  This  guarantees  that  any  involuntary  selections, 
such  as  could  be  caused  by  line  noise,  will  not  leave  the  program  in  an  indeterminate 
state.  Upon  any  timeout,  the  program  returns  control  to  the  main  menu  routine, 
maintaining  original  values  for  any  values  that  may  have  inadvertently  been  changed. 
Included  in  the  main  menu  are  options  that  allow  the  user  to  power  each  of  the 
instruments  interfaced  to  the  system  individually  and  to  subsequently  test  them  in  an 
interactive  test  mode,  displaying  results  to  the  user  in  real  time.  Additional  options  are 
used  to  offload  data  resident  on  the  systems  hard  drive  to  shore  using  the  Xmodem  data 
transfer  protocol  and  to  define  the  sampling  scheme  and  system  modes  of  operation. 

An  emergency  storm-mode  can  be  enabled  from  the  main-menu.  When  enabled, 
Monlog  will  detect  loss  of  communications  while  in  the  main  menu  routine,  and  should 
such  a  condition  be  detected,  the  program  will  begin  acquiring  data  using  a  conservative 
pre-defined  acquisition  scheme.  Since  loss  of  communication  most  likely  will  coincide 
with  loss  of  shore-supplied  power,  this  scheme  is  defined  for  collection  of  short  data  files 
at  regularly  spaced  intervals  in  an  attempt  to  conserve  battery  power  for  a  time  period 
roughly  equal  to  that  of  a  storm. 

Preparation  for  acquisition  begins  with  the  user  setting  the  time  and  date  and 
defining  the  sampling  scheme.  Parameters  such  as  the  number  of  points  per  ABS  profile, 
number  of  profiles  to  include  in  an  ensemble  RMS  profile,  final  sampling  rates  for  the 
ABS,  SPIs,  and  MTAs,  and  sampling  durations  are  all  defined  by  the  user.  Based  on  this 
information,  Monlog  allocates  all  available  data  memory  and  sets  up  pointers  to 
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appropriate  memory  locations  to  create  data  headers,  store  data,  create  look-up  tables  for 
use  by  the  squaring  and  RMS  algorithms,  and  temporarily  store  sums  during  data 
acquisition.  Figure  2-5  shows  the  structure  of  these  pointers  in  memory  based  on  values 
determined  from  the  user  configuration.  Note  that  the  ranges  labeled  TMU  and  PAD 
together  form  a  Mathwork's  Matlab  data  file  in  the  memory  of  the  model  7.  Such 
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Figure  2-5.  Memory  map  of  data  memory  on  the  model  7. 
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compact  structuring  beforehand  allows  one  to  offload  the  entire  range  of  memory 
sequentially  to  the  shore-based  computer  and  immediately  load  the  data  record  into  the 
analysis  software  for  data  inspection  and  processing.  Note  also  that  although  the  MTA 
data  are  not  collected  directly  by  the  model  7,  space  is  allocated  within  the  data  file  for 
each  of  the  MTAs.  MTA  data  are  offloaded  from  each  of  the  MTAs  after  the  rest  of  the 
data  file  has  been  filled.  This  data  format  is  convenient  since  the  data  from  all  of  the 
instruments  are  stored  in  the  same  file  along  with  both  the  starting  and  ending  times  and 
dates  of  the  collection.  In  prior  systems,  all  data  collected  by  the  system  resided  in  these 
data  files.  With  the  added  ability  to  collect  unbroken  data  files  of  length  greater  than 
what  could  possibly  be  stored  in  this  memory  segment,  only  the  first  portion  of  the 
collected  data,  whatever  portion  will  fit  within  the  2-megabyte  limit,  is  stored  in  this 
format. 

After  the  dynamic  memory  configuration,  all  that  remains  for  Monlog  to  do  is 
enable  the  interrupts  generated  by  the  TPU  channels  and  monitor  the  acquisition.  Four 
exception  routines  handle  all  of  the  data  acquisition.  Outside  of  these  four  exceptions, 
the  main  process  (that  process  running  prior  to  the  CPU's  receipt  of  any  interrupt) 
monitors  the  progress  of  the  exception  routines  and  at  the  appropriate  time,  performs  the 
data  analysis  and  storage.  Following  the  timing  diagram  shown  in  figure  2-4,  the  order  of 
events  following  the  memory  initialization  is  as  follows: 

1)  Interrupts  generated  by  TPU  channels  9  to  1 1  are  set  to  execute  the  same 
exception  routine,  entitled  NOTHING. 

2)  Interrupts  are  enabled. 
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3)  Upon  identifying  an  interrupt  from  TPU  channel  9,  that  interrupt  associated  with 
the  acquisition  of  transducer  F3,  the  NOTHING  routine  reassigns  each  interrupt 
its  own  exception  routine.  This  insures  that  a  profile  from  transducer  Fl  will  be 
the  first  to  be  sampled. 

4)  A  hardware  sync  pulse  is  sent  from  TPU  channel  12  to  the  ABS,  effectively 
resetting  the  ABS  to  sample  transducer  Fl  on  the  next  received  trigger  pulse. 

5)  A  hardware  trigger  is  sent  from  TPU  channel  13  to  the  ABS,  triggering  transducer 
Fl. 

6)  Shortly  after  sending  the  trigger  to  the  ABS,  TPU  channel  1 1  generates  an 
interrupt,  which  starts  the  exception  routine  responsible  for  sampling  transducer 
Fl. 

7)  As  data  are  continuously  being  sampled  from  the  analog  to  digital  converter  in 
relation  to  the  triggering  from  TPU  channel  14,  each  sample  is  acquired  from  the 
converter,  squared  by  use  of  a  previously  calculated  lookup  table,  and  added  to  a 
running  sum  appropriate  for  the  location  in  the  profile. 

8)  Transducer  Fl's  exception  routine  completes  and  acknowledges  the  interrupt, 
returning  control  back  to  the  main  process,  which  to  this  point,  continues  to 
monitor  the  progress  of  the  exceptions. 

9)  The  above  process  repeats  with  TPU  channel  1 3  triggering  transducer  F2,  and 
TPU  channel  10  generating  the  interrupt  that  executes  transducer  F2's  sampling 
routine. 

10)  Channel  10's  interrupt  is  acknowledged,  again  returning  control  to  the  main 
process. 
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1 1)  Again,  the  process  repeats,  with  TPU  channel  13  triggering  transducer  F3,  and 
TPU  channel  9  generating  the  interrupt  that  executes  transducer  F3's  sampling 
routine. 

1 2)  This  sampling  routine  additionally  selects  and  samples  each  SPI  attached  to  the 
system,  adds  the  result  to  an  appropriate  running  sum  for  the  particular  SPI,  and 
decrements  the  counters  being  monitored  by  the  main  process. 

1 3)  Once  the  main  process  has  detected  that  the  appropriate  number  of  profile 
acquisitions  have  occurred,  it  again  reassigns  the  interrupts  to  execute  the 
NOTHING  exception  routine.  This  routine  has  the  task  of  keeping  a  tally  of  the 
number  of  times  TPU  channel  9  interrupts  occur.  After  a  designated  number  of 
occurences,  this  routine  reassigns  the  interrupts  to  again  sample  analog  data. 

14)  Immediately  following  the  reassignment  of  the  interrupts  to  the  NOTHING 
exception,  the  main  process  performs  data  processing  and  storage.  Note  that  the 
NOTHING  exception  will  be  called  a  sufficient  number  of  times  to  allow  the  data 
processing  and  storage  tasks  to  complete  before  data  sampling  resumes. 

This  process  continues  until  the  time  limit  specified  in  the  sampling  scheme  definition  is 
surpassed. 

To  this  point,  the  details  of  the  data  analysis  have  been  neglected.  The  first 
analysis  algorithm  in  need  of  description  is  the  lookup  table  used  to  square  the  incoming 
samples  from  the  ABS.  Normal  multiplication  is  far  too  costly  in  processing  time,  so  a 
table  look-up  algorithm  is  employed.  This  table  is  simply  a  4096  element  (corresponding 
to  each  possible  value  from  the  12-bit  analog  to  digital  converter)  list  of  squared  values. 
Given  a  value  to  be  squared,  this  value  is  used  as  an  index  into  the  list,  or  table.  The 
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value  located  at  the  given  index  is  the  square  of  the  index.  Next,  the  second  algorithm  is 
responsible  for  providing  the  RMS  value  from  a  given  running  sum  of  squared  values. 
Again,  standard  algorithms  involving  multiplication  and  division  demand  too  much 
processor  time,  so  an  alternative  approach  is  used.  Since  the  number  of  elements  used  in 
the  calculation  of  the  RMS  value  is  known  in  advance  in  the  given  situation,  a  table  of 
sum  of  square  values  can  be  calculated  prior  to  data  acquisition.  If  N  elements  are  to  be 
used,  then  the  value  in  the  table  corresponding  to  the  ith  index  is  given  by  the  following: 
R,=  N{i +  0.5)1  (2-2) 

By  finding  the  lowest  number  in  the  table  greater  than  a  given  running  sum  of  squares 
value,  the  resulting  index  to  that  element  of  the  table  is  a  very  close  approximation  to  the 
RMS  value  for  the  given  sum.  For  example,  if  an  RMS  value  is  to  be  determined  from  23 
elements,  a  table  of  borderline  sums  of  23  squared  values  can  be  generated.  The  first 
element  in  the  table,  corresponding  to  index  0,  would  be  as  follows: 

ff„=23x(0.5)2=5.75  (2-3) 

If  a  given  sum  of  squares,  after  23  samples  have  been  added  to  it,  is  less  than  this  value, 
than  the  RMS  value  is  closely  approximated  by  the  index  into  the  table.  Finally,  since  the 
elements  of  the  table  are  ordered  by  value,  rather  than  searching  each  of  the  4096 
elements  of  the  table  for  the  first  value  greater  than  that  given,  a  binary  search  algorithm 
is  utilized.  The  search  begins  at  the  center  of  the  table,  and  the  values  are  compared. 
Based  on  the  results  of  this  comparison,  the  next  comparison  will  be  with  the  value 
centered  in  either  the  upper  or  lower  portion  of  the  table.  With  each  successive 
comparison,  the  number  of  table  elements  remaining  to  compare  is  cut  in  half.  Since  all 


23 

divisions  are  by  2,  binary  shift  operations  are  used  to  perform  the  division.  For  a  4096 
element  table,  only  12  comparisons  need  to  be  made.  The  assembly  code  used  for  this 
comparison  is  shown  in  figure  2-6.  Note  that  the  variables  used  in  the  beginning  of  the 
routine  as  memory  pointers  were  initialized  prior  to  this  segment  of  code. 
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;get  byte  from  sum(x"2) 
,-  initialize  sum(xA2)  storage 
; initialize  TEST 
; initialize  HALF 
; initialize  BASE 
; initialize  ANSWER 
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TEST=HALF  or  BASE 
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branch  back  1  final  time 
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Figure  2-6.  68000  series  assembly  necessary  for  binary  search  algorithm. 


2.5  Summary 
In  this  chapter,  a  system  of  oceanographic  instrumentation  capable  of  providing  a 
broad  set  of  measurements  relevant  to  sediment  transport  processes  is  presented.  All 
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instruments  are  in  close  proximity  to  the  data  system,  maximizing  data  integrity. 
Extremely  efficient  analysis  algorithms  written  in  C  and  assembly  language  provide  data 
processing  and  reduction  during  collection.  Use  of  multiple  processors  allows 
measurements  to  be  both  stored  internally  on  the  systems  internal  storage  and  exported  to 
shore  during  collection  through  serial  communications.  Data  stored  on  the  system  are 
offloaded  using  an  error  correction  protocol  to  a  compact,  binary  data  file  that  can  be 
immediately  viewed  and  interpreted  by  analysis  software. 


CHAPTER  3 
ACOUSTIC  BACKSCATTER  CONVERSION  TECHNIQUE 

Use  of  a  high  frequency  underwater  acoustic  transducer  to  profile  the  vertical 
distribution  of  sediment  concentration  has  been  demonstrated  under  both  laboratory  and 
field  conditions  by  several  investigators  (Hanes  et  al.,  1988;  Hanes  et  al.,  1993;  Hay  and 
Sheng,  1992;Thorneetal.,  1991;  Green  and  Vincent,  1991).  Basically,  profiles  of  the 
intensity  of  backscattered  sound  from  suspended  sediment  are  collected,  and  this  intensity 
is  related  to  the  suspended  sediment  concentration  at  each  point  in  the  profile.  A  typical 
technique  for  determination  of  concentration  from  acoustic  backscatter  data  requires 
inversion  of  the  acoustic  backscatter  equation  for  concentration,  which  yields  an  implicit 
equation  needing  an  iterative  technique  for  solution.  In  order  to  obtain  profiles  of 
concentration  from  a  single  transducer  operating  at  a  fixed  frequency,  knowledge  of  the 
sediment  grain  size  distribution  is  required  prior  to  applying  the  inversion  procedure.  In 
addition,  the  assumption  that  this  size  distribution  either  remains  constant  with  range  or 
with  a  predetermined  form  at  all  ranges  is  required,  as  the  scattering  and  absorption 
properties  of  the  sediment  are  dependent  upon  grain  size  (Thorne,  1993).  Recently, 
several  investigators  have  demonstrated  effective  measurement  of  both  the  profiles  of 
sediment  concentration  and  of  the  median  grain  size  of  the  distribution  (Hay  and  Sheng, 
1992;  Crawford  and  Hay,  1993).    The  technique  requires  the  use  of  multi-frequency 
acoustic  transducers.  Since  the  absorbing  and  scattering  properties  of  sediment  depend 
on  both  the  grain  size  and  upon  the  frequency  of  the  incident  sound,  each  unique 
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frequency  transducer  provides  an  independent  measurement  of  the  backscattered  intensity 
profile.  Hence,  each  point  in  the  profile  can  be  described  by  a  number  of  independent 
equations  equal  to  the  number  of  coincident  and  collocated  measurements  of  unique 
frequency.  Typically,  three  transducers  of  unique  frequency  are  used  to  collect 
coincident  intensity  profiles.  Although  in  theory,  use  of  three  frequencies  suggests  that  at 
each  measured  point,  concentration  and  two  parameters  of  the  grain-size  distribution  can 
be  determined,  typically,  only  concentration  and  one  parameter  of  the  distribution  are 
obtained.  The  sensitivity  of  the  equations  to  small  variations  in  intensity  and  also  the 
nonlinearity  of  the  sediment  size  response  functions  are  responsible  for  this  limitation. 

Lee  and  Hanes  (1995)  presented  an  explicit  solution  for  concentration,  to  be 
referred  to  in  this  paper  as  LH95,  from  the  acoustic  backscatter  equation,  significantly 
reducing  the  computational  effort  by  removing  the  need  for  iteration.  An  added  benefit 
derived  from  use  of  the  explicit  solution  is  removal  of  the  ambiguity  in  concentration 
solutions  obtained  by  the  iterative  solution  to  the  implicit  equation.  Figure  3-1  shows  the 
relation  between  concentration  and  transducer  voltage  (proportional  to  the  square  root  of 
intensity)  as  calculated  by  the  acoustic  backscatter  equation  (3-2)  assuming  a  constant 
concentration  profile.  From  this  figure,  it  is  apparent  that  a  single  voltage  value  from  the 
transducer  may  result  from  two  different  concentrations.  Physically,  this  can  be 
described  with  the  following  argument.  At  low  concentrations,  the  sound  attenuation  due 
to  sediment  in  the  sound  path  is  low,  resulting  in  an  increase  in  the  intensity  of  the 
backscattered  sound  with  increasing  concentration.  As  this  attenuation  becomes  more 
dominant,  the  intensity  of  the  backscattered  sound  begins  to  decrease  with  increased 
concentration.  So,  from  the  two  concentration  solutions  obtained  from  the  implicit 
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Theoretical  voltage  from  concentration;  5  MHz 


Dashed:  2  cm,  max  at  9.2  g/l 
Dotted:  16  cm,  max  at  1.3  g/l 
Solid:  31  cm,  max  at  0.7  g/l 


012345678 
Concentration  (g/l) 
Figure  3-1.  Expected  voltage  for  given  concentrations,  5  MHz  transducer. 
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equation,  the  choice  of  the  correct  solution  depends  upon  the  magnitude  of  the 
attenuation.  The  explicit  solution  for  concentration  yields  only  one  concentration,  as  this 
dependence  on  attenuation  is  accounted  for  by  integration  of  the  intensity  profile. 
In  spite  of  the  benefits  presented  by  use  of  the  LH95  explicit  solution,  it  is 
restricted  in  its  use  to  concentration  only.  Its  use  requires  that  the  grain  size  distribution 
be  known  and  constant  with  range.  Holdaway  and  Thorne  (1997)  extended  the 
functionality  of  the  solution  by  allowing  the  size  distribution  to  vary  but  to  retain  a 
predetermined  form  with  range  from  the  transducer.  In  the  following  derivation,  a 
similar  explicit  solution  for  concentration  retaining  the  dependence  on  the  grain  size 
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distribution  is  found  by  following  the  LH95  development.  Based  on  this  explicit  form, 
the  methods  used  for  evaluating  parameters  of  the  grain  size  distribution  are  re-examined, 
and  an  improved  method,  at  least  in  computational  effort,  is  introduced.  In  addition,  the 
existing  method  for  evaluating  the  grain  size  parameters  by  pairing  ACPs  of  different 
frequencies  is  extended  to  utilize  an  arbitrary  number  of  unique  frequency  transducers. 

3.1  Theory 
The  equation  that  relates  the  intensity  of  the  backscattered  acoustic  signal  to  the 
concentration  and  size  distribution  of  the  scatterers  in  suspension  is  referred  to  as  the 
acoustic  backscatter  equation.  This  equation  has  been  presented  in  several  forms,  each 
nearly  equivalent.  Presented  here  is  a  general  form  of  the  equation,  based  jointly  on  the 
form  presented  by  Hay  (1991)  and  Thorne(1993). 

,,j  2     lV2  _\     1      p/  srt  \  -4J(„„+„1)..sinh_B 

a,=-J^(z)C{z)dz  (3-2) 

The  variables  in  equations  (3-1  to  3-2)  are  defined  as  follows: 
V  =  voltage  read  from  transducer 
z  =  distance  from  transducer 
S  =  system  constant 

c  =  speed  of  sound,  assumed  constant  with  distance 
T  =  acoustic  pulse  width  (s) 
\jf  =  nearfield  correction  term  (see  appendix) 
F  =  backscatter  parameter  (see  appendix) 
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C  =  mass  concentration,  (sediment  density  in  F) 

a„.  =  water  attenuation  parameter  (see  appendix) 

as  =  sediment  attenuation  parameter 

B=(aw+f(z)c(?)) 

z0  =  range  from  transducer  at  which  first  concentration  and  size  is  evaluated. 

f  =  local  sediment  attenuation  proportionality  constant  (see  appendix) 
The  nearfield  correction  term,  i// ,  is  included  here  for  completeness,  and  it  can  easily  be 
included  in  the  following  derivation  by  temporarily  absorbing  it  in  the  backscatter 
parameter.  Since  it  is  simply  another  function  of  z,  it  does  not  complicate  the  derivation. 
It  has  not  been  included  in  the  following  derivation,  because  it  was  not  used  in  the 
subsequent  numerical  simulations  or  laboratory  tests.  The  complete  solution,  including 
this  term,  is  presented  in  the  appendix. 

In  equation  (3-1 ),  the  final  term  on  the  right  side  of  the  equation,  (sinh  B)/B  , 
presents  difficulty  when  trying  to  obtain  an  explicit  solution  for  concentration.  This  term 
accounts  for  the  difference  in  the  magnitude  of  the  sediment  attenuation  from  the  portion 
of  the  ensonified  volume  closest  to  the  transducer  to  the  portion  furthest  from  the 
transducer.  Figure  3-2  shows  the  magnitude  of  this  term  versus  the  term  B,  and  table  3-1 
shows  the  magnitude  of  this  term  for  several  cases  using  a  distribution  of  quartz  sand 
with  median  grain  diameter,  /it  =  2.65  and  standard  deviation,  o^  =0.25  .  Note  these 

parameters  are  given  in  units  of  phi,  defined  as  0  =  -log2  d  where  d  is  the  grain 
diameter  in  mm.  The  transducer  frequencies  listed  in  table  3-1  are  the  highest 
frequencies  used  in  this  study.  A  30  g/l  mass  concentration  corresponds  to  roughly  a  1% 
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Figure  3-2.  Significance  of  sinh(B)/B  term  in  acoustic  backscatter  equation. 


Table  3-1.  Values  of  sinh(B)/B  for  extreme  concentrations. 


Transducer 
frequency  (MHz) 

Concentration  (g/1) 

B 

sinhB 
B 

5 

30 

0.5591 

1.053 

5 

5 

0.0994 

1.002 

2.25 

30 

0.0956 

1.002 

2.25 

5 

0.0172 

1.000 
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concentration  by  volume,  the  approximate  upper  limit  of  concentration  before  multiple 
scattering  must  be  considered,  and  the  5  g/1  concentrations  are  the  highest  tested  in  this 
study.  From  these  results,  it  can  be  seen  that  for  the  highest  concentrations  for  which  the 
present  theory  applies,  this  term  can  be  significant  when  the  highest  frequency 
transducers  are  used.  When  operating  transducers  of  lower  frequency  with  suspensions 
of  lower  concentration,  this  term  is  very  close  to  unity.  Additionally,  since  the  sediment 
attenuation  term  within  the  exponential  function  is  an  empirically  determined  parameter, 
determination  of  its  value  without  including  the  (sinh  B)/B  term  may  compensate  for  its 
absence.  So,  for  the  remainder  of  this  text,  the  acoustic  backscatter  equation  will  be 
approximated  by  the  following  expression: 


VV=AF(z)C(z>^("-+"J 


(3-3) 


In  equation  (3-3),  the  system  sensitivity  constant,  S,  the  speed  of  sound,  c,  and  the  pulse 
width,  T,  have  been  combined  into  a  single  system  constant,  A.  The  system  constant  can 
be  later  seperated  back  into  these  constituents  in  order  to  correct  for  sound  speed 
variations. 

3.2  Development 
3.2.1   Concentration  inversion 

Beginning  with  the  general  form  of  the  acoustic  backscatter  equation  (3-4),  the 
concentration  dependence  is  removed  from  the  sediment  attenuation  term,  a  ,  giving  a 
form  in  which  the  multiple  term  dependence  on  concentration  is  more  obvious  (3-5). 

AF(z)c(z)  =  V2(z)z2exp(4z(a„,+a,(z))) 

(3-4) 
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Af(z)c(z)  =  V2(z)z2exp 


J4(a„,  +  £(z)c(z)>/z 


(3-5) 


Following  the  LH95  derivation,  first,  the  natural  logarithm  of  the  equation  is  found 

In  A  +  In  F  +  In  C  =  2 ln(Vz)  +  J(4a„,  +  AC,C)dz  (3_D) 


and  then  the  derivative,  denoted  by  (')■ 

'V'z+V 


F      C 


Vz 


+  4aH.  +  4fC 


(3-7) 


Upon  arranging  the  terms  of  equation  (3-7),  a  nonlinear  differential  equation  of  the 
Bernoulli  type  results. 


F'    2(^\-4a 


F  Vz 


C  =  4fC2 


(3-8) 


Rewriting  in  standard  form,  equation  (3-8)  becomes  the  following: 
C'+,,{z)C  =  Q{z)C" 


(3-9) 


with 
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F' 
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--At; 
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4a 

(3-10) 

(3-11) 
(3-12) 


Following  the  standard  method  of  solution  for  a  Bernoulli  equation,  the  following 
substitutions  can  be  made: 


t  =  C'~"=C';    c  =  r 


(3-13) 
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dC  _dC  dt 

dz      dt  dz 


(3-14) 


These  substitutions  result  in  a  readily  solved  first  order  linear  inhomogeneous  differential 
equation. 

(3-15) 

(3-16) 


-r2t'+  pt~l  =Qt" 
t'+pt  =  Q 


Where,  in  equation  (3-16),  the  following  apply: 

p  =  -p 
Q  =  -Q 


(3-17) 


Solution  to  equation  (3-16)  is  found  by  first  determining  the  integrating  factor. 
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M. 
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F 


(3-18) 
(3-19) 
(3-20) 


After  multiplying  the  equation  by  the  integrating  factor,  an  exact  differential  results, 
which  can  then  be  integrated  for  solution. 


d_ 
dz 

(Vz. 


- — — exp(4aaz)f 
F 


-4^    '  exp(4«„.z) 
F 


LrM 


exp(4a„z)r  =  -[4£- — — exp(4a„z)dz-l-y 
J         F 


(3-21) 
(3-22) 


Equation  (3-22)  is  then  solved  for t,  and  then  finally  for  the  concentration,  C. 


34 


M 


= a; Z_ (3-23) 

- — —  exp(4a„,z) 
F 

I  ^r-exp(4a,1z) 

'■-  =  — j — ^y (3-24) 

7-j4C^i-exp(4a„i)rfz 


Next,  the  boundary  conditions  (3-25)  are  applied  for  solution  of  the  integration  constant. 


C  =  C  I 

V  =  V0  \  at  z  =  z0 

F  =  Fa\ 


(3-25) 


7  =  ^f^-exp(4a,„z0)  (3_26) 


The  concentration  Co  at  the  nearfield  limit  z„ ,  which  is  the  closest  point  to  the  transducer 
not  in  the  nearfield,  may  be  estimated  in  many  applications  by  assuming  the 
concentration  and  size  are  constant  within  the  nearfield.  With  this  assumption,  the 
implicit  form  of  the  acoustic  backscatter  equation,  equation  (3-5),  simplifies  somewhat, 
giving  equation  (3-27). 

C0=%P*^>  (3-27) 

This  can  be  solved  iteratively  for  concentration  for  a  given  grain  size,  using  a  zero 
sediment  attenuation  form  of  the  equation  for  an  initial  estimate. 


C„ 


.  t^o)2  „.!«.. 


M- 


(3-28) 
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3.2.2  Size  determination 

The  technique  originally  introduced  by  Hay  and  Sheng  (1992)  for  determining  the 
median  size  of  particles  in  suspension  involved  first  approximating  the  acoustic 
backscatter  equation  (3-5)  by  neglecting  the  sediment  attenuation.  In  this  way  ratios 
could  be  constructed  from  the  approximate  equations  for  any  pair  of  unique  frequency 
transducers.  Since  the  concentration  dependence  of  the  attenuation  is  removed  by 
neglecting  the  sediment  attenuation,  the  remaining  concentration  terms  in  the  equation 
cancel  in  the  formation  of  the  ratio.  Hence,  the  only  remaining  unknowns  in  the  ratios 
are  functions  of  the  size  distribution  of  the  suspended  particles.  By  assuming  the  particle 
size  can  be  described  by  a  two-parameter  distribution,  the  log-normal  distribution,  and  by 
further  assuming  one  parameter  is  constant,  the  ratios  can  be  evaluated  over  a  range  of 
the  other  parameter.  The  standard  deviation  is  the  parameter  assumed  constant  and  the 
ratios  are  determined  over  a  range  of  median  particle  sizes.  Median  particle  size  is  then 
found  by  minimizing  the  difference  between  the  ratios  and  its  known  value  with  respect 
to  the  median  particle  size.  Crawford  and  Hay  (1993)  improved  the  technique  by  solving 
the  approximate  equations  first  for  those  terms  that  are  not  functions  of  the  size 
distribution  or  transducer  frequency.  These  terms  are  equal  in  all  of  the  equations, 
regardless  of  transducer  frequency,  so  equating  the  remaining  terms  in  the  other 
equations  eliminates  the  concentration  dependence.  Again,  the  minimization  technique  is 
applied  as  before  to  determine  the  median  particle  size. 

First  by  solving  each  transducer's  equation  only  for  concentration,  and  then  by 
minimizing  the  variance  in  the  concentrations  predicted  by  any  number  of  transducers 
with  respect  to  median  sediment  size,  Crawford  and  Hay's  technique  is  here  slightly 
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modified.  By  using  the  explicit  solution  for  concentration,  (3-24),  there  is  no  longer  the 
need  to  neglect  the  sediment  attenuation.  In  addition,  once  the  median  particle  size  is 
found  by  the  minimization  technique,  the  mean  concentration  is  readily  computed  from 
the  existing  concentration  solutions.  The  exact  procedure  is  as  follows. 

The  explicit  solution  for  concentration  derived  above  is  discretized  to  represent 
each  frequency  of  m  transducers  with  n+1  bins  outside  the  nearfield  region  of  each 
transducer. 

^^-exp(4o;,„z/)     .  (3-29) 

Fj  .  i  =  l,2,...,m 

yrl,j  ,/  =  l,2,...,« 

In  the  denominator  of  this  expression,  /  is  the  discrete  form  of  the  integral,  given  by  the 
following: 


,-4 


^,,%^exp(4«,,Zt)+C,JKt;Zt-|)2exp(4«^t.,) 

r.  ,  r,  ,_, 


The  integration  constant  is  defined  again  at  point  z.o- 

7,  =  <^exp(4a,Zo)  (3-31) 

It  should  be  mentioned  that  zero  concentration  in  suspension,  which  results  in  zero 
voltage  read  at  the  transducer,  will  result  in  an  indeterminate  value  of  the  integration 
constant.  The  location  of  z0  should  thus  be  set  at  the  first  range  with  non-zero 
concentration,  determined  by  the  following  iterative  technique. 
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The  initial  concentration  in  the  profile,  located  at  the  first  point  outside  all  of  the 
transducer  nearfields,  is  found  for  each  transducer  by  the  former  iterative  technique. 

Ci-M/**«  (3.32) 

Two  considerations  should  be  taken  when  determining  this  initial  concentration. 
First,  because  of  attenuation,  the  magnitude  of  the  voltage  read  from  the  transducer  is 
limited  in  magnitude  from  the  above  expressions,  as  is  apparent  in  examination  of 
figure  3-1.  Due  to  the  statistical  fluctuations  in  the  backscattered  signal  and  also  to 
instrument  noise,  it  is  possible  that  the  actual  signal  is  higher  in  magnitude  than  this 
theoretical  limit.  In  such  cases,  the  iterative  technique  will  not  converge  to  a  solution.  A 
simple  divergence  test  in  the  iteration  algorithm  will  reveal  this  condition.  Minimization 
of  the  difference  between  the  initial  concentration  guess  and  that  returned  by  equation 
(3-32)  will  produce  a  good  concentration  approximation  in  such  cases.  Second,  it  is 
important  to  realize  that  equation  (3-32)  will  produce  two  concentration  solutions  for  the 
reasons  discussed  previously.  For  the  lowest  frequency  transducers  typically  used,  the 
higher  magnitude  solution  is  regularly  above  the  expected  range  of  applicability  of  the 
present  theory,  and  can  safely  be  ignored.  For  the  higher  frequency  transducers,  the 
decision  of  which  concentration  to  use  must  be  based  on  physical  arguments  or  by 
comparison  with  the  results  from  lower  frequency  transducers.  In  field  measurement  of 
suspended  sediment,  the  transducer  is  usually  a  sufficient  distance  from  the  bed,  such  that 
the  higher  magnitude  solution  is  again  outside  the  expected  range  and  can  be  safely 
ignored.  For  example,  a  typical  5.0  MHz  transducer  with  a  16  cm  nearfield  will  give  two 
solutions  atz0=  16  cm.  As  seen  in  figure  3.1,  the  lower  magnitude  solution  will  fall 
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between  0  g/1  and  1 .6  g/1.  and  the  higher  above  1 .6  g/1.  Should  the  transducer  be  located 
a  half  meter  from  the  bed,  concentrations  above  1 .6  g/1  will  not  likely  be  found  at  zo  =  16 
cm  (34  cm  above  the  bed),  and  the  lower  magnitude  solution  is  the  most  probable 
concentration. 

The  complete  technique  for  solution  is  as  follows.  Equation  (3-32)  is  solved 
iteratively  for  the  initial  concentration  for  each  transducer  over  a  range  of  median  grain 
diameters.  Note  that  both  F  and  £  are  functions  of  the  grain  size  distribution,  typically 
assumed  to  be  log-normal.  Calculation  of  these  parameters  first  requires  one  to 
determine  the  distribution  based  on  the  given  median  grain  diameter.  The  form  of  the 
log-normal  distribution  is  given  in  the  appendix.  The  standard  deviation  of  the  grain  size 
distribution  is  assumed  constant  and  is  determined  by  other  physical  arguments.  In  the 
case  of  field  measurement  of  sediment  suspension,  the  standard  deviation  is  generally 
assumed  equal  to  that  of  the  distribution  of  sediment  in  the  seabed  below  the  transducer. 
For  each  median  grain  diameter,  the  mean  concentration  and  the  variance  in  the 
concentration  between  transducers  is  calculated. 

A*c=^|c,  (3-33) 

<Tr2=-ic,'-(lf;c,T  (3-34) 

The  median  grain  diameter  is  recognized  as  that  with  the  minimum  concentration 
variance,  and  the  concentration  is  given  by  the  corresponding  mean  concentration.  If 
only  two  transducers  are  used,  more  than  one  solution  for  the  median  grain  diameter  is 
possible.  In  this  case,  determination  of  size  is  still  possible  if  the  range  of  grain  sizes  is 
restricted  and  appropriate  transducer  frequencies  are  selected  in  advance. 
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After  the  initial  concentration  is  found,  the  integration  constant,  y,  can  be  found 
for  each  transducer  from  equation  (3-31).  The  solution  for  the  remainder  of  the  bins  in 
the  profile  proceeds  by  solving  equation  (3-29)  for  each  transducer  for  a  range  of  median 
grain  diameters,  and  then  by  selecting  the  correct  grain  diameter  by  minimization  of  the 
concentration  variance  between  transducers,  given  by  equation  (3-34).  Again,  the 
concentration  is  given  by  the  corresponding  mean  value  from  all  transducers. 


3.3  Verification  of  Technique 
As  an  initial  test  of  the  calibration  technique,  ideal  voltage  profiles  for  1 .0,  2.25, 
and  5.0  MHz  transducers  were  simulated  using  the  acoustic  backscatter  equation  (3-4) 
(a)  (b) 
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Figure  3-3.  Numerically  generated  profiles  for  (a)  1 .0,  (b)  2.25,  and  (c)  5.0  MHz 
transducers. 
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with  a  predetermined  log-normal  distribution.  For  simplicity,  the  generated  profiles 
contained  no  nearfield,  or  in  other  words,  the  acoustic  backscatter  equation  is  assumed  to 
be  valid  at  the  face  of  the  transducer  and  beyond.  Figure  3-3  shows  these  numerically 
generated  profiles.  The  concentrations  used  were  0.01,  0.02,  0.04,  0.08,  0.16,  0.32, 
0.63.1 .25,  2.5,  and  5.0  g/1,  and  the  grain  size  distribution  was  assumed  to  have  a  median 
grain  diameter  fi^  =  2.66  and  a  standard  deviation  c4  =  0.25.  Since  in  the  absence  of 

attenuation,  the  voltage  read  from  the  transducer  increases  with  increasing  suspended 
particle  concentration,  the  concentrations  are  easily  distinguishable  at  zero  distance. 
Higher  attenuation  with  higher  operational  frequency  and  with  higher  concentration  is 
apparent  from  the  figures.  Particularly,  the  profiles  from  the  highest  concentrations  in  the 
5  MHz  simulation  are  attenuated  so  heavily  that  away  from  the  transducer  little,  if  any, 
signal  remains, 
a)  b)        


Figure  3-4.  Resulting  (a)  concentration  and  (b)  median  grain  size  from  inversion 
technique  using  2.25  and  5.0  MHz  profiles.  Asterisks  indicate  known  values. 


Shown  in  figure  3-4  are  the  concentration  and  size  profiles  resulting  from 
applying  the  new  technique  using  double  precision  calculations  (64  bit)  with  only  the 


41 


2.25  and  5.0  MHz  voltage  profiles  from  figure  3-3.  Since  only  two  frequencies  were 
used  in  this  test,  the  range  of  median  grain  sizes  was  restricted  within  +/-  2  standard 
deviations  of  the  initial  known  distribution.  In  figure  3-4b,  the  size  profiles  overlap  for 
all  cases  except  the  cases  involving  the  highest  two  concentrations.  Not  surprisingly,  the 
technique  accurately  produced  the  initial  concentration,  as  in  Holdaway  and  Thome's 
(1997)  simulations,  and  size  in  the  majority  of  the  cases.  This  test  does  illustrate  the 
(a)  (b) 


Figure  3-5.  Resulting  (a)  concentrations  and  (b)  median  grain  size  from  inversion 
technique  using  1 .0  and  2.25  MHz  profiles.  Asterisks  represent  known  values. 


(a) 


Figure  3-6.  Resulting  (a)  concentrations  and  (b)  median  grain  size  from  inversion 
technique  using  1.0,  2.25,  and  5.0  MHz  profiles.  Asterisks  represent  known  values. 
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difficulty  though  in  using  the  highest  frequency  transducer  through  significant  ranges  of 
high  concentrations.  As  noted  before  in  figure  3-3,  at  the  highest  concentrations,  the  high 
sediment  attenuation  results  in  large  signal  loss  away  from  the  transducer.  In  figure  3-4, 
the  error  induced  in  evaluating  the  concentration  and  size  from  this  small  signal  is 
apparent.  Figure  3-5  shows  the  results  from  use  of  the  1 .0  and  2.25  MHz  signals.  The 
results  in  this  test  showed  excellent  agreement  with  both  the  original  size  and 
concentrations,  even  through  a  one-meter  range  of  relatively  high  concentration.  The 
final  inversion,  shown  in  figure  3-6,  uses  all  three  of  the  simulated  signals.  Again, 
because  of  the  signal  loss  in  the  high  concentration  5  MHz  data,  the  results  exhibit  similar 
behavior  to  those  in  figure  3-4.  For  this  reason,  it  is  important  to  be  aware  of  signal  loss 
when  working  with  the  highest  frequency  transducers,  such  that  this  can  be  considered  in 
the  inversion  algorithm. 

Laboratory  tests  were  performed  in  a  recirculating  calibration  chamber  which 
produces  a  suspension  of  sediment  of  approximately  constant  concentration  and  constant 
grain  size  distribution  (Hanes  et  al.,  1988).  Transducer  frequencies  of  1.0,  2.2,  and  5.0 
MHz  were  used  for  the  measurements.  Backscattered  intensity  signals  were  collected  at 
100  Hz  with  each  transducer,  and  the  ensemble  RMS  was  determined  from  1  minute's 
data  for  each  concentration.  As  in  the  numerical  simulations,  the  distribution  parameters 
for  the  sediment  were  a  median  grain  diameter  fit  =  2.66  and  a  standard  deviation  a^  = 

0.25.  The  concentrations  used  were  0,  0.1,  0.2,  0.3,  0.4,  and  0.5  g/1.  Figure  3-7  shows 
the  agreement  between  the  known  concentration  and  that  determined  with  the  new 
inversion  technique.  Several  factors  explain  the  form  of  the  error  curve  with  range.  First, 
concentrations  in  the  chamber  were  determined  by  adding  calculated  dry  masses  of 
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Figure  3-7.  Comparison  of  known  and  calculated  concentrations  (a)  shown  at  distinct 
concentrations  and  ranges  of  (*)  40  cm,  (o)  50  cm,  (x)  60  cm,  and  (+)  70  cm;  and  (b) 
shown  as  the  mean  error  of  all  concentrations  with  range. 


Figure  3-8.  (a)  Comparison  of  known  (broken  line)  and  determined  (solid  line)  median 
grain  size  and  (b)  the  resulting  error  with  range. 


sediment  to  the  known  volume  of  water.  Due  to  hindered  settling  within  the  funnel  at  the 
base  of  the  chamber,  actual  concentrations  in  the  tube  may  be  slightly  lower  than  those 
calculated.  Next,  the  initial  concentration  measurement  is  located  40  cm  from  the 
transducer  face  due  to  both  nearfield  effects  from  the  transducers  and  from  complications 
introduced  by  amplifier  saturation  at  shorter  ranges.  Determination  of  the  initial 
concentration  must  therefore  be  done  with  a  signal  that  has  already  experienced  40  cm  of 
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water  and  sediment  attenuation  through  significant  concentrations.  The  sensitivity  of  the 
concentration  measurement  to  attenuation  increases  the  likelihood  of  error,  particularly 
for  high  concentrations,  long  attenuation  paths,  and  high  operational  frequencies. 
Finally,  for  calibration  of  the  acoustic  transducers,  the  error  between  known  and 
calculated  concentrations  was  minimized  in  the  range  from  40  cm  to  90  cm.  In  this 
minimization  technique,  approximately  half  of  the  calculated  concentration  profile 
typically  falls  below  the  known  value  and  half  above.  This  effect  is  apparent  in  the  error 
profile  of  the  present  concentration  evaluation,  Figure  3-7b,  which  shows  the  best 
agreement  in  the  center  portion  of  the  profile.  If  the  transducer  calibration  were 
performed  at  just  a  single  range,  as  is  often  described  in  the  literature,  and  the 
measurements  presented  in  this  text  evaluated  using  the  single  point  calibration 
information,  the  error  in  the  determination  of  concentration  is  less  than  5%. 

Figure  3-8  compares  the  known  median  grain  size  with  that  measured  in  the 
circulation  chamber.  Again,  if  the  calibration  is  performed  at  a  single  range  and  the  size 
determined  from  the  single  point  calibration,  the  error  in  evaluation  of  median  grain  size 
is  less  than  10%.  Even  in  this  case,  the  evaluated  median  grain  size  is  slightly  higher 
than  the  known  value.  Errors  result  in  this  evaluation  from  use  of  somewhat  low 
concentrations  for  determination  and  from  differences  in  various  sediments  not  accounted 
for  explicitly  in  the  empirical  form  function.  Use  of  low  concentrations  was  made 
necessary  by  operation  of  the  5.0  MHz  transducer  with  an  initial  concentration  evaluation 
point  located  a  significant  distance  from  the  transducer.  Presently,  one  form  function  is 
said  to  describe  noncohesive  quartz  sediment  (see  appendix),  but  it  is  expected  that  grain 
properties  of  a  given  sediment  sample  will  modify  the  form  function  slightly.  Empirical 
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evaluation  of  the  form  function  for  a  given  sediment  type  would  likely  improve  the  error 
in  determination  of  the  median  grain  size. 

It  should  be  noted  that  evaluating  concentration  and  size  with  a  constant 
concentration  profile  is  actually  a  more  demanding  application  of  the  technique  and 
system  than  is  typically  experienced  in  field  measurements,  due  to  the  propagation  of 
error  through  the  profile.  In  measurements  of  sediment  suspension  above  the  seabed,  the 
transducer  is  typically  far  enough  from  the  seabed  such  that  the  concentrations  near  the 
transducer  are  low. 

3.4  Summary 
A  new  technique  of  determining  both  concentration  and  the  median  grain 
diameter  of  suspended  particles  has  been  presented.  The  significant  advantage  of  the 
technique  is  that  by  using  an  explicit  solution  for  concentration,  the  median  grain 
diameter  can  be  found  without  having  to  neglect  sediment  attenuation.  In  addition, 
because  incorporating  the  correct  median  grain  diameter  in  the  explicit  solution  will 
produce  an  identical  concentration  regardless  of  the  operational  frequency  of  the 
transducer,  the  concentration  variance  between  any  number  of  transducers  can  be 
minimized  to  find  the  median  grain  diameter.  Numerical  simulations  show  the  technique 
produces  both  the  expected  concentrations  and  grain  diameters.  In  addition,  laboratory 
results  from  a  recirculating  calibration  chamber  verify  that  the  technique  applies  well  in 
determining  sediment  size  and  concentration  from  measurements  of  backscattered 
acoustic  intensity. 


CHAPTER  4 
SUSPENSION  TIME  SCALES 

In  the  previous  chapters,  a  new  system  of  instrumentation  capable  of  accurate 
depiction  of  sediment  suspension  processes  with  high  spatial  and  temporal  resolution  has 
been  described.  In  addition,  a  new,  robust  process  of  data  conversion  from 
multifrequency  acoustic  backscatter  data  to  concentration  and  median  sediment  size  has 
been  introduced.  In  the  fall  of  1996,  a  system  similar  to  that  described  in  the  first  chapter 
was  deployed  from  the  Sensor  Insertion  System  (SIS)  at  the  Army  Corps  of  Engineers 
Field  Research  Facility  (FRF)  in  Duck,  North  Carolina.  In  this  chapter,  the  time  series  of 
concentration  profiles  obtained  from  the  acoustic  backscatter  measurements  collected  at 
this  project  are  examined  in  relation  to  the  instantaneous  hydrodynamic  measurements. 

It  is  common  in  the  study  of  sediment  suspension  in  a  wave  environment  to 
decompose  the  concentration  into  steady  and  fluctuating  components  (Nielsen,  1992). 
The  significance  of  each  component  in  a  sediment  transport  calculation  depends  on  the 
relative  importance  of  the  two  transport  mechanisms:  transport  by  currents  or  transport  by 
waves.  In  the  longshore  direction,  the  sediment  flux  computed  from  the  product  of  the 
mean  concentration  and  steady  current  velocity  has  been  used  successfully  for 
determining  the  mean  longshore  rate  of  transport  (Hanes  and  Huntley,  1986).   In  general, 
determination  of  the  mean  vertical  suspended  sediment  concentration  profile  typically 
involves  computation  of  the  near  bed  concentration  by  use  of  a  reference  concentration 
model  and  computation  of  concentrations  above  by  a  vertical  distribution  model.  Of  the 
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many  models  of  reference  concentration,  a  simple  linear  relation  between  the  bed  shear 
stress  and  reference  concentration  is  shown  to  work  best  (Smith  and  McLean,  1977; 
Thosteson,  1995).  The  mean  vertical  concentration  distribution  is  best  described  by  a 
model  incoiporating  both  turbulent  diffusion  and  vertical  convection  due  to  vortex  ripples 
(Nielsen,  1992;  Lee,  1994).  In  examination  of  cross-shore  transport,  the  fluctuating 
component  of  concentration  becomes  more  important  (Huntley  and  Hanes,  1987). 

There  is  significant  evidence  of  the  importance  of  low  frequency  water  wave 
motion  in  the  process  of  cross-shore  sediment  transport  (Huntley  and  Hanes,  1987;  Beach 
and  Sternberg,  1991;  Osborne  and  Greenwood,  1992).  This  has  been  attributed  to  the 
fixed  phase  difference  between  the  components  of  velocity  and  concentration,  where  the 
low  frequency  velocities  result  from  free  and  group-bound  infragravity  waves.  Low 
frequency  variation  in  concentration  has  been  shown  to  be  associated  with  wave  groups 
(Hanes,  1991).  The  aim  of  the  present  investigation  is  to  examine  the  significance  of  the 
concentration  variation  at  various  frequencies  and  to  further  examine  suspension  by  wave 
groups.  It  is  hoped  that  this  will  aid  in  development  of  future  models  that  predict  the 
fluctuating  components  of  concentration. 

4.1  The  SIS96  Project 
Shown  in  figure  4-1  is  the  SIS  on  the  FRF's  pier.  All  instrumentation  is  deployed 
from  the  SIS,  which  consists  of  a  crane  mounted  on  tracks  that  extend  along  the  length  of 
the  pier.  Instruments  attached  to  an  arm  (called  the  "bah")  at  the  end  of  the  crane's  boom 
can  be  positioned  with  reasonable  precision  at  locations  near  the  seabed.  Four  bayonets 
located  at  the  end  of  arm  closest  to  the  pier  are  forced  into  the  seabed  by  the  weight  of  the 
crane  to  stabilize  the  bah.  Cross-shore  movement  of  the  SIS  allows  data  collection  to  be 
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Figure  4-1.  Sensor  Insertion  System  (SIS) 


performed  in  areas  with  varying  sediment  composition  and  wave  conditions.  The 
bulkiness  of  this  arrangement  immediately  suggests  that  the  framework  will  interfere 
with  the  processes  to  be  measured.  While  there  is  certainly  an  effect  from  the  presence  of 
the  structure,  precautionary  measures  are  taken  to  minimize  this  impact.  The  arm  itself 
consists  of  pipe  of  smaller  diameter  than  the  main  structure,  and  is  just  massive  enough 
not  to  flex  by  wave  forcing.  Next,  the  arm  is  distanced  as  far  from  the  bed  as  possible 
considering  the  range  limitations  of  the  acoustic  instruments.  Orientation  of  the  arm  is 
longshore,  such  that  it  will  have  the  minimum  influence  in  the  vertical  and  cross-shore 
directions.  In  the  presence  of  long-shore  currents,  measurements  are  taken  on  the 
appropriate  side  of  the  pier  to  be  upstream  of  the  pier  and  the  main  structure  of  the  crane. 
Finally,  the  instruments  themselves  are  located  at  the  end  of  the  arm  farthest  from  the 


49 

main  structure  and  pier.  Though  these  measures  minimize  the  influence  of  the  supporting 
structure  and  the  pier  itself,  it  is  expected  that  some  influences  endure.  This  is  accepted 
as  a  trade-off  for  the  ease  of  repositioning  the  test  site  and  reconfiguring  instrumentation. 
It  should  be  mentioned  that  even  after  the  longest  deployments  of  the  instrumentation 
(approximately  1  hour),  no  scour  hole  could  be  seen  with  the  MTA  measurements, 
suggesting  that  the  influence  of  the  arm  itself  was  minimal. 

As  mentioned  previously,  the  system  of  instrumentation  used  in  this  project  was 
slightly  different  from  that  described  in  chapter  one.  In  fact,  this  project  inspired  many  of 
the  new  features  of  the  newer  system.  The  fundamental  limitation  of  the  system 
deployed  in  this  project  was  the  constraint  on  the  duration  of  data  collection.  As  will  be 
recognized  later,  the  analysis  of  the  data  from  these  experiments  indicates  that  long 
records  of  suspension  must  be  examined  to  capture  low  frequency  events  that  can 
dominate  the  record.  The  other  difference  in  the  systems  is  in  the  instruments  utilized. 
Figure  4-2  shows  the  instruments  as  they  were  positioned  on  the  arm  of  the  SIS  during 
the  project.  Of  the  instruments  shown  in  figure  4-2,  the  following  instruments  were  a 
part  of  the  system  and  used  in  this  investigation:  2  Simrad  Mesotech  model  810  ACPs,  1 
Sontek  ADV,  1  TransMetrics  P21LA-25  PSIS  pressure  sensor,  3  Seatek  MTAs,  1 
DeepSea  MicroSeaCam  1050,  and  1  D&A  OBS  3.  The  remaining  instruments  shown  in 
the  figure  were  owned  and  operated  by  the  FRF.  Note  that  two  individual  Simrad 
Mesotech  ACPs  with  frequencies  of  2.25  and  5.0  MHz  were  used  instead  of  the  ABS 
system  described  in  the  previous  chapters.  As  will  be  seen  later,  this  introduced 
complication  into  trying  to  determine  size  from  the  backscatter  data. 


50 


SIS96  Instrument  Arrangement 


View  from  South 
when  on  pier 


Figure  4-2.  Instrumentation  used  in  project. 


Table  4-1.  Calibration  constants  for  instruments  used  in  present  analysis. 


Instrument 

Gain 

Offset 

System  constant 

Pressure  sensor 

8.35x10"'  W 

/count 

(salt  water) 

-6.77  m 
(salt  water) 

- 

Current  X 

1.187xl0-1n/ 

/s  count 

-2.457  m/ 

- 

Current  Y 

1.203x10  '  W 

/s- count 

-  2.492  W 
/  s 

- 

OBS 

1.409x10-'^/ 

/I  count 

-0.164  y( 

- 

2.25  MHz  ACP 

- 

11  mV 

0.464 

5.00  MHz  ACP 

- 

2mV 

0.929 
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Calibration  of  the  instruments  which  measure  from  only  a  single  point  in  space, 
termed  single  point  instruments  (SPIs),  is  straight  forward.  Hydrostatic  measurements 
are  used  for  calibration  of  the  pressure  sensor,  and  the  ADVs  are  calibrated  using  a 
moving  cart  in  a  tank  of  still  water.  The  OBS  sensor  is  calibrated  in  a  sediment- 
recirculating  calibration  tank  with  the  ACPs.  Table  4- 1  shows  the  calibration  constants 
for  the  various  instruments.  As  described  in  the  previous  chapter,  the  only  undetermined 
parameter  for  the  ACPs  is  the  system  constant.  Due  to  a  slight  DC  offset  in  the  output 
signals  from  the  ACPs,  table  4-1  also  lists  the  optimum  offset  for  each  transducer.  An 
optimization  process  is  utilized  to  determine  both  the  system  constant  and  DC  offset 
which  produces  the  concentration  profiles  closest  to  the  known  concentrations  in  the 
sediment-recirculating  tank.    The  predicted  concentration  profiles  versus  the  known 
values  of  concentration  for  the  two  ACPs  are  shown  in  figure  4-3.  Across  all  locations 
and  concentrations,  the  mean  errors  for  the  2.25  MHz  and  5.0  MHz  transducers  are  8.8% 
and  17.9%  respectively. 


(a) 


2.25  MHz  S=0.464   DC=0.01lvi 


0.2  0.3 


0.7  0.8 


Figure  4-3.  Calculated  (solid  lines)  concentration  profiles  from  (a)  2.25  MHz  and  (b) 
5.00  MHz  calibration  data  using  optimum  system  constant  and  DC  offset  versus  known 
concentration  (dashed  lines). 
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For  four  days,  from  October  29  until  November  1,  1996,  experiments  were 
performed  from  various  locations  on  the  pier,  ranging  from  water  depths  of  1 .4  to  7.0 
meters.  Measurements  were  attained  with  wave  conditions  ranging  from  //„,„  =  0.35  to 
1 .0  m,  due  to  the  capability  of  moving  the  SIS.  Waves  were  a  mixture  of  locally 
generated  seas  with  peak  periods  near  6  seconds  and  an  underlying  swell  component  near 
1 1  seconds.  The  local  component  was  most  prevalent  at  the  start  of  the  week,  with  the 
swell  component  becoming  more  dominant  in  the  latter  days.  Grab  samples  of  sediment 
were  collected  at  the  start  of  each  data  collection  run,  and  sieve  analysis  indicated  the 
median  grain  sizes  ranged  from  120  to  200  microns  (3.06  and  2.32  respectively  on  the  0 
scale).  Sieve  analysis  also  indicated  that  the  grain  size  distributions  of  the  samples  from 
the  runs  investigated  in  the  present  study  (identified  in  section  4.3)  were  described  well 
by  the  log-normal  distribution.  Shown  in  figure  4-4  is  the  best  fit  normal  cumulative 
distribution  function  to  the  sieve  data  from  run  number  1 8.  Figure  4-5  shows  the 


Mean  phi  =  2.242  Std.  dev.  =  0.4482 


Figure  4-4.  Best  fit  normal  cumulative  distribution  function  to  sieved  grain  size  data; 
sample  taken  at  location  of  and  prior  to  run  18. 
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Figure  4-5.  (a)  Profile  taken  from  pier  on  north  (solid)  and  south  (dashed)  side  with 
superimposed  experimental  water  depths  (*)  and  (b)  corresponding  sieved  median  grain 
size. 


variation  of  the  median  sediment  grain  size  with  pier  location  and  the  beach  profile  as 
taken  over  the  edge  of  the  pier  with  a  plumb  bob.  Elevation  in  this  figure  is  relative  to 
the  mean  water  level  at  the  time  of  the  survey.  Jette  presents  additional  information 
regarding  the  SIS96  project  (1997). 


4.2  Size  Determination 
As  described  in  the  previous  chapter,  the  median  grain  diameter  of  the  particles  in 
suspension  can  be  determined  by  computing  the  concentrations  from  each  frequency 
transducer  over  a  range  of  sizes.  Again,  the  size  is  determined  to  be  that  which 
minimizes  the  variance  in  the  computed  concentrations  across  transducers. 
Determination  of  size  from  field  data  is  complicated  by  several  factors  that  sometimes  act 
in  conjunction.  These  complicating  factors  include  the  following:  multiple  zeros  in  the 
variance  versus  size  relation,  poor  measurement  resolution,  statistical  fluctuation  in  the 
backscattered  signal  returns,  and  spatially  separated  acoustic  beams. 
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Because  of  the  nonlinearity  of  the  form  function,  it  is  possible  that  the  variance  is 
zero  at  more  than  one  grain  size.  Increasing  the  number  of  unique  sound  frequencies 
decreases  the  likelihood  of  the  variance  having  multiple  roots.  In  this  project,  only  two 
transducer  frequencies  were  used,  making  the  determination  of  the  correct  root  difficult  at 
times.  It  is  not  uncommon  to  obtain  a  variance  versus  size  relation  as  that  shown  in 
figure  4-6.  In  such  a  case,  a  physical  argument  based  on  the  sizes  that  a  bed  sample 
contains  is  used  to  choose  the  more  likely  sediment  size.  If  the  backscattered  signal  is 
mainly  from  a  washload  component  in  the  suspended  sediment  (small  suspended  particles 
x10"8 


80  100         120         140         160         180         200         220         240         260         280 

Median  grain  diameter  (microns) 

Figure  4-6.  Multiple  roots  in  variance-size  relation. 
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not  found  locally,  but  instead  advected  from  another  region),  than  certainly  this  choice 
will  be  in  error. 

Next,  resolution  of  the  transducers  and  sampling  resolution  must  be  considered. 
For  instance,  in  the  present  experiments,  the  transducers  are  sampled  with  a  12-bit  analog 
to  digital  converter.  At  small  amplitudes,  a  one-count  change  in  the  measured 
backscattered  signal  can  result  in  a  change  of  order  in  determined  concentration.  Since 
ultimately  it  is  the  difference  in  concentration  measurements  across  transducers  which 
determines  the  median  sediment  size,  this  large  change  in  concentration  drastically 
changes  the  resulting  size  evaluation.  For  this  reason,  a  minimum  value  of  the 
backscattered  signal  strength  is  required  before  the  size  evaluation  can  be  trusted.  In 
practice,  only  sizes  obtained  from  concentrations  greater  than  50  mg/1  are  used. 

Because  the  measurement  of  the  returned  signal  from  suspended  sediment  is  a 
random  process,  many  instances  or  profiles  must  be  collected  to  obtain  statistically 
meaningful  results.  As  described  in  chapter  2,  the  root  mean  square  (RMS)  of  a 
predetermined  number  of  profiles  is  generated  by  the  acquisition  system.  The  number  of 
profiles  included  in  this  RMS  profile  is  chosen  such  that  the  error  in  the  concentration 
measurements,  proportional  to  the  reciprocal  of  the  square  root  of  the  number  of  profiles, 
is  minimal.  Again,  because  of  the  sensitivity  of  the  size  evaluation  to  small  concentration 
differences,  a  small  error  in  concentration  may  result  in  a  relatively  large  error  in  the 
evaluation  of  size.  For  this  reason,  an  RMS  computed  from  a  significant  number  of 
measured  signals  (averaging  of  the  measured  intensities)  must  be  used  in  order  to  obtain  a 
reliable  size  evaluation.  It  is  essential  that  the  sea  bottom  location  in  the  profile  be 
determined  prior  to  the  RMS  process,  since  movement  of  the  bottom  with  time  will  make 
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it  difficult  to  determine  its  location  from  an  RMS  profile.  For  this  same  reason,  with 
regard  to  bed  location,  the  closest  reliable  size  estimate  is  that  just  above  the  highest  bed 
location  over  the  averaging  period.  In  the  present  experiment,  there  was  significant 
movement  of  the  bed  location  in  the  profile,  as  indicated  in  figure  4-7.  The  change  in  bed 
location  was  due  not  to  accretion,  but  instead  to  settling  of  the  framework  into  the  seabed. 
This  was  verified  by  examination  of  the  MTA  bottom  profile,  which  showed  a  uniform 
movement  of  the  profile  with  little  change  in  the  ripple  field. 


Bottom  location  relative  to  that  at  start 


0  200         400         600         800        1000       1200       1400       1600       1800       2000 

Elapsed  time  (s) 
Figure  4-7.  Change  in  bottom  location  due  to  settling  of  instrument  framework. 
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Finally,  it  is  often  difficult  in  the  field  to  obtain  acoustic  backscatter 
measurements  from  transducers  with  collocated  beams.  Physical  restrictions  due  to  the 
size  of  the  transducers  require  that  the  beams  be  spatially  separated.  Since  a  significant 
amount  of  temporal  averaging  is  performed  before  the  analysis,  this  introduces  no 
difficulty  should  the  seabed  be  flat.  If  ripples  exist,  as  they  did  in  nearly  every 
experiment  in  this  project,  then  it  is  likely  that  one  transducer  will  receive  backscatter 
from  a  concentration  profile  which  is  higher  or  lower  spatially.  This  spatial  separation 
corresponds  to  the  transducers  being  positioned  above  different  areas  of  the  ripple 
wavelength.  In  such  a  case,  the  differences  in  concentrations  between  the  two 
transducers  results  not  only  because  of  sediment  size  variations,  but  also  because  the 
measurements  are  located  at  different  elevations  from  the  bed.  Away  from  the  bed,  the 
resulting  error  in  size  is  probably  small,  as  advection  and  mixing  will  remove  horizontal 
gradients  in  concentration.  However,  near  the  bed,  the  comparison  of  concentrations 
measured  from  different  elevations  from  the  bed  will  result  in  notable  error  in  the  size 
evaluation.  The  backscatter  profiles  can  be  aligned  based  on  bed  elevation  as  opposed  to 
alignment  by  transducer  face  locations,  but  this  will  result  in  the  neglect  of  any 
differences  in  the  true  vertical  concentration  and  size  profiles  over  a  ripple  wavelength. 

Despite  careful  consideration  of  each  of  the  above  complexities,  a  consistent 
estimate  of  the  median  sediment  size  could  not  be  obtained  from  the  data  collected  in  this 
project.  In  figure  4-8,  the  resulting  size  profiles  from  two  separate  data  runs  are  shown 
together  with  completely  different  results.  The  data  for  both  runs  were  collected  at  the 
same  location,  with  one  hour  between  the  start  times  of  the  runs.  Both  size  profiles  result 
from  RMS  profiles  of  the  entire  16  minute  runs.  Inspection  of  the  time  series  of  the  bed 
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(b) 
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Figure  4-8.  Perceived  median  grain  size  profiles  from  (a)  run  19  and  (b)  run  20  and 
corresponding  near  bed  concentration  profiles. 


elevations  from  each  transducer  indicates  the  bottom  location  was  the  same  for  each 
transducer.  A  sediment  sample  taken  from  the  bed  indicates  a  median  grain  diameter  of 
200  microns  was  present  at  this  location.  Although  neither  of  the  results  seems 
unreasonable,  examination  of  the  concentration  variance  across  transducers  indicates 
agreement  in  concentration  is  never  achieved.  Over  the  entire  range  of  the  trial  grain 
sizes,  selected  to  be  from  plus  or  minus  three  standard  deviations  of  the  local  median 
grain  size,  no  applied  size  will  result  in  equal  concentration  readings  across  transducers. 
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As  to  which  transducer  yields  the  higher  concentration,  there  is  no  consistency.  Either 
transducer  is  just  as  likely  to  respond  with  a  higher  concentration  measurement. 

Recall  that  the  beams  of  the  two  transducers  used  in  this  project  were  close 
together  (within  4  inches  of  one  another)  but  not  collocated.  For  this  reason,  it  is 
possible,  and  the  results  seem  to  indicate  that  a  horizontal  gradient  in  concentration  exists 
between  the  transducers.  Although  this  result  seems  unlikely  with  time  averaged  data, 
the  likelihood  of  occurrence  increases  if  the  time  series  of  concentration  is  dominated  by 
a  few  infrequent  suspension  events  -  a  hypothesis  which  is  verified  by  the  following 
analysis.  Furthermore,  this  result  indicates  the  importance  of  having  truly  collocated 
beams  if  an  evaluation  of  median  grain  size  is  to  be  performed. 

4.3  Dominant  Frequency  Band  of  Suspension  Events 
In  order  to  study  the  most  dominant  time  scales  of  sediment  suspension,  four  runs 
of  the  30  runs  collected  will  be  examined.  The  investigation  is  limited  to  these  four  runs, 
because  these  were  the  only  runs  of  sufficient  length  to  provide  confidence  in  the  low 
frequency  portions  of  the  spectral  analysis  to  be  presented.  Recall  that  in  the  system  used 
for  this  project  the  duration  of  data  collection  is  limited  by  the  available  memory  in  the 
data  logger,  unlike  the  newer  system  described  in  chapter  1 .  Table  4-2  shows  the 
conditions  under  which  the  experiments  were  performed.  In  this  table,  the  wave  height  is 
determined  by  correcting  the  pressure  time  series  for  depth  attenuation  using  linear  wave 
theory  and  is  then  verified  using  that  obtained  by  correcting  the  velocity  time  series, 
again  using  linear  wave  theory  (Dean  and  Dalrymple,  1984).  In  this  calculation, 
contributions  from  wave  periods  less  than  3  seconds  and  greater  than  20  seconds  are 
removed,  albeit  examination  of  the  spectra  prior  to  removal  shows  little  energy  in  these 


Table  4-2 
study. 
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Conditions  at  the  measurement  site  during  experiments  examined  in  this 


Run 

Date 

Time 
(EDT) 

Duration 
(MM:SS) 

Pier  side 

Location 
(m) 

djo 
(urn) 

Depth 
(m) 

18 

10/31/97 

12:39:36 

32:14 

South 

207 

211 

2.41 

23 

11/01/97 

09:01:08 

43:44 

South 

226 

194 

3.09 

25 

1 1/01/97 

10:57:54 

37:56 

North 

226 

187 

3.20 

27 

11/01/97 

12:35:24 

37:56 

North 

238 

194 

3.60 

Run 

HmO 

(m) 

(») 

Opeak 

o 

(cm/s) 

u  direction 

o 

l\  ripple 

(cm) 

(cm) 

18 

0.53 

10.0 

97 

16 

350 

0.5 

13 

23 

0.57 

10.9 

88 

10 

11 

0.5 

15 

25 

0.39 

11.7 

88 

9 

3 

1.5 

14 

27 

0.51 

10.1 

88 

5 

5 

1.2 

11 

portions.  The  angles  given  are  in  a  reference  frame  so  0°  is  directed  longshore  to  the  left 
of  an  onshore  observer  looking  out  at  sea,  and  they  increase  in  the  clockwise  direction. 
Wave  direction  is  found  using  the  maximum  entropy  method  with  the  pressure  and 
horizontal  velocity  data,  and  is  given  as  the  direction  from  which  the  waves  propagate. 
In  all  cases  listed,  the  instruments  were  placed  offshore  of  the  breakpoint  such  that  the 
instruments  remained  submerged  for  the  duration  of  each  run.  In  addition,  collection 
away  from  breaking  waves  avoids  contamination  of  the  concentration  measurements  by 
bubbles.  Due  to  the  difficulties  mentioned  above  in  evaluating  the  grain  size  profiles,  the 
assumption  of  constant  grain  size  is  utilized  in  the  conversion  from  the  backscattered 
signal  to  concentration.  The  grain  size  used  in  the  analysis  is  that  obtained  from  a  grab 
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Figure  4-9.  Surface  elevation  spectrum  for  run  18  with  80%  confidence  intervals. 


0.25 


sample  taken  in  the  vicinity  of  the  instruments  prior  to  data  collection.  Figure  4-9  shows 
the  surface  elevation  spectrum  from  run  number  18.  In  this  case  and  in  each  case  that 
follows,  the  spectrum  is  found  using  the  entire  record  of  surface  elevation.  A  Bartlett 
spectral  window  is  then  used  to  smooth  the  spectrum,  resulting  in  a  spectral  estimate  with 
approximately  6  degrees  of  freedom.  Note  that  the  most  significant  portion  of  the  energy 
is  found  in  the  incident  band  of  the  surface  elevation  spectrum.  Also,  note  the  lack  of 
energy  in  the  lowest  frequency  band  -  that  at  the  extreme  left  in  the  figure.  Logged 
observations  indicate  that  the  waves  consisted  of  swell  from  an  offshore  low-pressure 
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Figure  4-10.  Plot  of  the  100  mg/1  contour  for  whole  time  series  of  run  1  i 
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Figure  4- 1 1 .  Plot  of  the  1  g/1  contour  for  whole  time  series  of  run  1 8. 
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system  with  little  or  no  locally  generated  wind  waves,  in  agreement  with  the  measured 
spectrum. 

A  qualitative  observation  of  the  concentration  profile  time  series  indicates  that 
suspension  events  in  general  tend  to  be  intermittent  with  only  infrequent  occurrences  of 
high-concentration  events,  as  described  by  Hanes  (1988).  This  can  be  seen  in  figure  4- 
10,  where  the  100  mg/l  contour  as  determined  from  the  2.25  MHz  transducer  is  plotted. 
For  comparison,  the  same  time  series  is  shown  in  figure  4-11,  but  the  one  gram  per  liter 
contour  is  instead  plotted.  Inspection  of  figure  4-1 1  shows  that  the  high  concentration 
events  are  indeed  less  frequent  than  the  lower  concentration  events  seen  in  figure  4-10. 
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Figure  4-12.  Time  series  of  the  squared  bottom  velocity  magnitude  and  vertically 
integrated  concentration. 
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In  addition,  these  infrequent  high  concentration  events  stay  confined  to  the  region  very 
close  to  the  bed.  Although  these  high  concentration  events  are  few  in  number,  they 
generate  concentrations  which  greatly  outweigh  the  typical  concentrations  found  between 
events.  The  infrequency  of  the  high  concentration  events  is  again  apparent  in  the  bottom 
time  series  of  figure  4-12.  Shown  in  the  upper  time  series  for  comparison  is  the  square  of 
the  bottom  velocity  magnitude.  This  comparison  will  be  examined  in  more  detail  in 
section  4-4. 

Next,  the  concentration  time  series  is  brought  into  the  frequency  domain  by  use  of 
the  Fast  Fourier  Transform  (FFT).  Transformation  of  the  concentration  time  series  to  the 
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Figure  4- 1 3.  Concentration  spectrum  for  run  1 8. 
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frequency  domain  results  in  what  will  be  termed  the  concentration  spectrum.  This 
concentration  spectrum  reveals  the  relative  importance  of  each  frequencies  contribution 
to  the  total  variation  in  concentration.  Based  on  the  observations  just  made,  the  lowest 
frequencies  should  show  the  highest  contribution  to  the  total  variation  in  concentration. 
In  figure  4-13,  the  concentration  spectrum  is  plotted  from  the  time  series  of  vertically 
integrated  concentration  profiles  of  run  18,  so  the  relative  magnitudes  of  the  total 
concentration  contributions  at  each  frequency  can  be  examined.  As  expected,  the  largest 
portion  of  the  suspended  sediment  concentration  fluctuation,  or  variance,  is  accounted  for 
by  variation  at  the  lowest  frequencies. 

In  order  to  examine  this  further  and  to  look  at  the  significance  of  this  outcome  at 
different  elevations  from  the  bed,  the  cumulative  variance  function  (CVF)  is  introduced. 
The  CVF,  denoted  by  Z„,  indicates  at  a  particular  frequency  the  portion  of  the  total 
variance  accounted  for  by  lower  frequencies.  In  analytical  form,  Zxx  is  given  by  the 
following  expression: 

\sM'W       , 

z,Sf)=z =  \)sa(f'W  (4-D 

Js.W   °  ° 

0 

In  discrete  form,  the  lower  limits  of  integration  are  replaced  by  the  lowest  resolvable 
frequency,  1/7  where  7  is  the  record  duration.  Likewise,  the  upper  limit  in  the 
expression  in  the  denominator  is  replaced  by  the  Nyquist  frequency.  For  both  the 
concentration  spectrum  and  the  bottom  velocity  spectrum,  the  CVF  is  determined,  and  in 
each  case,  the  percentage  contribution  to  the  total  fluctuation  overeach  band  is  easily 
determined.  Since  the  concentration  can  be  expected  to  be  more  closely  related  to  the 
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Figure  4-14.  Bottom  velocity  spectrum  for  run  18. 
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near  bed  velocity  than  to  the  surface  elevation,  the  bottom  velocity  spectrum,  shown  in 
figure  4-14,  is  used  instead  of  a  surface  spectrum.    It  is  found  by  attenuating  the 
measured  velocity  spectrum  to  the  bed  by  use  of  linear  wave  theory  (Dean  and 
Dalrymple,  1984).  A  comparison  of  the  CVFs  generated  for  the  bottom  velocity  in  figure 
4-14  and  the  concentration  spectrum  from  figure  4-13  is  shown  in  the  first  plot  of  figure 
4-15.  In  addition,  the  second  plot  makes  the  same  comparison,  but  with  the  CVF 
generated  from  the  concentration  time  series  measured  1  cm  from  the  bed.  Because  the 
highest  concentrations  are  found  near  the  bed,  the  integrated  concentration  is  dominated 
by  the  contribution  from  the  near  bed  concentrations.  So,  it  comes  as  no  surprise  that  the 


67 


1 
0.9 

eg 

§0.7 
ni 
T5  0.6 

0 

a. 
I  0-5 

CO 

o 

r0.4 

I  0.3 
o 

CD 

u^O.2 
0.1 

0 


i 

/ 

/ 

- 

- 

- 

/  1 
^/     1 

I 

Concentration  (integrated) 

Surface  elevation 

- 

- 

1 

J 

0 


0.05 


0.1  0.15 

Frequency  (hz) 


0.2 


0.25 


0.05 


0.2 


0.1  0.15 

Frequency  (hz) 
Figure  4-15.  Cumulative  variance  functions  from  concentration  and  bottom  velocity 
spectra  for  run  18. 
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two  plots  indicate  nearly  the  same  behavior  for  the  integrated  and  near-bed 
concentrations.  In  both  cases,  a  significant  portion  (nearly  one  third)  of  the  total  variation 
in  concentration  is  accounted  for  in  the  low  frequency  band.  In  comparison,  very  little  of 
the  energy  of  the  surface  spectrum,  figure  4-9,  or  the  bottom  velocity  spectrum,  figure 
4-14,  is  found  at  low  frequencies.  In  order  to  clearly  show  the  significance  of  lower  and 
incident  frequencies,  the  plots  are  cut  off  at  0.25  Hz.  Within  the  frequency  range  from 
0.25  to  1 .00  Hz,  the  variation  in  velocity  and  in  concentration  is  uniformly  distributed. 

These  results  indicate  the  near-bed  suspended  sediment  concentration  time  series 
has  a  very  significant  low-frequency  component.  Furthermore,  the  forcing  mechanism  is 
not  apparent  at  the  low  frequencies,  since  this  region  is  poorly  represented  in  the  surface 
elevation  spectrum.  Although  figure  4-15  indicates  that  the  depth-integrated  suspended 
sediment  concentration  has  similar  behavior  to,  and  is  likely  dominated  by  the  near-bed 
concentration,  it  is  still  instructive  to  examine  the  behavior  further  from  the  bed.  Figure 
4-16  makes  the  same  comparison  as  in  figure  4-15,  but  instead  uses  concentrations 
measured  5  and  10  cm  from  the  bed.  Note  that  the  low  frequency  contribution  diminishes 
and  the  variation  in  concentration  is  more  uniformly  distributed  across  the  spectrum. 
Also,  with  increasing  height  from  the  bed,  the  contribution  to  the  total  concentration 
variation  by  frequencies  greater  than  0.25  Hz  becomes  more  relevant.  This  and  the 
diminishing  low  frequency  contribution  both  indicate  that  the  variation  becomes  more 
uniformly  distributed  with  frequency  as  the  distance  from  the  bed  increases.  Above  10 
cm,  the  concentration  variation  becomes  still  more  uniform  with  frequency,  but  the 
concentrations  become  so  small  that  signal  to  noise  ratio  of  the  concentration 
measurement  becomes  too  low. 
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Figure  4-16.  Cumulative  variance  functions  from  concentration  and  bottom  velocity 
spectra  for  run  18. 
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Figure  4-17.   Tso  indicating  less  lower  frequency  (high  period)  relevance  with  increased 
distance  from  bed  for  run  18. 
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To  further  examine  this  dependence  on  the  distance  from  the  bed,  the  frequency  at 
which  the  CVF  of  the  concentration  spectrum  is  equal  to  0.5  is  determined  at  each 
measurement  elevation  above  the  bed.  The  corresponding  period  is  that  at  which  50%  of 
the  variation  in  concentration  occurs  above  and  below,  and  is  designated  by  Tx.  A  high 
value  of  Tx  therefore  indicates  that  the  low  frequency  variation  prevails  over  the  higher 
frequencies  in  the  suspension  time  series.  A  plot  of  7"50  shown  in  figure  4-17  reinforces 
previous  observations  showing  the  diminishing  contribution  to  the  total  concentration 
variation  by  lower  frequencies  with  distance  from  the  bed. 
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4,4  Correlation  with  Velocities 


Long  waves  have  been  shown  to  drive  sediment  transport  in  the  nearshore 
environment  (Beach  and  Sternberg,  1991).  In  the  present  experiments,  the  bottom 
velocity  spectrum  and  the  surface  elevation  spectrum  show  very  little  energy  at  the  low, 
or  long  wave,  frequencies  where  significant  concentration  variation  exists.  In  this  case, 
the  mechanism  resulting  in  the  low  frequency  concentration  variation  is  not  apparent. 
Simple  examination  of  the  wave  energy  spectrum  from  a  wave  record  can  reveal  the 
frequency  components  contributing  to  the  sea  state,  but  will  not  provide  information  on 
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Figure  4-18.  Surface  displacement  time  series  from  run  18  with  envelope  determined  by 
Hilbert  transform. 
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amplitude  modulation  resulting  from  the  interaction  of  waves  at  various  frequencies.  The 
interaction  between  components  at  the  surface  results  in  low  frequency  amplitude 
modulation  of  the  wave  record,  termed  'groupiness'  due  to  the  'groups'  of  waves  formed. 
By  use  of  the  Hilbert  transform,  the  envelope  of  the  wave  record  can  be  found  (Haller 
and  Dalrymple,  1995).  Such  an  envelope  is  shown  superimposed  on  a  portion  of  a 
surface  elevation  time  series  in  figure  4-18.  Wave  groups  are  generally  believed  to 
contribute  to  the  forcing  of  long  waves  (Longuet-Higgins  and  Stewart,  1964)  and  the 
suspension  of  sediment  in  the  nearshore  region  (Hanes,  1994).  Because  visual 
observations  at  the  time  of  the  experiments  indicated  that  wave  groups  were  present,  it  is 
speculated  that  the  same  interactions  that  result  in  wave  groups  result  also  in  the  low- 
frequency  suspension  events  seen  in  the  previous  section.  In  the  following,  the  source  of 
these  interactions  is  examined. 

Begin  by  considering  the  sum  of  only  two  components  of  slightly  different 
frequency. 

T)  =  a  cos(&,  x  -  CT,r)  +  a  cos(k2x  -  G2t)  (4-2) 

The  wave  numbers  and  angular  frequencies  alternatively  can  be  represented  by  the 
following: 

k  -F-M 

*,=*  +  **  (4"3) 

2 


and 
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_    Act 
ct,  =  a  +  — 

2 

Substituting  these  expressions  and  simplifying  gives 

(Ak       Act  ^      ,-      _> 
rj  =  2acos  —x — —t  \coS(kx  -  at  j 

Squaring  this  expression  and  simplifying  shows  each  interaction  term. 


7J-  =a- 


+  cos(Afoc  -  Aor)  +  cos(2(fct  -  at))  +  4-cos((2£  +  Ak)x  -  (2ct  +  Act). 


(4-4) 


(4-5) 


(4-6) 


+  icos((2t  - M).v-(2ct  -  Act>) 

Examining  the  terms  within  brackets  one  at  a  time,  the  first  term,  being  independent  of 
frequency,  is  simply  a  'DC  offset.  The  second  term,  of  greatest  interest  in  this  study,  can 
be  rewritten  as 

cos^-^-fo-o-Jr)  (4-7) 

which  clearly  shows  this  term  results  from  the  the  difference  of  the  original  frequency 
components.  Similarly,  the  second  and  third  term  can  be  rewritten  as 

cos((jfc;  +  fe,  )x  -  (a,  +a,)t)  (4"8) 

which  again  clearly  shows  this  term  results  from  the  sum  of  the  original  frequency 
components.  The  final  two  terms  are  harmonics  of  the  two  original  frequency 
components,  respectively,  as  can  be  seen  by  rewriting  the  fourth  term  as  follows: 

icos(2£2;t-2cv)  (4-9) 


Of  the  four  components  mentioned  above,  the  frequency  difference  between  components 
results  in  wave  groups.  In  order  to  examine  the  frequencies  of  the  wave  groups,  the 
envelope  presented  in  figure  4-18  is  transformed  into  the  frequency  domain, 
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Figure  4- 1 9.  Spectrum  of  envelope  from  run  1 8. 
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resulting  in  the  envelope  spectrum  shown  in  figure  4-19.  From  figure  4-19,  it  is  apparent 
that  the  frequencies  best  represented  in  the  envelope  spectrum  are  the  same  as  those 
represented  in  the  lowest  portion  of  the  concentration  spectrum,  figure  4-13.  Re- 
examination of  figure  4-13  also  reveals  an  active  range  of  concentration  variation  at 
frequencies  higher  than  the  incident  wave  frequencies.  The  components  in  this  range 
likely  result  from  the  harmonics  and  the  frequency  sum  terms  discussed  above.  This 
further  supports  the  conjecture  that  interaction  terms  contribute  significantly  to  sediment 
suspension. 

Considering  that  all  of  the  interactions  are  seen  in  the  concentration  variation,  it  is 
worthwhile  to  examine  the  relation  between  the  square  of  the  bed  velocity  magnitude  and 
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the  near  bed  concentration.  As  before,  the  bed  velocity  magnitude  is  determined  by 
combining  both  horizontal  velocity  components,  which  are  found  by  attenuating  the 
velocity  measurements  using  linear  wave  theory.  Concentration  measurements  made  1 
cm  above  the  bed  are  utilized  in  the  following  analysis.  To  this  point,  the  results  shown 
in  the  figures  have  been  restricted  to  those  from  run  18  for  clarity.  It  should  be  noted  that 
the  results  from  the  other  3  runs  used  in  this  investigation  show  still  greater  significance 
in  the  lowest  frequencies  of  the  concentration  spectra.  Application  of  the  following 
methods  of  analysis  will  verify  this,  and  this  analysis'  results  from  data  runs  23,  25,  and 
27  will  be  shown  in  figures  4-24  through  4-26  respectively. 

In  order  to  examine  the  relation,  the  coherency  function  is  utilized.  The 
coherency  function  indicates  whether  one  signal  can  be  expressed  as  a  linear  function  of 
another  signal.  It  is  calculated  by  performing  auto  and  cross-spectral  analysis  over 
sections  of  the  signals,  and  then  determining  the  linear  relation  between  sections. 

{c>)}2+{e>)?  (4-'°) 


y(co) 


S„(ft))Sn.(<o) 


In  this  expression,  Sxx  and  S„.  are  the  autospectral  densities  of  the  respective  signals,  and 
Cxy  and  Q„,  are  the  cospectrum  and  quadrature  spectrum  -  the  real  and  imaginary 
components  of  the  cross-spectral  density  function.  Ochi  (1990)  provides  thorough 
explanations  concerning  the  development  and  application  of  each  of  these  functions  as 
well  as  the  additional  spectral  analysis  techniques  presented  in  this  dissertation.  If  the 
same  linear  relation  holds  between  various  sections,  the  coherency  function  will  return  a 
value  of  one.  If  the  sections  are  related  by  nearly  linear  relations,  the  coherency  function 
will  still  return  a  value  close  to  one.  Should  the  sections  have  completely  different  linear 
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relations  or  should  no  linear  relation  exist  for  certain  sections,  the  coherency  function  will 
be  zero.    Ideally,  the  signal  should  be  broken  into  as  many  sections  as  possible, 
providing  many  degrees  of  freedom  in  the  analysis.  In  addition,  it  is  desirable  to  obtain 
good  frequency  resolution  at  the  lower  frequencies  in  the  spectral  analysis,  requiring  that 
the  sections  be  sufficiently  long.  In  the  present  case  the  length  of  the  data  files  was 
limited  by  the  available  memory  in  the  data  logger.  Consequently,  the  number  of  degrees 
of  freedom  in  the  analysis  is  limited  to  approximately  10,  which  allows  a  maximum 
period  of  3  minutes  to  be  resolved.  Coherency  function  values  of  better  than  60%  were 
found  in  most  of  the  comparisons  between  the  square  of  the  velocity  and  the  near  bed 
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Figure  4-20.  Coherency  function  between  square  of  bed  velocity  magnitude  and  near  bed 
concentration  for  run  18. 
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concentrations,  particularly  at  the  lowest  frequencies  in  the  spectrum,  as  seen  in  figure 
4-20  and  plot  (c)  of  figures  4-24  through  4-26.  There  is  significant  variation  in  the 
magnitude  of  the  coherency  function  across  frequencies,  due,  in  part,  to  the  limited 
degrees  of  freedom  in  the  analysis.  In  addition,  it  is  difficult  to  determine  a  linear 
relation  in  the  regions  of  the  spectra  in  which  there  is  little  energy,  due  to  limited 
instrument  resolution.  Therefore,  the  value  of  the  coherency  function  in  these  regions  is 
not  as  meaningful.  For  this  reason,  the  results  of  the  coherency  analysis  are  presented  in 
another  form.  In  figure  4-21,  the  concentration  spectrum  for  the  nearbed  concentration  is 
shown.  The  curve  is  shown  as  solid  in  those  regions  in  which  the  coherency  function, 
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Figure  4-2 1 .  Concentration  spectrum  (dotted)  with  areas  of  coherency  function  >  60% 
indicated  (solid)  for  run  18. 
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determined  from  the  near  bed  concentration  and  square  of  the  bed  velocity  magnitude, 
returned  a  value  greater  than  60%.  These  marked  regions  are  the  most  dominant  regions 
of  the  concentration  spectrum,  covering  a  total  of  76%  of  the  total  variation  of  the 
spectrum.  This  indicates  that  a  linear  relation  between  the  velocity  squared  and  near  bed 
concentration  is  likely.  In  each  of  figures  4-24  through  4-26,  plot  (b)  verifies  this  result, 
showing  high  coherence  in  the  portions  of  the  spectrum  with  the  highest  magnitudes  of 
variance.  Though  only  run  18  showed  significant  variation  in  concentration  at  incident 
wave  band  frequencies,  it  should  be  noted  that  significantly  higher  coherence  was  found 
in  this  region  as  the  number  of  degrees  of  freedom  in  the  analysis  was  increased.  Since 
increasing  the  number  of  degrees  of  freedom  simply  means  increasing  the  number  of 
instances  used  in  the  evaluation,  this  could  suggest  that  the  statistical  means  of  the 
processes  are  linearly  related  at  these  frequencies.  More  runs  with  significant  variation  in 
this  frequency  band  need  to  be  analyzed  to  confirm  such  a  hypothesis. 

Since  the  coherency  analysis  indicates  a  linear  relation  exists  between  the  square 
of  the  bed  velocity  magnitude  and  the  near  bed  concentration,  it  is  interesting  to  examine 
the  transfer  function. 

Y{o))=X(w)H{(o)  (4-11) 

The  transfer  function,  H(co),  is  simply  a  function  which  relates  a  linear  system's  input, 
X(co),  the  square  of  the  bed  velocity  magnitude  in  this  case,  to  the  system's  output,  Y(m), 
the  near  bed  concentration  (Ochi,  1990).  In  the  frequency  domain,  the  input  signal 
multiplied  at  each  frequency  by  the  transfer  function  produces  the  output  signal. 
Considering  that  the  input  and  output  signals  can  vary  in  both  magnitude  and  phase,  the 
transfer  function  evaluated  at  a  particular  frequency  in  general  produces  a  complex 
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number.  The  phase  variation  is  indicative  of  the  time  lag  between  the  processes.  Figure 
4-22  shows  the  phase  of  the  transfer  function  calculated  for  run  18.  Examination  of 
figure  4-22  and  plot  (d)  of  figures  4-24  through  4-26  show  a  relatively  small,  and 
typically  negative,  transfer  function  phase.  There  is  a  slight  trend  to  decrease  in  phase  as 
frequency  increases  in  the  low-frequency  band  of  the  relation.  Such  a  trend  indicates  that 
the  time  lag  between  the  signals  is  nearly  constant.  It  is  expected  that  this  time  lag  is 
related  to  the  upward  sediment  flux  from  mixing  and  to  the  fall  velocity  of  the  sediment. 
Future  investigations  which  include  measurements  of  sediment  size  will  aid  in  verifying 
such  a  relation. 
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Figure  4-22.  Phase  of  transfer  function  for  run  18.  Negative  phase  indicates 
concentration  lags  square  of  velocity  magnitude. 
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Figure  4-23.  Magnitude  of  transfer  function  for  run  18. 


Finally,  in  figure  4-23  and  in  plot  (e)  of  figures  4-24  through  4-26,  the  magnitude 
of  the  transfer  function  is  shown.  For  each  particular  run,  the  magnitude  of  the  transfer 
function  at  the  frequencies  in  which  the  coherence  was  high  varies  little  with  frequency. 
The  relative  difference  in  the  transfer  function  magnitudes  for  different  runs  is  attributed 
to  the  difficulty  in  quantifying  the  bed  location  in  the  concentration  profiles.  In  other 
words,  because  the  spatial  resolution  of  the  concentration  measurement  was  limited  to 
0.75  centimeters  by  the  sampling  rate  of  the  system,  so  was  the  resolution  in  determining 
the  bed  location.  With  an  exponential  height  variation  in  concentration,  a  small  error  in 
determining  the  bed  location  can  result  in  a  significant  change  in  concentration.  For 
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example,  referring  back  to  figure  4-8,  the  concentration  changes  from  1 .5  g/1  at  the 
closest  measurement  to  the  bed  to  approximately  0.5  g/1  at  the  next  closest  measurement 
0.75  cm  higher.  Regardless,  it  is  promising  to  see  that  the  value  is  nearly  constant  across 
frequency  for  a  particular  run,  and  that  this  value  is  of  similar  magnitude  between  runs.  It 
is  expected  that  this  parameter  too  is  related  to  sediment  parameters  and  local 
hydrodynamic  conditions,  such  as  bed  roughness  or  bedforms.  In  relation  with  the 
calculation  of  the  time  mean  reference  concentration,  future  research  could  attempt  to 
relate  the  magnitude  of  the  sediment  resuspension  coefficient,  since  a  linear  relation 
seems  likely  (Smith  and  McLean,  1977;  Thosteson,  1995). 

4.5  Discussion  of  Results 
Although  nearly  every  possible  complexity  preventing  accurate  evaluation  of  the 
median  grain  diameter  of  suspended  sediment  was  considered  there  were  insurmountable 
difficulties  in  determining  grain  size  from  the  data  collected  in  this  experiment. 
Horizontal  gradients  in  concentration  between  measurements  are  believed  to  be 
responsible  for  the  difficulties  in  determining  grain  size.  Vortex  ripples  were  measured 
in  every  data  run,  and  may  be  responsible  for  this  horizontal  variation  in  concentration. 
Because  of  the  difficulties  encountered,  the  constant  grain  size  assumption  was  made  in 
the  remainder  of  the  analysis.  Concentration  measurements  from  the  2.25  MHz 
transducer  were  used  as  opposed  to  those  from  the  5.0  MHz  transducer,  since  for  the 
grain  sizes  encountered  in  the  region,  the  2.25  MHz  unit  shows  less  sensitivity  to  grain 
size.  It  should  be  noted  that  errors  resulting  from  this  assumption  are  likely  small  in  the 
present  analysis,  because  the  active  portion  of  the  profile  remained  small.  Namely,  the 
concentration  profile  stayed  confined  to  the  area  very  close  to  the  bed,  so  analysis  errors 
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Figure  4-24.  Results  from  run  23:  (a)  bottom  velocity  spectrum;  (b)  concentration 
spectrum;  (c)  coherence;  (d)  transfer  function  phase;  (e)  transfer  function  magnitude. 
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Figure  4-25.  Results  from  run  25:  (a)  bottom  velocity  spectrum;  (b)  concentration 
spectrum;  (c)  coherence;  (d)  transfer  function  phase;  (e)  transfer  function  magnitude. 
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Figure  4-26.  Results  from  run  27:  (a)  bottom  velocity  spectrum;  (b)  concentration 
spectrum;  (c)  coherence;  (d)  transfer  function  phase;  (e)  transfer  function  magnitude. 


85 

which  can  propagate  down  through  the  profile  didn't  have  far  to  propagate.  Analysis  of 
the  concentration  spectrum  shows  that  low  frequency  variation  is  as  well  represented  as 
variation  in  the  incident  wave  band.  Failing  to  account  for  the  low  frequency  variation  in 
concentration  will  result  in  neglecting  a  large  percentage  of  the  variation.  Considering 
that  the  surface  elevation  spectrum  and  near  bed  velocity  spectrum  show  no  infragravity 
component,  the  sediment  suspension  was  not  driven  by  long  waves.  The  frequencies  of 
the  variation  coincide  well  with  the  incident  band  frequencies  and  the  resulting 
frequencies  from  interactions  between  components  of  the  incident  band.  This  indicates 
that  the  groupiness  of  the  waves  is  responsible  for  the  low  frequency  variations  in 
concentration.  Physically,  several  explanations  are  possible.  Sediment  may  be  more 
readily  suspended  after  several  consecutive  large  waves  than  after  only  a  single  large 
wave,  perhaps  due  to  group  enhanced  fluidization  of  the  bed.  Possibly  the  time  required 
for  the  sediment  to  settle  introduces  a  sort  of  memory  effect,  such  that  each  consecutive 
large  wave  of  a  group  will  add  more  sand  to  that  currently  suspended.  Likewise, 
turbulence  generated  at  the  bed  may  build,  changing  the  boundary  layer  structure  and  the 
settling  properties  of  the  sediment. 

Linear  coherence  between  the  square  of  the  velocity  and  the  near  bed 
concentration  is  found  across  most  of  the  frequency  spectrum.    Because  many  attempts 
at  modeling  sediment  transport  treat  the  wave  conditions  by  a  single  representative  wave 
height  and  period,  no  low  frequency  components  from  interactions  will  result.  Not 
surprisingly,  there  has  been  little  success  in  predicting  an  instantaneous  concentration 
time  series.  Should  linear  interactions  between  spectral  components  be  considered,  the 
present  results  seem  to  indicate  prediction  of  the  concentration  time  series  is  promising. 
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Analysis  of  the  transfer  function  in  this  investigation  was  somewhat  limited,  but 
the  results  show  a  linear  trend  in  phase,  that  being  an  increasing  phase  difference  with 
increasing  frequency.  A  transfer  function  with  a  linear  phase  shift  indicates  a  constant 
time  lag  between  the  input  and  output  signals.  Because  the  length  of  the  data  records 
were  limited  in  this  study,  the  exact  magnitude  of  this  time  lag  can  not  be  determined 
with  confidence.  Finally,  the  magnitude  of  the  transfer  function  showed  little  frequency 
dependence  for  a  particular  ran,  but  different  runs  showed  different  transfer  function 
magnitudes.  These  differences  are  likely  related  to  changes  in  the  hydrodynamic 
conditions,  but  are,  in  part,  explained  by  limited  measurement  resolution.  Further 
research  will  be  necessary  to  draw  more  definite  conclusions  about  the  variation  in  this 
parameter. 


CHAPTER  5 
CONCLUSIONS 

In  this  dissertation,  a  system  of  oceanographic  instrumentation  capable  of 
providing  a  broad  set  of  measurements  relevant  to  littoral  sediment  transport  processes  is 
presented.  Instrumentation  includes  a  multiple  frequency  acoustic  backscatter  system 
(ABS),  multiple  frequency  multiple  transducer  arrays  (MTAs),  acoustic  Doppler 
velocimeters  (ADVs),  an  optical  backscatterance  sensor  (OBS),  pressure  sensor, 
temperature  sensor,  tilt  and  compass  sensors,  and  an  underwater  video  camera.  All 
instruments  are  in  close  proximity  to  the  data  system,  maximizing  data  integrity.  The 
system  includes  several  self-monitoring  functions,  including  power  monitoring,  leak 
detection,  internal  temperature  measurement,  and  communications  monitoring.  The 
entire  unit,  which  is  deployed  underwater,  is  capable  of  self-operation  should  the  external 
power  and  communication  connection  be  compromised.  Although  externally  powered, 
the  system  contains  a  reserve  of  battery  power  for  current  intensive  operations  and  for 
operation  in  the  event  of  the  loss  of  external  power.  A  user  interface  and  dynamic 
memory  configuration  allow  both  real-time  testing  of  the  instrumentation  and  broad 
flexibility  in  the  definition  of  sampling  schemes.  Data  collected  by  the  system  are 
reduced  through  highly  efficient  analysis  algorithms  during  data  collection,  and  are  both 
stored  internally  on  the  systems  internal  storage  and  exported  in  nearly  real  time  to  shore 
during  collection  through  serial  communications.  Those  data  stored  on  the  system  are 
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offloaded  using  an  error  correction  protocol  to  a  compact  binary  data  file  that  can  be 
immediately  viewed  and  interpreted  by  analysis  software. 

In  order  to  interpret  the  acoustic  backscatter  measurements  obtained  by  the 
system,  a  new  technique  of  determining  both  concentration  and  the  median  grain 
diameter  of  suspended  particles  has  been  presented.  The  significant  advantage  of  the 
technique  is  that  by  using  an  explicit  solution  for  concentration,  the  median  grain 
diameter  can  be  found  without  having  to  neglect  sediment  attenuation.  In  addition, 
because  incorporating  the  correct  median  grain  diameter  in  the  explicit  solution  will 
produce  an  identical  concentration  regardless  of  the  operational  frequency  of  the 
transducer,  the  concentration  variance  between  any  number  of  transducers  can  be 
minimized  to  find  the  median  grain  diameter.  Previous  techniques  have  typically  paired 
transducers  and  applied  an  empirical  procedure  on'the  results  from  the  pairs  to  obtain  the 
grain  diameter.  Numerical  simulations  were  performed  to  verify  the  capability  of  the 
technique  to  produce  both  the  expected  concentrations  and  grain  diameters.  In  addition, 
laboratory  results  from  a  recirculating  calibration  chamber  verify  that  the  technique 
applies  well  in  determining  sediment  size  and  concentration  from  measurements  of 
backscattered  acoustic  intensity. 

Finally,  a  system  similar  to  that  described  in  the  text  was  used  in  a  data  collection 
project  utilizing  the  SIS  of  the  Army  Corp's  Field  Research  Facility.  Difficulties 
encountered  in  trying  to  obtain  the  median  grain  size  of  the  suspended  sediment, 
particularly  due  to  the  spatially  disjoint  measurements  of  the  acoustic  transducers,  led  to  a 
discussion  of  field  measurement  of  suspended  sediment  grain  size.  In  the  present 
analysis,  the  spatially  separated  beams  resulted  in  inconsistent  results  when  attempting  to 
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determine  grain  size.  Because  of  the  sensitivity  of  any  size-evaluating  technique  to  small 
concentration  changes,  it  is  believed  the  inconsistencies  were  due  to  varying  horizontal 
gradients  in  concentration  between  the  beams  of  the  transducers. 

The  remainder  of  the  analysis  is  devoted  to  the  determination  of  the  relevant  time 
scales  of  sediment  suspension.  Near  the  bed,  variation  in  concentration  at  low 
frequencies  is  as  common  as  variation  at  the  frequencies  of  the  incident-band  waves,  but 
the  near  bed  velocity  spectrum  indicates  little  infragravity  wave  forcing  at  these 
frequencies.  Concentration  variation  is  found  to  occur  at  frequencies  consistent  with 
those  predicted  by  the  linear  interactions  of  incident-band  components.  At  the  lowest 
frequencies,  these  linear  interactions  result  in  groupiness  of  the  surface  elevation  time 
series  as  well  as  groupiness  in  the  velocity  time  series  at  the  bottom.  The  high  coherency 
function  values  obtained  across  most  of  the  frequency  spectrum  when  relating  the  near 
bed  concentration  to  the  square  of  the  bottom  velocity  indicate  good  linearity  between  the 
two  signals.  Inclusion  of  more  details  of  the  hydrodynamics  such  as  boundary  layer 
structure  and  bottom  roughness  will  likely  improve  the  coherence  and  make  possible 
instantaneous  prediction  of  concentration.  Because  a  significant  portion  of  the  near  bed 
variation  in  concentration  seems  to  be  related  to  the  linear  interaction  between  incident 
band  components,  parameterization  of  the  incident  band  waves  by  a  single  frequency  will 
certainly  result  in  tremendous  error  in  any  attempt  to  predict  the  instantaneous 
concentration. 

Limitations  in  the  present  work  were  encountered  in  each  step  of  the 
investigation.  First,  the  system  used  in  the  SIS96  experiments  limited  the  duration  of 
data  collection,  limiting  the  confidence  in  the  spectral  analysis  techniques  presented  in 
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chapter  4.  Since  the  duration  of  data  collection  is  much  longer  in  the  newer  system, 
described  in  chapter  2,  future  data  sets  should  provide  higher  confidence  and  improved 
results.  Next,  concentrations  in  the  recirculating  calibration  chamber  used  for  calibration 
of  the  acoustic  transducers  and  OBSs  were  not  verified  by  suction  sampling.  It  is 
expected  that  the  errors  in  both  concentration  and  median  grain  size  determination  will 
decrease  when  suction  sampling  becomes  a  routine  part  of  this  calibration.  Furthermore, 
the  empirical  form  function  is  assumed  to  describe  any  noncohesive  quartz  sediment. 
This  assumption  should  be  verified  by  empirical  determination  of  the  form  function  for 
the  sediment  found  on  the  site  of  the  experiments.  Regarding  the  difficulties  encountered 
in  the  size  evaluation,  collocated  acoustic  beams,  perhaps  in  the  transducer  design  itself, 
will  simplify  the  evaluation  and  improve  the  results.  Due  to  the  limited  time-span  of  the 
SIS96  project  and  the  few  long-duration  data  files  collected,  the  results  presented  in 
chapter  4  are  based  on  a  restricted  set  of  hydrodynamic  conditions.  Continued 
experimentation  and  analysis  of  other  data  sets  offering  a  wider  range  of  conditions  are 
needed  for  verification  of  these  results.  Finally,  the  results  should  be  verified  at  other 
locations  to  ensure  effects  of  the  pier  and  supporting  structure  are  minimal. 


APPENDIX  A 
ACOUSTIC  PARAMETERS 

In  the  following,  for  completeness,  the  expressions  used  for  calculation  of  various 
parameters  present  in  the  backscatter  equations  are  presented. 

A.  1   Log-normal  Distribution 
If  the  median  grain  size,  /ii  ,  and  standard  deviation,  <T4 ,  are  expressed  in  units  of 

<p,  the  probability  distribution  in  terms  of  the  grain  radii,  a,  is  given  in  the  subsequent 
formula. 


p{a)= 


(A-l) 


A.2  System  Constant 
The  system  constant,  A ,  is  determined  by  calibrating  each  transducer  in  a 
recirculating  calibration  chamber  with  a  series  of  suspensions  of  constant  concentrations 
of  a  known  grain  size  distribution  of  sediment. 


A. 3  Water  Attenuation 
The  water  attenuation  parameter  is  determined  in  seawater  by  the  expression  of 
Fisher  &  Simmons  (1977). 
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a   =  (55.9-2.377  +  0.0477r2-0.000384r,)/M2xlO";' 

(A-2) 


In  equation  (A-2),  fu  is  the  frequency  in  mega-Hertz  and  T  the  temperature  in  degrees 
Celsius. 

A.4  Sediment  Attenuation 
The  sediment  attenuation  parameter  is  determined  following  the  method  outlined 
by  Thome  (1993). 

i=^y}x)  (a-3) 

The  pointed  brackets  represent  a  value  representative  for  a  distribution  of  sediment  with 
probability  density  function,  p(a,j,  rather  than  for  a  single  sediment  size,  so  the 
representative  value  for  a, ,  the  particle  radius,  is  given  in  equation  (A-4). 


(fl»)  =  Ja»i'{«.)ia. 


(A-4) 


In  equation  (A-3),  %  is  the  normalized  total  scattering  cross  section,  shown  in  figure 
A-l,  and  given  by  the  following: 


(A-5) 
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x  is  the  dimensionless  particle  radius,  given  as  the  product  of  the  particle  radius  and  the 
acoustic  wave  number.    Ka  is  a  constant  based  on  the  particle  density  and 
compressibility  and  is  equal  to  0.18  for  quartz  sediment. 
10° 


10  10  10 

Dimensionless  particle  radius,  x 

Figure  A- 1 .  Normalized  total  scattering  cross  section. 


Finally,  equation  (A-6)  gives  the  representative  value  of  %  for  a  distribution  of  sediment. 


<*>-*■ 


\a,p(at  )ia,\a?xp(a,  Ya, 


I  a,  p(aK)dat 


(A-6) 
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A.4  Backscatter 
Equation  (A-7)  gives  the  backscatter  function. 


(A-7) 


In  this  equation,  the  form  function,  as  defined  by  Thorne  (1993),  is  given  by 

K,x-    1 


/  =  C„ 


\  +  Kfx' 


(A-8) 


with 


l-v,e  L  * 


1+V,e  L  ,J 


(A-9) 


Figure  A-l  shows  the  form  function  as  given  in  equation  (A-8).  Note  the  form  function 
is  largest  in  magnitude  when  the  acoustic  wavelength  is  equal  to  the  particle  diameter,  or 
when  the  nondimensional  particle  radius  is  equal  to  n.  In  this  equation, 
v  ,=  0.37, 77  ,=  0.5,  x,  =  1.4,  v  ,=  0.28, 77  2=  2.2,  x2  =  2.8  based  on  an  empirical  fit  to 
measured  data.  And  based  on  the  compressibility  and  density  of  the  sediment  relative  to 
water,  K,  =  1.1 4  for  non-cohesive  sedimentary  material.  The  representative  value  of  the 
form  function  for  a  distribution  of  sediment  is  given  by  the  following: 


(/>«*- 


I  a,  p(fl,  pa,  J  a,  fp(a,  \lat 


\as~  p{as)da, 


(A- 10) 
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Figure  A-2.  Form  function. 
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A.5  Nearfield  Correction 
Equation  (2-1)  presents  the  acoustic  backscatter  equation  with  the  nearfield 
correction  term,  \jf .  This  term,  as  defined  by  Thorne  et  al  (1991),  is  given  as  follows: 


V/  =  l 


for  z  >  £  z„(farfield) 
for  z  <  £  z„  (nearfield) 


(A-ll) 


In  the  original  text,  £  is  set  equal  to  2.  z„  is  the  theoretical  limit  of  the  nearfield,  and  is 
given  by  equation  (A-12). 
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K  a 

z„  =-T-  (A- 12) 


In  (A-12),  X  is  the  acoustic  wavelength  and  a,  is  the  radius  of  the  transducer's 
piezoelectric  crystal.  The  explicit  solution  for  concentration  including  the  nearfield 
correction,  equation  (A- 13),  is  found  following  the  same  derivation  presented  above. 

i — —>—  exp(4av„z) 

C  = ' (V~)V (A"13) 

y-  J4C    z    V   exp(4aj)dz 

t 
In  this  solution,  the  integration  constant  is  also  slightly  modified. 

r=(v^Mexp(4fVo)  (Ai4) 

Like  the  other  subscripted  variables,  \f/Q  is  determined  at  zo- 


APPENDIX  B 
LIST  OF  VENDORS 


Product:  2.25  MHz,  5.00  MHz  810  Special  Echo  Sounders 

Model  900  rotating  scanning  sonar  system 

Simrad  Mesotech  Systems  Ltd. 

1598  Kebet  Way 

Port  Coquitlam,  B.C.  V3C  5W9 

Phone:(604)464-8144 

Fax:(604)941-5423 


Product:  Acoustic  Backscatter  System  (ABS) 

Centre  for  Environment,  Fisheries  and  Aquaculture  Science  (CEFAS, 

formerly  known  as  The  Fisheries  Laboratory) 
N.  Pearson 
Lowestoft 
Suffolk 
NR33  OHT 
UK 

Phone: +44  1502  524319 
Fax: +44  1502  524318 


Product:  Acoustic  Doppler  Velocimeter  (ADV) 

SonTek,  Inc. 

6837  Nancy  Ridge  Drive,  Suite  A 
San  Diego,  CA  92121 
Phone:(619)546-8327 
Fax:(619)695-8131 


Product:  Micro-SeaCam  1050 

DeepSea  Power  &  Light 

48 1 9  Ronson  Court 

San  Diego,  CA  921 11-1803 

Phone:(619)576-1261 

Fax:(619)576-0219 
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Product;  Multiple  Transducer  Array  (MTA) 

Seatek 

5108  NW  27lh  Drive  32605 
Gainesville,  FL  32605 
Phone:(352)337-1393 


Product:  MUX32  Multiplexor  board 

Daedulus  Research,  Inc. 
1780  Research  Parkway,  Suite  108 
Logan,  UT  84321 
Phone:(801)752-2141 
Fax:(801)753-1971 


Product:  Optical  Backscatterance  Sensor  (OBS)  model  3 

D&A  Instruments  and  Engineering 
40-A  Seton  Rd. 
Port  Townsend,  WA  98368 
Phone:  (800)  437-8532 


Product:  Pressure  Sensor 

TransMetrics,  Inc. 
5325  Naiman  Parkway 
Solon,  OH  44139 
Phone:(216)248-2229 
Fax:(216)248-7780 


Product:  Tattletale  model  7,  model  8 

Onset  Computer  Corporation 
536  MacArthur  Blvd 
P.O.  Box  3450 
Pocasset,  MA  02559-3450 
Phone:  (508)  563-9000 
Fax: (508) 563-9477 


Product:  Underwater  cable 

Mercury  Wire  Products,  Inc. 
Mercury  Dr. 
Spencer,  MA  01562 
Phone:  (508)  885-6363 
Fax:(508)885-3316 


APPENDIX  C 
MONLOG  1.0  PROGRAM  LISTING 


MONLOG  1.0   -   Revised  version  of  DuckLog  1.1,  with  addition  of 

MTA  communications  routines  and  parallel  port  output  of  data  for 
real  time  data  output  to  shore.   In  addition,  analog  sampling  is 
modified  appropriately  to  match  new  ADV's  requirements  as  well 
as  additional  hardware  modifications. 

Corrections  have  been  made  to  the  memory  allocation  routine  in 
agreement  with  Express  3.0 

Initial  release  -  Eric  Thosteson  -  started  9/17/96  */ 


Include  files 


ttinclude  <atod.h> 
#include  <datafile.h> 
ttinclude  <drive.h> 
ttinclude  <offload.h> 
#include  <pario.h> 
ttinclude  <pins.h> 
ttinclude  <psram.h> 
ttinclude  <rtc.h> 
ttinclude  <time.b> 
ttinclude  <timing.h> 
ttinclude  <tserial.h> 
ttinclude  <tpu.h> 
ttinclude  <userio.h> 

*  Defines 


ttdefine  CLEAR  "\xlB[2J" 
ttdefine  UNDON  "\xlB[4m" 
ttdefine  UNDOFF  "\xlB[0m" 
ttdefine  FLSHON  "\xlB[5m" 
ttdefine  FLSHOFF  "\xlB[0m" 

ttdefine  VBUSCAL  7.292e-3 
ttdefine  VREGCAL  7.282e-3 
ttdefine  VUNREGCAL  1.856e-2 

ttdefine  MTA_MAXSCANS  3406 


Global  variable  declarations 


bool  d3;  /*  mux  select  -  LSB  */ 

bool  d4;  /*  mux  select  */ 

bool  d5;  /*  mux  select  -  MSB  */ 

bool  e0;  /*  adv  #1  power  on/off  */ 

bool  el;  /*  camera  power  on/off  */ 

bool  e2;  /*  tpu  buffer  enable  */ 

bool  e5;  /*  charge  on/off  */ 

bool  e6;  /*  tilt,  compass,  temp  on/off 
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boo 1  e7 ; 

bool  fl; 
bool  £6; 
bool  £7; 

bool  df7; 
boo 1  firs  t_t ime ; 
bool  storm; 
bool  tpuO ; 

char  acquisition_time[12] ; 
char  ch; 

char  instring[201 ; 
char  mta_prog[35] ; 
char  mta_time [20] ; 
DriveErr  de; 
int  data_file_pointer; 
int  num_to_f ill; 
int  samp_mode; 
int  stringerror; 
int  string/point; 
short  inchan=8; 
short  outchan-7; 
ucpv  pioport; 
unsigned  long  tmal28; 
unsigned  long  'sqbase; 
unsigned  long  *rtbase; 
unsigned  long  *flsum; 
unsigned  long  *flsqr; 
unsigned  long  *spsum; 
unsigned  long  *spsqr; 
unsigned  long  *f2sum; 
unsigned  long  *f2sqr; 
unsigned  long  *f3sum; 
unsigned  long  *f3sqr; 
unsigned  long  *flhed; 
unsigned  short  *flrms; 
unsigned  long  *f2hed; 
unsigned  short  *f2rms; 
unsigned  long  *f3hed; 
unsigned  short  *f3rms; 
unsigned  long  *sphed; 
unsigned  short  *spmen; 
unsigned  long  "mtlhed; 
unsigned  short  'mtlbas; 
unsigned  long  *mt2hed; 
unsigned  short  *mt2bas; 
unsigned  long  *mt3hed; 
unsigned  short  *mt3bas; 
unsigned  long  *tmhed; 
unsigned  short  *tmbas; 
unsigned  long  *dmhed; 
unsigned  short  "dmbas 
unsigned  long  *tflbas 
unsigned  long  *tf2bas 
unsigned  long  *tf3bas 
unsigned  long  *tspmen 
unsigned  short  mode-3 
unsigned  long  nosp; 
unsigned  long  perm_npp; 
unsigned  long  npp; 
unsigned  long  perm_prof_pe 
unsigned  long  prof_period; 
unsigned  long  perm_npa; 
unsigned  long  npa; 
unsigned  long  perm_nupc,- 
unsigned  long  nupc; 
unsigned  long  n£p; 
unsigned  long  shore_nfp; 
unsigned  long  total_nfp; 
unsigned  long  mem„nfp; 
unsigned  long  dif_nfp; 
unsigned  long  storflag; 


/*  obs  &  pressure  power  on/off  */ 

/*  ADV  #2  power  on/off  */ 

/*  3  freq.  ACP  power  on/off  */ 

/*  MTA  power  on/off  */ 

/*  MTA  disable  */ 

/*  first  time  through  data  collection  flag  */ 

/*  storm  mode  flag  -  1-enabled  */ 

f*  ABS  gain  select  */ 

/*  string  containing  total  time  to  fill  file  * 

/*  location  for  tt4  character  input  */ 

/*  input  string  for  tt4  communications  */ 

/*  program  string  for  mta  *7 

/*  time  string  for  syncing  tt4  */ 

/*  disk  drive  error  flag  */ 

/*  disk  data  file  pointer  */ 

/*  number  of  data  files  to  fill  */ 

/*  timed /continuous  sampling  mode  */ 

/*  error  code  of  mta  com  error  */ 

/*  current  position  of  in  &  outgoing  tt4  com  * 

/*  TPU  channel  for  mta  com  input  */ 

/*  TPU  channel  for  mta  com  output  */ 

/*  Reserve  location  for  par.  i/o  pointer  */ 

/■"  tma  -  128  */ 

/*  base  address  of  x"2  lookup  table  */ 

/*  base  address  of  (x+0.5)*2  lookup  table  */ 


base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  address  of  sum 
base  addres: 
base  addres 


x)  for  Fl 

x"2)  for  Fl  */ 

x)  for  SPI  channels  */ 

x"2)  for  SPI  channels 

x)  for  F2  */ 

x"2)  for  F2  */ 

x)  for  F3  */ 

x"2)  for  F3  */ 
for  Fl  header  */ 
for  Fl  rms  data  */ 


/*  base  address  for  F2  header  */ 

/*  base  address  for  F2  rms  data  */ 

/*  base  address  for  F3  header  */ 

/*  base  address  for  F3  rms  data  */ 

/*  base  address  for  SPI  header  */ 

/*  base  address  for  SPI  mean  data  */ 

/*  base  address  for  mta  #1  header  (32  x  5  MHz) 

/*  base  address  for  mta  #1  data  */ 

I*  base  address  for  mta  #2  header  (16  x  2.25) 

/*  base  address  for  mta  #2  data  */ 

/*  base  address  for  mta  #3  header  (16  x  2.25) 

/*  base  address  for  mta  #3  data  */ 

/*  base  address  for  time  header  */ 

/*  base  address  for  time  data  "/ 

/*  base  address  for  dummy  header  */ 

/*  base  address  for  dummy  data  */ 

/*  storage  base  for  Fl  rms  profiles  */ 

/*  storage  base  for  F2  rms  profiles  */ 

/*  storage  base  for  F3  rms  profiles  */ 

/*  storage  base  for  SPI  mean  profiles  */ 

/*  sampling  mode  -  init.  to  3  frequency  mode  * 

/*  number  of  SPI ' s  to  sample  */ 

/*  permanent  npp  storage  */ 

/*  number  of  points  in  profile  */ 
iod;  /*  permanent  prof_period  storage  */ 

/*  final  profile  storage  period  in  hundredths 

/*  permanent  npa  storage  */ 

/*  number  of  profiles  to  average  */ 

/*  permanent  nupc  storage  */ 

/*  number  of  unaveraged  profiles  remaining  */ 

/*  number  final  data  file  profiles  to  collect 

/*  total  number  of  shore  profiles  */ 

/*  sometimes  shore  profiles,  sometimes  1  */ 

/*  number  of  profiles  data  file  will  hold  */ 

/*  dif  between  shore  and  mem  profile  numbers  * 
/*  storage  flag  -  zero  while  storing  */ 
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unsigned  long  perm_sdt;     /* 

unsigned  long  nfpc;         /* 

unsigned  long  wordprof;     I* 

unsigned  long  spichans;     /* 

unsigned  long  offset;       /* 

unsigned  long  hours_between; 

unsigned  long  num_collections, 

unsigned  short  bytecomb;    /* 

unsigned  char  *offstart;    /* 

unsigned  char  *offend;      /* 

unsigned  long  offlen;      /* 

unsigned  long  pad;  /* 

unsigned  short  vreg;       /* 

unsigned  short  vunreg,-     /* 

unsigned  short  vbus;        /* 

unsigned  short  volt_leakl;  /* 

unsigned  short  volt_leak2r  /* 

unsigned  short  tempin;     /* 

unsigned  short  max_files 

unsigned  long  start_time 

unsigned  short  mta_blank 

unsigned  short  mta_delay 

unsigned  short  mta_gain; 

unsigned  long  mta_numscans;    /* 

unsigned  long  shore_mta_numscans; 

unsigned  long  temp_mta_numscans,- 

unsigned  short  mta_pulse;    /*  mta  acoustic  pulse  width  microsecs 

unsigned  short  mta_thresh[4] ;   /*  mta  threshold  in  millivolts  */ 

unsigned  char  prof_index,-     /*  shore  data  header  index  */ 


Number  of  hours  to  collect  shore  data  */ 
number  of  final  profiles  remaining  */ 
number  of  bytes  in  a  profile  of  words  *7 
number  of  bytes  in  all  SPI ' s  */ 

.01's  of  hours  from  12  AM  for  1st  file-timed  */ 
/*  ,01's  of  hours  between  df ' s  -  timed  mode  */ 
/*  total  data  files  to  collect  in  timed  */ 
location  to  combine  A-D's  lsb  &  msb  */ 
starting  address  of  offload  */ 
ending  address  of  offload  */ 
length  of  file  to  offload  */ 
number  of  pad  values  to  use  in  dummy  */ 
location  for  regulated  voltage  sampling  */ 
location  for  unregulated  voltage  sampling  */ 
location  for  bus  voltage  sampling  */ 
location  for  leak  detector  #1  sampling  */ 
location  for  leak  detector  #2  sampling  */ 
internal  temperature  sample  */ 
maximum  number  of  data  files  on  disk.  */ 
starting  time  of  data  acquisition  */ 
mta  blanking  interval  in  microsecs  */ 
time  between  mta  scans  in  0.01's  of  sees  */ 
mta  gain  */ 
/*  number  of  mta  scans  to  make  */ 

number  of  mta  scans  -shore  collection  */ 


*  Function  declarations 

int  menu_input { int , int , int ) ; 
void  an_power (int) ; 
void  clear_mem(void) -, 
void  collect_mta (void) ; 
void  data_f ile_res (void) ; 
void  def_samp_scheme (void) ; 
void  delay_no_clock(void) ; 
void  display_n(int, int, int) ; 
void  f lexc (void) ; 
void  f 2exc (void) ; 
void  f 3exc (void) ; 
void  grabdata (char) ; 
void  grabstring (char) ; 
void  initialize (void) ; 
void  inst_test  (void)  ,- 
void  main_sample (void) ; 
void  mat_headers (void) ; 
void  mem_setup (void) ; 
void  nothing (void) ; 
void  of f load(void) ; 
void  ret_f rom_reset (void) ; 
void  set„time (void) ; 
void  start_collection(void) ; 
void  5tart_mta (void) ; 
void  storm_mode (void) ; 
void  talk_thru (void) ; 
void  tpu_setup (void) ; 
void  voltage (void) ; 


Main  program 


/*  variable  declarations  */ 
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bool  reset; 
char  m,ack; 

float  fvreg, fvunreg, fvbus; 
int  i; 
time_t  timer; 

program  */ 

/*  reset  check  goes  here  */ 

if  (reset==l)  ret„f rom_reset ( ) ; 

initialize ( ) ; 

in_menu : 

/*  Get  voltages  from  charging  circuitry  */ 
voltage ( ) ; 
fvunreg  -  VUNREGCAL  *  ( float) (vunreg) ; 
fvreg  =   VREGCAL  *  { float) (vreg) ; 
fvbus  =  VBUSCAL  *  (float) (vbus) ; 
if  (  (storm==l)  &&  (fvunreg<5  .0)  )  storm_mode  ( )  ,- 
if  (storm==l)  voltage ( ) ; 
printf ("\n\nMonsterLog  1.00.005  Next  data  file  -  %03dH 

, data_f ile_pointer) ; 
timer=RtcToCtm() ; 
printf (" \nSampling  mode  is  "); 
if  (samp_mode==0) 

printf ( "CONTINUOUS    " ) ; 
else 

printf ("TIMED         " ) ; 
printf ( "  4s" , ctime l&timer) ) ; 

printf ( "Charging  is  " ) ; 
if  <e5==0> 

printf ("ON  -  Full  charge  ■); 

else 

printf  ("OFF  "  )  ,- 

if  ( (volt_leakl>250Q)6&(volt_leak2>2500) ) 

printf ("No  leak  detected\n" ) ; 
else 

( 

printf (FLSHON); 
printf ( "LEAK  DETECTED\n" ) ; 
printf (FLSHOFF) ; 
} 
printf ("Package  temperature  =  %04 . If ",[ tempin*50Q. 0/4096 . 0) -273 . 0) ; 
printf!"  Emergency  storm  mode  is  "); 

if  {storm==0) 

printf ("DISABLED\n") ; 
else 

printf ( "ENABLED\n" } ; 
printf ("Voltages:  Shore  =  %05.2f   Battery  =  %05.2f   Bus  =  %05.2f\n\n" 
, fvunreg, fvreg, fvbus) ; 


printf ("A  -  Analog  power  menu\n"); 
printf ("D  -  Define  sampling  scheme\n"); 
printf  ("E  -  Toggle  storm  mode\n" )  ,- 
printf ( "I  -  Instrument  tests \n" ) ; 
printf ("M  -  MTA  talk~through\n" ) ; 
printf ("0  -  Offload  data  files\n"); 
printf ( "Q  -  Quit\n") ; 

printf ("R  -  Reset  next  data  file  number\n"l; 
printf ("S  -  Start  data  collectionSn" ) ; 
printf ("T  -  Set  time  and  date\n"); 
printf ( "ENTER  -  Redisplay  this  menu\n\n" ) ; 
printf ( " 

printf ( "Selection  =>"); 
selection: 

m- (char )menu_input (0,0,0) ; 
if  (m==(-l)) 

{ 

printf  (CLEAR)  ,- 

printf ( "Timeout\x07" ) ; 

goto  main_menu; 
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} 
if   (ro«lQ) 
{ 

printf (CLEAR) ; 

goto  main_menu; 

} 
if     (!  t»=='S'  |  |ra=='I'  |  |m=='R'  |  |m=='T'  |  |m=='A'  |  |m==*M' 
| |m=='0' 1 |m=='Q' | |m==-D' | |m=='E' ) ) 

goto  selection; 
printf("%c   RETURN  to  accept,  any  other  key  to  reject .%c" ,m, 13) ; 
ack=(char)menu_input (0, 0,0); 
if  (ukM(-l)) 

{ 

printf ! " \n\nTimeout\x07\n\n" ) ; 

goto  main_menu; 

) 
if  (ack!=10) 

{ 

printf (" 

,13)  j 

printf ( "Selection  =>"); 

goto  selection; 

) 
printf ( "\n" ) ; 

if    (m--'Q')    goto   main_bail; 
if    (m== ' I ' )    inst_test ( ) ; 
if    (m== ' A' )    an_power (0) ; 
if    (m== ' D' )    def_samp_scheme ( ) ; 
if    (m== 'R'l    data_file_res ( ) ; 
if    (m== 'T ' )    set_time ( ) ; 
if    lm=='S' )    start_collection( ) ; 
if    (m=='0'l    offloadd; 
if    (m=='M')    talk_thru(); 
if    (m=='E'l 

{ 

if  (storra==0)  storm=l; 

else  storm-0; 

printf (CLEAR) ; 

) 
goto  raain_menu; 
ma  in_ba  i 1 : 

printf ( "VnDropping  into  monitor\nH ) ; 


*  Analog  power  routine 


void  an_power(int  a) 
{ 

char  m; 
float  fvreg, fvunreg, fvbus; 

power_menu: 

if  (a==0) 
{ 

printf (CLEAR) ; 
printf (UNDON) ; 

printf ( "Power  selection\n\n" ) ; 
printf (UNDOFF) f 

voltage ( ) ; 
fvunreg  -  VUNREGCAL  *  ( float) (vunreg) ; 
fvreg  -   VREGCAL  *  ( float) {vreg) ; 
fvbus  =  VBUSCAL  *  (float) (vbus) ; 


printf ("Voltages:  Shore  =  %05.2f   Battery  ■  %05.2f   Bus  ■  %05.2f\n\n" 
,  fvunreg,  fvreg,  fvbus)  ,- 

printf  ( "C  -  Charging " )  ; 

if  (e5==0) 
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printf ( "ON\n") ; 
else 

printf ("OFF\n") ; 
printf ("A  -  ACP  |3  frequency)..."); 
if  (f6«0) 

printf {"OFF\n") [ 
else 

printf ("ON\n") ; 

printf  ("0  -  OBS/pressure "); 

if  (e7==0) 

printf ("OFF\n") ; 
else 

printf <"ON\n") j 

printf ("1  -  ADV  #1 "); 

if  (eO==0) 

printf {"OFF\n") ; 
else 

printf ("ON\n") ; 

printf { "2  -  ADV  #2 " )  ; 

if  (fl»0) 

printf ("OFF\n"> ; 
else 

printf { "ON\n") ; 
printf {"T  -  Tilt,  compass,  temp."); 
if  (e6==0) 

printf ( "OFFSn" ) ; 
else 

printf ("ON\n") ; 

printf  {"M  -  MTA's "); 

if  (df7==l) 
{ 

printf  ("OFF  -   Disabled\n" )  ,- 
} 

else 
{ 

if  (f7==0) 

printf ("OFF  -  Enabled\n"); 
else 

printf ( "ON\n" ) ; 
} 

printf  ("V  -   Video  camera "); 

if  (el==0) 

printf ("OFF\n") ; 
else 

printf ( "ON\n" } ; 
printf ( "ENTER  -  return  to  main  menu"); 
m=  { char )  menu_input  (0,0,0)  ,- 
if  (m==(-l))  goto  power_bail; 
} 
else 

m=0; 

if  (m=='C")  /*  Charging  on/off  -  unaffected  by  'a'  value 

( 
if  (e5--Q)  /*  Turn  charge  off  */ 

{ 

e5  =  l; 

PSet(E,5)  ,- 

} 
else 

t 

e5  =  0; 

PClear(E,5)  ,- 

} 
goto  power_menu; 
} 
if  ( (m=='A' ) | | (a!=0) )    f*    ACP  (3  frequency)  on/off  */ 
{ 
if  ( ( (f6==0)&&(a==0) ) | | (a==l) )     (*   ACP  on  */ 

{ 

f6  =  l; 
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PSet(F,6) ; 

} 
else 

t 

£6  =  0; 

PClear  {F,  6) ; 

) 
if  (a==0)  goto  power_menu; 
) 
if  { (m=='0' ) | | (a!=0) )    /*  OBS/pressure  on/off  V 
( 
if  ((  (e7==0)s,Ma==0))  I  I  <a==l) )     /•  OBS/pressure  on  */ 

( 

e7  =  l; 

PSet(E,7) ,- 

) 
else 

{ 

e7  =  0; 

PClear(E,7)  ; 

) 
if  (a==0)  goto  power_menu; 
) 
if  llm'l'!  |  |  (a!  =0)1    /•  ADV  #1  on/off  */ 
( 
if  ( ( (e0==0)&&{a==0) ) | | (a==lt )     /*  ADV  #1  on  «/ 

{ 

e0  =  l; 

PSet(E,0) | 

) 
else 

( 

e0  =  0; 

PClear (E,0) j 

) 
if  (a==0)  goto  power„menu; 
) 
if  IIVS'l  |  |  (a!=0))    /•  ADV  »2  on/off  */ 
( 
if  ( ((fl==0)ss(a==0) ) | | (a==l) )     /*  ADV  »2  on  */ 

( 

fl  =  l; 

PSet(F.l)  ; 

) 
else 

( 

£1  =  0; 

PClear  (F,  1)  ; 

} 
if  (a==0)  goto  power_menu ; 
) 
if  { (m=='T' ) | | (a!=0) )    /•  Tilt,  compass,  temp  on/off  */ 
< 
if  ( ( (e6==0)&t(a==0) ) | | (a==l) }      /*  Tilt,  compass,  temp  on  * 

{ 

e6  =  l; 

PSet(E,6)  ; 

} 
else 

( 

e6=0; 

PClear  (E,  6)  ; 

) 
if  (a==0)  goto  power_menu; 
) 
if  (  (m=='M')  |  |  (a!=0))  /•  MTA  on/off  */ 

( 
if  (  (  (a==0)&Md£7==l)&s(f7«0)  )  |  |  (  (a==l)  Si  (d£7==0)  )  ) 

{  /*  In  menu,  disabled,  off  */ 

f7  =  l;      /*  Turn  on  */  /*  OR  function  call  on,  */ 

PSet(F,7) ;  /»  enabled  '/ 


106 


df7=0; 


/*  Enable  */ 


else  if  ({a==0)&&(df7==0)&&{f7==0)) 
{ 

df7=l;      /*  Disable  */ 
} 
else  if  ((a==0)fcfc(f7==l)) 
{ 

f7=0;       /*  Turn  off  */ 
PClear(F,7) ; 
df7=0;      f*    Enable  */ 
} 
if  (a==0)  goto  power_menu; 
} 
if  (am'V) 
{ 
if  (el==0) 


/*  In  menu,  enabled,  off 


In  menu,  on 


/*  Video  on/off  -  unaffected  by  'a'  value 


el=l; 

PSet(E.l) ; 
} 


el=0; 

PClear (E,l) ; 
} 

goto  power_menu; 

} 
if  ( (m!=10)&&(a==0) ) 

{ 

printf (CLEAR) j 

goto  power_menu; 

} 
power_bail : 
if  U(m==10)  |  |  (m==(-l)  >)&&(a==0)>  printf  (CLEAR)  ; 


Clear  memory  &  rms  lookup  table  initialization  routine 


void  clear_mem(void) 
{ 
int  i  ; 


machine  mc68020    ;Allows  advanced  addressing  modes 
movem.l  d0-d7/a0-a6, - (a7)   ;Push  registers  on  stack 


,-Data  memory 

move . 1 

_f lsum, aO 

move . 1 

_f lsqr, al 

move . 1 

_f2sum,a2 

move . 1 

_f2sqr, a3 

move . 1 

_f 3sum, a4 

move . 1 

_t  3sqr ,  a5 

move . 1 

_npp , d3 

subq. 1 

§$01,d3 

iloop  move.l  #0, (a0}+ 

; clear  eac 

move . 1 

#0, (al)+ 

move . 1 

#0, (a2)+ 

move . 1 

#0,  <a3)  + 

move . 1 

#0, (a4)+ 

move . 1 

#0,  [a5)  + 

dbf  d3 

iloop:         , 

move . 1 

_spsum, aO 

move. 1 

_spsqr, al 

move . 1 

_nosp,d3 

subq . 1 

#$01, d3 

jloop  move. 1  #0, (aO) + 

; clear  eac 

x)  memory 
xA2)  memory 
x)  memory 
x'2)  memory 
x)  memory 
x"2 )  memory 


use  aO  for  address  of  Fl  sum 
use  al  for  address  of  Fl  sum 
use  a2  for  address  of  F2  sum 
use  a3  for  address  of  F2  sum 
use  a4  for  address  of  F3  sum 
use  a5  for  address  of  F3  sum 
use  d3  for  index 
reduce  by  1  for  proper  indexing 
each  memory  location 


continue  until  all  memory  clear 
use  aO  for  address  of  SPI  sum(x)  memory 
use  al  for  address  of  SPI  sum(x"2)  memory 
use  d3  for  index 
reduce  by  1  for  proper  indexing 
each  memory  location 
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movG.l  #0, (al)  + 

dbf  d3,jloop:        .-continue  until  all  memory  clear 

movem.l  (a7) +, d0-d7/a0-a6   ;Pull  registers  from  stack 
#endasm 

/*  xA2  lookup  table  */ 
#asm 

movem.l  d0-d7/a0-a6, - (a7)   ;Push  registers  on  stack 

;  x-^2  table 

move.l  _sqbase,a0  ;use  aO  for  base  of  x"2   lookup  table 

move.l  #$0FFF,d0  ;start  with  highest  value 

s21oo  move.l  dO.dl  ;transfer  value  to  dl  for  squaring 

mulu  d0,dl  ; square  it 

move.l  dl, ($0000, a0,d0. 1*4)    ;store  it  using  root  as  index 

dbf  d0,s21oo:  ,-continue  for  entire  table 

movem.l  (a7) +,d0-d7/a0-a6      ;Pull  registers  from  stack 


(x+0 . 5) "2   lookup  table  initialization  */ 

C  (i=0;i<4096;i++) 
{ 

* (rtbase+i)= (unsigned  long) ( (float) (i+0.5)*(i+0.5) *npa) ; 
) 


*  Collect  MTA  data 

void  collect_mta (void) 
{ 

char  mta_id; 
char  comp_string[10]  ; 
int  i , j ; 
long  timeout; 

/*  Initialize  TPU  communications  */ 

PClear(E,2) ;  /*  Enable  TPU  buffers  */ 

i=TSerOpen(inchan, INP, 100, 9600, 'N' , 8, 1) ;     /*  input  channel  */ 
printf ( " \nOpen  channel  %d  -  status  =  %d,  INP  =  %d\n" , inchan, i, INP) ; 
i=TSerOpen(outchan,OUTP, 100, 9600, 'N' , 8, 1) ;   /*  output  channel  */ 
printf ("Open  channel  %d  -  status  ■  %d  OUTP  ■  %d\n\n" , outchan, i,OUTP) ; 

for  <j=0;j<3;j++) 
( 

/*  Collect  MTA  data  -  R2  command,  MTA  #  mta_id  */ 

mta_id= (char)  j+49; 

request_data: 

printf ( "Sending  data  request  (%cR2)  .  .  .  \n"  ,mta_id)  ,- 

TSerPutByte  ( outchan ,  0 )  ,- 
TSerPutByte(outchan,mta_id)  ,- 
TSerPutByte  (outchan,  'R'  )  ,- 
TSerPutByte (outchan, '2'); 
TSerPutByte ( outchan , 13 ) ; 

printf! "Awaiting  acknowledge  (%cA2) . . . \n" ,mta_id) ; 
grabstring(mta_id) ; 
if  (stringerror ! =0) 


108 


goto  collect_mta_bail; 
sprintf (comp_string, "%cA2 " ,mta_id) ; 
if  (strcmp(instring,  comp_string)  !  =0) 

sprintf (comp_string, "%cN" ,mta_id) ; 
if  (strcmp (instring, comp_string) ==0) 

( 

printf  ( "Device  not  ready.  ..  \n" )  ,- 

timeout=0; 

while (timeout<200000) 
timeout+=l; 

goto  request_data; 

} 
goto  collect_mta_bail; 
} 

printf ( "Acknowledged  -  beginning  receipt  of  data...\n"); 

grabdata(mta_id) ; 
if  (stringerror ! =0) 

goto  collect_rata_bail; 

printf ("Data  received  &  residing  at"); 
if  (j==0)  printf  ("  %p\n"  ,mtlbas)  ,- 
if  (j==l)  printf("  %p\n" ,rat2bas) 
if  (j==2)  printf!"  %p\n" ,mt3bas) 

ollect_mta_bail ■ 

if  (stringerror ! -0) 

printf ( "Communications  error. \n\n" ) ; 


TSerClose(inchan) ; 

TSerClose(outchan) ; 

PSet(E,2);        /*  Disable  TPU  buffers 


Data  file  reset  routine 


void  data_f ile_res (void) 

{ 

int  n,m; 

printf (CLEAR) ; 

printf (UNDON) ; 

printf ( "Reset  data  file  pointer"); 

printf (UNDOFF) | 

printf ("\n\nEnter  new  value  for  disk  file  pointer  (l-%d).\n" 

,max_f iles)  ; 
get_disk_pointer : 

n-data_f ile_po inter; 

m=inenu_input  (3  ,  0,n)  ; 

if    (m--(-l))    goto  data_bail; 

if    ((m<l) | | (m>500))    goto   get_disk_pointer; 

data_f ile_pointer=m; 


data_bail : 
printf (CLEAR) ; 


*  Sampling  scheme  definition  routine 
void  def„samp_scheme (void) 
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char  am; 

int  m,n; 

int  min_hours; 

int  total_time_hours; 

int  total_time_mins; 

int  total_time_secs; 

int  total_time_hunds; 

float  total_secs; 

time„t  timer; 

struct  tm  time_struct, *vtime_struct; 

printf (CLEAR) ; 

printf (UNDON) ; 

printf ( "Define  sampling  scheme"}; 

printf (UNDOFF) ; 

printf (" \n\nNumber  of  points  per  ACP  {3  freq. )  prof ile. \n" ) ; 
get_npp: 

n=perm_npp; 

m=menu_input ( 3 , 0 ,  n) ; 

if  (a«*(-l))  goto  scheme_bail; 

if  ( (m<10) | | (m>200) )  goto  get_npp; 

perm_npp-m; 

printf (" \n\nAveraged  profile  &  SPI  sample  period  in  seconds. \n" ) ; 
get_period: 

n=perm_prof_period; 

m=menu_input ( 3 , 1 , n ) ; 

if  (m==(-l))  goto  scheme_bail; 

if  ( (m<3) | | (m>200) )  goto  get_period; 

perm_prof_period=m; 

printf ( N\n\nSeconds  between  MTA  profiles  (minimum  2.00).\n"); 
get_mta_del: 

n-mta_delay,- 

m-menu_input { 4 , 1 , n ) ; 

if  (%»■(-!))  goto  scheme_bail; 

if  ( (m<100) | | (m>6000) )  goto  get_mta_del; 

mta_delay=m; 

npp-perm_npp  ,- 
nfpc=0; 
mem_setup ( ) ; 

printf ( "\n\nNumber  of  profiles  for  TT  data  file  -  default  is  maximum. \n" ) ; 
get_nf pc : 

n=nfp; 

m=menu_input (4, 0, n) ; 

if  (m--(-l))  goto  scheme_bail; 

if  ( (m<l) | | (m>nfp) )  goto  get_nfpc; 

shore_nfp-ra; 

dif_nfp=l; 

if  (m-=nfp) 

{ 

mem_nfp-0; 

} 
else 

{ 

mem_nfp=m; 


printf ( "\n\nNumber  of  profiles  to  include  in  average"), - 
printf ( "  -  default  is  maximum. \n" ) ; 
get_npa : 

n=perra_prof_period-3; 
m-menu_input ( 3 , 0 , n) ; 
if  (m==(-l))  goto  scheme_bail; 

if  ( (m<l) | | (m>(perm_prof_period-3) ) )  goto  get_npa; 
perm_npa=m,- 
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npp-pernwipp ; 

pro f _per i od-perm_pro f _per iod ; 

npa=perm_npa  ,- 

nfpc=mem_nfp,- 

mem_setup( ) ; 

total_secs=nfp*prof_period/10Q . 0; 

total_time_hours=total_secs/3600,- 

total_time_mins= { total_secs-3600 . O*total_time_hours ) / 60 . 0 ; 

total_time_secs= { total_secs-3600 . O*total_time_hours 

-60 . O*total„time_mins) ; 
total_time_hunds= (total_secs- (int) total_secs) * 100.0; 
min_hours=total_time_hours+l; 
sprintf ( acquis! tion_time, "%02d: %02d:%02d. %02d" , total_time_hours 

, total_time_mins, total_time_secs, total_time_hunds) ; 

printf (CLEAR) ; 

printf ( "Total  sampling  time  ■  %s" , acquisition_time) ; 

printf ( "\n\nTotal  hours  acquisition  for  shore  data  file.  (0-TT  file)\n"); 
get_sdt : 

n-perm_sdt; 
m=menu_input ( 3 ,  0 , n) ; 
if  (m==-l)  goto  scheme_bail ; 
if  ( (m<0) | | (m>48) )  goto  get_sdt; 
perm_sdt-m; 

shore_nfp- (100 . 0/perm_j?rof_period) *3600*perm_sdt; 
if  (shore_nfp==0) 
{ 
shore_nfp=nfp; 

dif_nip=l; 
) 
else 

dif_nfp=shore_nfp-nfp+l  ; 

i  f    (penn„sdt>rain_hours ) 

min_hours=perm_sdt; 

printf C'\n\nSet  MTA  parameters?  (Y/N) \n" ) ; 
printf  ( "'Selection  ->"}; 
get_mtapar : 

am=menu„input (0,0,0); 

if  (am==(-l))  goto  scheme_bail; 

if  ( ! ( (am==(  Y' ) ) | | (am==( 'N' ) } ) )  goto  get_mtapar; 

printf("%c   ENTER  to  accept,  any  other  key  to  reject .", {char Jam) ; 

m-menu_input (0,0,0) ; 

if  (m==(-l))  goto  scheme_bail; 

if  (m!=10) 
{ 
printf ( "%c 

,13,13) ; 
printf ( "Selection  ->"); 
goto  get_mtapar; 
} 
if  (aapuB'N1) 

goto  nomta,- 

setmta: 

printf ( "\n\nMTA  blanking  in  millimeters  (all  MTA's)\n"); 
get_mta_blank: 

n-mta_blank; 

m-menu_input (3 , 0, n) ; 

if  (m=={-l)}  goto  scheme_bail; 

if  ( (m<0) | | (m>325) )  goto  get_mta_blank; 

mta_blank=m; 

printf ("\n\nMTA  #1  threshold  in  millivolts.  \n" )  ,- 
get_mta_threshl : 
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n=mta_thresh[l]  ,- 

m=menu_input (4, 0,n) ; 

if  (m==(-l)}  goto  scheme_bail; 

if  ((m<0) | | (m>9999))  goto  get_mta_threshl; 

mta„thresh[l] =m; 

printf  ( "\n\nMTA  #2  threshold  in  millivolts.  \n" }  ,- 
get_mta_thresh2 : 

n=mta_thresh[2] ; 

m=mGnu_input (4, 0,  n)  ; 

if  (m==(-l))  goto  scheme_bail; 

if  ( (m<0) | | (m>9999) )  goto  get_rata_thresh2; 

mta„thresh[2] =ra; 

printf  ( "  \n\nMTA  #3  threshold  in  millivolts.  \n" )  ,- 
get_mta_thresh3 : 

n=mta_thresh[3] ; 

m=menu_input (4 , 0,n) ; 

if  (m==(-l))  goto  scheme_bail; 

if  ({m<0) | | (m>9999))  goto  get_mta_thresh3 ; 

mta_thresh [ 3 ] -m; 

nomta : 

printf (" \n\nSampling  mode  -  currently  "); 
i  f  ( samp_mode==0 ) 

printf ("COKTINUOUSSn") ; 
else 

printf ("TIMED\n") ; 
printf ( "   C  -  ContinuousSn" ) ; 
printf("   T  -  Timed\n"); 
printf!"  %c 

,13)  ; 
printf ! "Selection  ->"); 
get_mode : 

am=menu_input (0,0,0)  ; 

if  (am==(-l))  goto  scheme_bail; 

if  (! ( (am==( 'C ) ) | | <am==( 'T' ) ) ) )  goto  get_mode; 

printf("%c   ENTER  to  accept,  any  other  key  to  reject .",  (char)  am)  ,- 

m=menu_input (0,0,0); 

if  (m==(-l})  goto  scheme_bail; 

if  (m!=10) 
{ 
printf ( "%c 

,13,13) ; 
printf ( "Selection  =>"}; 
goto  get_mode; 
} 
if  (am=='C') 

samp_mode=0; 


else 


{ 

samp_mode=l; 

f irst_time-l; 

timer=RtcToCtm ( ) ; 

vtime„struct=localtime (&timer) ; 

time_struct .  tm_year-vtime_struct->tm_year,- 

time_struct. tm_mon-vtime_struct->tm_mon; 

time_struct .  tm_n\day=vtime_struct->tm_mday; 

time__struct.  tm_hour=0; 

time_struct . tm_min=0 ; 

time_struct . tm_sec-0; 

start_time=mktime(&time_struct)+( (RtcToCtmO -mktime (&time_struct) ) 

/ (3600*hours_between) } * (3600*hours_between) 

+3600*hours_between,- 


if  (samp_mode--0)  goto  scheme_bail; 

printf ( "\n\nHours  between  collection  starting  times. \n" ) ; 
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get_between: 

n=hours_between ; 

m=mGnu_input ( 2 , 0 , n) ; 

if  (m==(-U)  goto  scheme_bail; 

if  (  (in<min_hours)  |  |  (m>24)  )  goto  get_between,- 

hours_between=m; 

timer=RtcToCtm( )  j 

vtime_struct=localtime(&timer) ; 

time_struct. tm_year=vtime_struct->tm_year; 

t  ime_s  true  t .  tm_mon=vt  ime_s  true  t  -  >  tm_mon  ; 

tirae_struct.  tm_mday=vtime_struct->tm_mday; 

time_struct. tm_hour=0 ; 

time_struct. tm_min-0; 

time_struct. tm_sec-0; 

start_time=mktime(&time_struct)  +  (  (RtcToCtmO  -mktime (&time_struct)  ) 

/ (3  600*hours_between) ) * (3600*hours_between) 

+3600*hours_between; 


printf ( " \n\nNumber  of  collections  (1  data  file  per  collection) \n" ) ; 
get_collect: 

n-num_col lections; 

m=menu_input ( 3 , 0 , n ) ; 

if  (m=={-D!  goto  scheme_bail; 

if  ( (m<l) | | (m>(max_files-data_file_pointer+l) ) )  goto  get_collect; 

num_col lections =m; 

scheme_bail: 

printf (CLEAR) ; 
} 


Delay  for  ~1  second  without  using  clo 


void  delay_no_clock(void) 
{ 
long  waste_time; 

waste_time=500000; 
while (waste_time>0) 

waste_time-=: 
} 


Display  value  of  numeric  input  for  menu_input ( ) 


void  display_n ( int  a,int  d, int  n) 
{ 
float  f; 

if  (d»*lj  f=(float)n/100; 
printf ("%c" ,13) ; 
switch (a) 
{ 
case  1 : 

if  (d==0) 

printf ( "=>%d",n) ; 
else 

printf ("=>%. 2 f",f) ; 
break; 
case  2: 

if  (d==0) 

printf ("=>%02d",n) ; 
else 

printf ("=>%. 2f", f ) ; 
break; 
case  3 : 
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if  (d==0) 

printf ("=>%03d",n) ; 
else 

printf ("=>%04.2f",f) ; 
break; 
case  4 : 

if  (d==0) 

printf ("=>%04d" ,n) ; 
else 

printf ("=>%05. 2 f",f); 
break ; 
} 


Fl  exception 


/*  flexc  is  the  exception  which  executes  on  an  interrupt  from  TPU  11.  1 
transducer  is  sampled  in  this  routine.  */ 

void  flexc (void) 

{ 
#asm 

; f lentry  is  the  entry  point  for  Fl ' s  sampling  routine  exception 
flentry  movem.l  d0-d7/a0-a6,  -  (a7)      ,-Push  registers  to  stack 

; Initialize  address  registers  for  profile  sampling 


move.l  _flsum,. 
move.l  _flsqr,al 
lea  _bytecomb,a4 
move.l  _sqbase,a5 
move.l  „npp,d3 
subq.l  #$01, d3 


use  aO  for  address  of  Fl ' s  sum(x) 
use  al  for  address  of  Fl ' s  sum{xA2) 
location  for  byte  combining  in  a4 
use  a5  for  base  of  x"2  lookup  table 
use  d3  for  index  (points  in  profile) 
reduce  by  one  for  proper  indexing 


f  lwpa 

move.b  $F8COOA,dO 

bp 1  f lwpa ; 

f  lwpb 

move.b  $F8CO0A,dO 

bmi  flwpb: 

;At  this  point,  both  the  channel  number  and  port  number  are 

;Set  the  multiplexer  and  as  soon  as  the  A->D's  not  busy,  set  the 
; channel  number  for  the  next  point  collected.  This  also  provides 
;time  for  the  MUX  to  settle  down  after  switching  ports. 

bclr  #$03, $FFFFFC15     ;Set  port  0  on  multiplexer  board 
bclr  #$04,$FFFFFC15 
bclr  #$05,$FFFFFC15 

Read  SFR  to  check  busy  flag 

;Stay  here  until  busy  flag  goes  low 

SFR  -  busy  goes  high  when  conversions  done 

;Stay  here  until  conversion  complete 
move.b  #O,$F8C0O0        ;Set  channel  0  for  next  sample 

,-Channel  0  and  port  0  are  now  set  for  the  next  sample 

flloopa   move.b  $F8C0OA,dO        ;Read  SFR  to  check  busy  flag 

bpl  flloopa:  ;Stay  here  until  busy  flag  goes  low 

flloopb  move.b  $F8C0OA,d0       ;SFR  -  busy  goes  high  when  conversions  done 

bmi  flloopb:  ;Stay  here  until  conversion  complete 

move.b  $F8C002, (a4}+     ;Get  MSB  from  A->D 

move.b  $F8C0OO, (a4) +     ;Get  LSB  from  A->D 

move.b  #0,$F8C0OO       ; Set  channel  0  for  next  sample 

move.w  -(a4),dl         ;Bring  entire  word  into  dl 

add.l  dl,(a0)+  ;Adds  to  current  sum  for  profile  point 

move.l  ($0000,a5,dl.l*4t ,d2    ; Find  square  of  dl  from  lookup  table 

add.l  d2,(al)+  ;Adds  to  current  sum(x'v2)  for  profile  point 

dbf  d3, flloopa:         ;Jump  back  to  main  loop  until  profile's 


flbail   eori.w  #$0800,  $FFFE20   ,-Acknowledges  and  turns  off  interrupt 

movem.l  (a7) +,d0-d7/a0-a6   ;Pull  registers  from  stack 
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F2  exception 


f2exc  is  the  exception  which  executes  on  an  interrupt  from  TPU  10. 
F2  transducer  is  sampled  in  this  routine.  */ 


void  f2exc(void) 


;f2entry  is  the  entry  point  for  F2 ' s  sampling  routine  exception 
f2entry  movem.l  d0-d7/a0-a6,  -  (a7)      ,-  Push  registers  to  stack 


; Initialize  address  registers  for  profile  sampling 


move . 1  _f2sum,a0 
move.l  _f2sqr,al 
lea  _bytecomb,a4 
move.l  _sqbase,a5 
move . 1  „npp , d3 
sufaq.l  #501, d3 


use  aO  for  address  of  F2 ' s  sum(x) 
use  al  for  address  of  F2 ' s  sum(x"2) 
location  for  byte  combining  in  a4 
use  a5  for  base  of  x"2  lookup  table 
use  d3  for  index  (points  in  profile) 
reduce  by  1  for  proper  indexing 


;At  this  point,  both  the  channel  number  and  port  number  are 

;Set  the  multiplexer  and  as  soon  as  the  A->D"s  not  busy,  set  the 
,-channel  number  for  the  next  point  collected.  This  also  provides 
;time  for  the  MUX  to  settle  down  after  switching  ports. 


bclr  #$03,$FFFFFC15 
bclr  #$04,$FFFFFC15 
bclr  #$05,$FFFFFC15 


;Set  port  0  on  multiplexer  board 


f2wpa    move.b  SF8C00A,d0 

bpl  f2wpa: 
f2wpb    move.b  SF8CO0A,d0 

f2wpb: 


move.b  #O,$F8C000 


,-Read  SFR  to  check  busy  flag 

,-Stay  here  until  busy  flag  goes  low 
,-SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 


; Set  channel  0  for  next  sample 


; Channel  0  and  port  0  are  now  set  for  the  next  sample 


f 21oopa  move. 
f21oopb  move. 


$F8C0OA,dO        ;Read  SFR  to  check  busy  flag 

bpl  f21oopa:  ;Stay  here  until  busy  flag  goes  low 

$F8C00A,d0       ;SFR  -  busy  goes  high  when  conversions  done 


bmi  f21oopb: 
move.b  $F8C002, (a4)+ 
move.b  $F8COO0, (a4)+ 
move.b  #0,$F8COOO 
move.w  - (a4) , dl 
add.l  dl, (a0)+ 
move.l  ($0000,a5,dl. 
add.l  d2, (al)+ 
dbf  d3, f21oopa: 


Stay  here  until  conversion  complete 

Get  MSB  from  A->D 

Get  LSB  from  A->D 

Set  channel  0  for  next  sample 

Bring  entire  word  into  dl 

Adds  to  current  sum  for  profile  point 

d2    ;Find  square  of  dl  from  lookup  table 

Adds  to  current  sum(xA2)  for  profile  point 

Jump  back  to  main  loop  until  profile's 


*  F3  exception 


eori.w  #$0400, $FFFE20   ; Acknowledges  and  turns  off  interrupt 
movem.l  (a7) +, d0-d7/a0-a6   ;Pull  registers  from  stack 


/*  f3exc  is  the  exception  which  executes  on  an  interrupt  from  TPU  9.   The 

F3  transducer  is  sampled  in  this  routine,  and  the  number  of  unaveraged 
profiles  remaining  to  collect  is  determined.  */ 
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void  f3exc(void) 


,-f3entry  is  the  entry  point  for  F3  '  s  sampling  routine  exception 
f3entry  movem.l  d0-d7/a0-a6, - (a7)      ; Push  registers  to  stack 


; Initialize  address  registers  for  profile  sampling 


move.l  _f3sum,a0 
move. 1  _f 3sqr , al 
lea  _bytecomb,a4 
move.l  _sqbase,a5 
move . 1  _npp , d3 
sutaq.l  #$01,(13 


use  aO  for  address  of  F3 ' s  sum(x) 
use  al  for  address  of  F3 '  s  sum(x',2) 
location  for  byte  combining  in  a4 
use  a5  for  base  of  x^2  lookup  table 
use  d3  for  index  (points  in  profile) 
reduce  by  1  for  proper  indexing 


;At  this  point,  both  the  channel  number  and  port  number  are 


; Set  the  multiplexer  and  as  soon  as  the  A->D's  not  busy,  set  the 
; channel  number  for  the  next  point  collected.  This  also  provides 
;time  for  the  MUX  to  settle  down  after  switching  ports. 


bclr  #$03,$FFFFFC15 
bclr  #$04, $FFFFFC15 
bclr  #$05,$FFFFFC15 


; Set  port  0  on  multiplexer  board 


f3wpa    move 
f  3wpb 


$F8CO0A,d0 

bp 1  f  3wpa : 
move.b  $F8CO0A,dO 

bmi  f3wpb: 


move.b  #O,$F8C0OO 


Read  SFR  to  check  busy  flag 

; Stay  here  until  busy  flag  goes  low 
SFR  -  busy  goes  high  when  conversions  done 

; Stay  here  until  conversion  complete 


,-Set  channel  0  for  next  sample 


; Channel  0  and  port  0  are  now  set  for  the  next  sample 


f31oopa  move.b  $F8C00A,d0 


;Read  SFR  to  check  busy  flag 


bpl  f 31oopa: 
f31oopb  move.b  $F8CO0A,d0       ;SFR  -  bus? 
bmi  f31oopb: 
move.b  $F8C002, (a4)+ 
move.b  $F8C0OO, (a4) + 
move.b  #0,$F8COOO 
move . w  - ( a4 } , dl 
add.l  dl, (a0)+ 
move.l  {$0000,a5,dl.l 
add.l  62,  (al)+ 
dbf  d3,f31oopa: 
done 


;Stay  here  until  busy  flag  goes  low 
goes  high  when  conversions  done 

; Stay  here  until  conversion  complete 
;Get  MSB  from  A->D 
;Get  LSB  from  A->D 
; Set  channel  0  for  next  sample 
; Bring  entire  word  into  dl 
,-Adds  to  current  sum  for  profile  point 
r4),d2    ,-  Find  square  of  dl  from  lookup  table 
,-Adds  to  current  sum(xA2)  for  profile  point 
;Jump  back  to  main  loop  until  profile's 


,-Slow  channel  sampling  will  be  handled  from  within  this  exception 
;only,  resulting  in  100  Hz  sampling  from  the  SPI ' s . 


move.l  _spsum,a0 
move.l  #$00, d3 


,-use  aO  for  base  address  of  SPI  sums 
,-initialize  loop,  using  d3  for  index 


This  gets  a  bit  strange.   What  it  amounts  to  is,  if  the  TT  is  changing 
channels,  add  a  5  sample  delay  to  allow  things  to  settle. 

Check  if  sampling  from  channel  1,  jump  to  channel  3  portion  if  not 


dlywpa 
dlywpb 


move.l  #$05,  dl; 
move.b  $F8CO0A,d0 
bpl  dlywpa: 
move.b  $F8CO0A,d0 
bmi  dlywpb: 
move.b  #1, $F8CO00 
subq.l  #$01, dl 
bne  dlywpa 

btst.l  #3,d3 
bne  del4c3 : 


Set  index  for  delay  purposes 

Read  SFR  to  check  busy  flag 

Stay  here  until  busy  flag  goes  low 

SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 

Set  channel  1  for  next  sample 

Decrement  index 

Loop  through  5  times 

;Test  if  channel  3  should  be  used 
.-Jump  to  delay  for  channel  3  if  bit  set 
,-Otherwise,  flow  continues  with  channel  1 
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btst.l  #0,d3  ;see  if  bit  3  should  be  set 

beq  bit3zero:  ;branch  to  clear  portion  if  not 

bset.b  #3,$FFFFFC15     ; set  bit  3  -  mux  select  (lsb) 
bra  test4: 

bit3zero  bclr.b  #3,$FFFFFC15      ;clear  bit  3  -  mux  select  (lsb) 

test4     btst.l  #l,d3  ; see  if  bit  4  should  be  set 

beq  bit4zero:  jbranch  to  clear  portion  if  not 

bset.b  #4,$FFFFFC15      ;set  bit  4  -  mux  select 
bra  test5 : 

bit4zero  bclr.b  #4,$FFFFFC15     ;clear  bit  4  -  mux  select 

tests    btst.l  #2,d3  ; see  if  bit  5  should  be  set 

beq  bit5zero:  ;branch  to  clear  portion  if  not 

bset.b  #5,$FFFFFC1S      ;set  bit  5  -  mux  select  (msb) 
bra  spwpa:  ; branch  to  watch  busy  flag 

bit5zero  bclr.b  #5,$FFFFFC15     ;clear  bit  5  -  mux  select 

;Port  is  now  set  -  wait  until  A->D's  not  busy  &  set  channel 
;number.   Again,  this  gives  time  for  the  MUX  to  settle  down. 
;But,  then  again,  it  looks  like  the  mux  settles  just  fine, 
;but  the  Tattletale  is  not  so  well  behaved!  (ie.  P.O.S.) 

spwpa    move.b  $F8C00A,dO       ,-Read  SFR  to  check  busy  flag 

bpl  spwpa:  ,-Stay  here  until  busy  flag  goes  low 

spwpb    move.b  $F8C00A,dO       ;SFR  -  busy  goes  high  when  conversion's  done 

bmi  spwpb:  ;Stay  here  until  conversion  complete 

move.b  #l,$F8C00O        ;Set  channel  1  for  next  sample 
Appropriate  channel  and  port  are  now  set  for  the  next  sample 


sploopa   move.b  $F8c00A,dO        ;Read  SFR  to  check  busy  flag 

bpl  sploopa:  ;Stay  here  until  busy  flag  goes  low 

sploopb  move.b  $F8C00A,d0       ,-SFR  -  busy  goes  high  when  conversion's  done 


Stay  here  until  conversion  complete 
Get  MSB  from  A->D 
Get  LSB  from  A->D 


bmi  sploopb: 

move.b  $F8C002, (a4) 

move.b  $F8C00O, (a4) 
;move.b  #l,$F8CO00        ,-Set  channel  1  for  next  sample 

move.w  -(a4),dl         ;Bring  entire  word  into  dl 

add.l  dl,(aO)+  ;Adds  to  current  sum  for  SPI  point 

addq.l  #$01,  d3  ,- increment  loop  index 

bra  spmain  .-Unconditional  branch 


CHANNEL  3  PORTION 


del4c3 

move.l  #$05, dl; 

;del4c3 

addq.l  #$01, d3; 

move.l  #$05, dl; 

dl2wpa 

move.b  $F8COOA,dO 

bpl  dl2wpa: 

dl2wpb 

move.b  $F8COOA,d0 

bmi  dl2wpb: 

jnkch3 

move.b  #3, $F8COO0 

; Set  index  for  delay  purposes 
; TEMPORARY! ! ! 

;Set  index  for  delay  purposes 
;Read  SFR  to  check  busy  flag 
;Stay  here  until  busy  flag  goes  low 
;SFR  -  busy  goes  high  when  conversions  done 
;Stay  here  until  conversion  complete 
;Set  channel  3  for  next  sample 


jnklp     subq.l  #$01, dl 
bne  dl2wpa 


; Decrement  index 
,-Loop  through  5  times 


spmain2  btst.l  #0,d3 

beq  bit3zer2 : 

bset.b  #3, 

bra  test4b: 
bit3zer2  bclr.b  #3,$FFFFFC15 
test4b   btst.l  #l,d3 

beq  bit4zer2 : 

bset.b  #4, 

bra  test5b: 
bit4zer2  bclr.b  #4,$FFFFFC15 
test5b   btst.l  #2,d3 

beq  bitSzer2 : 

bset.b  #5, 

bra  spwpa2 : 
bit5zer2  bclr.b  #5,$FFFFFC15 


;see  if  bit  3  should  be  set 
,-branch  to  clear  portion  if  not 
$FFFFFC15      ;set  bit  3  -  mux  select  (lsb) 

,- clear  bit  3  -  mux  select  (lsb) 
;see  if  bit  4  should  be  set 
,-branch  to  clear  portion  if  not 
$FFFFFC15      ;set  bit  4  -  mux  select 

,- clear  bit  4  -  mux  select 
,-see  if  bit  5  should  be  set 
,-branch  to  clear  portion  if  not 
5FFFFFC1B      ;set  bit  5  -  mux  select  (msb) 
,-branch  to  watch  busy  flag 
,- clear  bit  5  -  mux  select 
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,-Port  is  now  set  -  wait  until  A->D's  not  busy  &  set  channel 
,-number.   Again,  this  gives  time  for  the  MUX  to  settle  down. 
,-But,  then  again,  it  looks  like  the  mux  settles  just  fine, 
;but  the  Tattletale  is  not  so  well  behaved!  (ie.  P.O.S.) 


spwpa2  move.b  $F8CO0A,d0 
bpl  spwpa2 : 

spwpb2  move.b  $F8CO0A,dO 
bmi  spwpb2 : 


Read  SFR  to  check  busy  flag 

Stay  here  until  busy  flag  goes  low 

SFR  -  busy  goes  high  when  conversion's  done 

Stay  here  until  conversion  complete 


move.b  #3,$F8C00O       ;Set  channel  3  for  next  sample 
;Appropriate  channel  and  port  are  now  set  for  the  next  sample 


sploopa2  move.b  $F8C00A,dO 
bpl  sploopa2 : 

sploopb2  move.b  $F8C0OA,d0 
bmi  sploopb2 


Read  SFR  to  check  busy  flag 

Stay  here  until  busy  flag  goes  low 

SFR  -  busy  goes  high  when  conversion's  done 

Stay  here  until  conversion  complete 

move.b  $F8C002, <a4) +     ;Get  MSB  from  A->D 

move.b  $F8CO0O,  (a4)+     ,-Get  LSB  from  A->D 
;move.b  #l,$F8COO0        ;Set  channel  1  for  next  sample 

move.w  -{a4),dl  ;Bring  entire  word  into  dl 

add.l  dl,  (a0)+  ,-Adds  to  current  sum  for  SPI  point 

addq. 1  #$01, d3  ; increment  loop  index 

move.l  _nosp,d2  ;get  value  of  nosp  for  comparison 

cmp.l  d2,d3  ,-check  if  all  SPI '  s  have  been  sampled 

bne  spmain2:  .-branch  back  if  not 


; Since  F3  is  the  final  transducer  from  which  profiles  are 

;the  number  of  unaveraged  profiles  still  remaining  will  be 
.-determined  from  within  this  exception. 


subq.l  #l,_nupc         .-Decrement  number  of  profiles  collected 

eori.w  #$0200 , $FFFE20   ; Acknowledges  and  turns  off  interrupt 

movem.l  (a7) +,  d0-d7/a0-a6   ,-Pull  registers  from  stack 
rte 


ttendasm 
} 


*  MTA  data  retrieval  routine 

/*  grabdata  gets  the  data  from  the  MTA.  Routine  completes  when  number 
of  bytes  received  =  either  72*mta_numscans  or  40*mta_numscans. 
If  the  routine  times  out,  stringerror=2 ,  no  error,  stringerror=0.  */ 

void  grabdatafchar  mta_id) 
{ 

int  stringerror=0; 
long  limit=0; 
long  timeout=0; 
long  stringpoint=0; 
unsigned  char  "mtabas; 

if  (mta_id=='l' ) 
{ 
mtabas= (unsigned  char  M  mtlbas; 
limit=72*mta_numscans; 
} 
if  (mta_id== ' 2 ' ) 
{ 

mtabas= (unsigned  char  *)  mt2bas; 
limit-40*mta_numscans; 
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} 
if  (nita_id« '  3 ' ) 

{ 

mtabas= (unsigned  char  *)  mt3bas; 
limit=40*mta_numscans; 

} 

while (timeout<2 000 000) 
{ 

if  (TSerByteAvail (inchan) ) 
{ 

timeout=0; 

ch-TSerGetByte ( inchan) ; 
* (mtafaas+stringpoint) =ch; 
stringpoint+=l; 
if  (stringpoint== (limit) ) 
{ 

stringerror-0; 
goto  grabdata_bail; 
) 
} 
timeout+-l  ,- 
1 
stringerror=2 ; 

grabdata_bai 1 : 

printf  ( " " )  ,- 

} 


*  HTA  string  retrieval  routine 

/*  grabstring  gets  a  string  from  the  MTA.   If  the  string  overruns  (>20  bytes) 
stringerror=l,  if  the  routine  times  out,  stringerror=2 ,  no  error, 
stringerror=0 .  */ 

void  grabstring(char  mta_id) 
{ 

long  timeout=0; 
stringpoint-0; 
stringerror=0; 

while (timeout<600000) 
{ 

if  (TSerByteAvail (inchan) ) 
{ 

timeout=0; 

ch-TSerGetByte{  inchan)  ,- 
instring[stringpoint] =ch; 
if  (ch!^0) 

{ 

printf ("%c" ,ch) ; 

stringpoint+-l ; 

} 
else 

{ 

stringpoint-0; 

) 
if  (stringpoint>19) 

{ 

stringerror=l; 

goto  grabstring_bail; 

} 
if  (ch-=13) 

{ 

stringpoint-=l; 

instring[stringpoint] =0; 

printf <"\n") ; 

goto  grabstring_bail; 


119 


/*  Mux  select  LSB  */ 

/*  Mux  select  */ 

/*  Mux  select  MSB  */ 

/*  TPU  buffer  disable  <0=enable) 


timeout+-l ; 
} 

grabstring_bail : 

printf ("" ) ; 
} 


*  Initialization  routine 

void  initialize (void) 
( 

char  m; 

long  baud_rate; 
long  sys_freq; 
long  tma; 
time_t  timer; 
struct  tm  time_struct, *vtime_struct; 

printf (UNDOFF) ; 

printf (CLEAR) ; 

printf (UNDON) ; 

printf { " Initialization. " ) ; 

printf (UNDOFF) ; 

/*  Digital  I/O  configuration  */ 
/*  Multiplexor  */ 
PConfOutp(D, 3) 
PConf0utp(D,4) 
PConfOutp(D,5) 
PConfOutp(E,2) 
PSet(E,2) ; 
e2  =  l; 

/*  Charging  */ 

PConfOutpfE, 5) ;  /*  Battery  charge  (l=off) 

PClear(E,5} ; 
e5=0; 

/*  Instrument  power  */ 
e6=0;     /*  Tilt,  compass,  temp  off  */ 
PConf0utp(E,6) ; 
PClear(E,6) ; 

e7=0;     /*  OBS/pressure  off  */ 
PConfOutpfE, 7) ; 
PClear(E,7) ; 

eO-0;     /*  ADV  #1  off  "/ 
PConfOutp(E,0) ; 
PClear(E,0) ; 

fl=0;     /*  ADV  #2  off  */ 
PConfOutp(F,  1)  ,- 
PClear(F,l) ; 

el=0;     /•  Video  off  */ 
PConfOutpfE, 1) ; 
PClear(E,l); 

f6-0;     /*  ACP  (3  frequency)  off  */ 
PConf0utp(F,6) ; 
PClear(F,6)  ,- 

f7=0;     /*  MTA'S  off  */ 
PConfOUtpfF, 7) ; 
PClear(F,7) ; 
df7=0;    /*  MTA's  enabled  */ 

/*  ABS  Gain  -  default  to  low  gain  slope  */ 

tpuO=0; 

TPUSetPin(0,OJ ; 
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/*  Disk  drive  power 
DriveOff 0 ; 


/*  Software  flags  */ 

/*  Sampling  mode 
samp_mode=0  ,- 
storm-0; 


/*  Disk  data  file  pointer  */ 

data_f ile_pointer=l; 

max_f iles=249;        /*  Maximum  number  of  disk  data  files  */ 

num_to_fill=l;       /*  Number  of  data  files  to  fill  */ 

/*  Default  sampling  scheme  definition  */ 

nosp=14;  /*  Number  of  SPI ' s  to  sample  */ 

perm_npp=120;  /*  Number  of  points  per  profile  (3  freq. )  */ 

perm_npa=22;  /*  Number  of  profiles  to  include  in  average  */ 

perm_prof_period=25;  /*  Storage  period  in  hundredths  of  seconds  */ 

shore_nfp=0;  /*  Number  of  final  profiles  to  collect  =  0=max  */ 

total_nfp-l;  /*  Either  #  of  shore  profs  or  1  for  inst_test  */ 

dif_nfp=l;  /*  Shore  prof iles -memory  profiles 
mein_nf  p=0 ; 

offset=100;  /*  Hours  from  midnight  if  timed  mode  */ 

perm_sdt=0;  /*  Shore  data  time  in  hours  */ 

hours_between=400;  /*  .01's  hours  between  collections  if  timed  mode  */ 
num_collections=max_f iles;  /*  Number  of  collections  if  timed  mode  */ 

/*  MTA  parameters  */ 

mta_blank-325; 
mta__delay=200; 

mta_gain-123 ; 

mta_pulse=12; 
mta_thresh [  1 ] =300 ; 
mta_thresh[2] =700; 
mta^thresh [ 3 ] =3  000 ; 

/*  Default  time  &  date  */ 

time„struct. tm_sec=0; 
time_struct. tm_min=0; 
time_struct. tm_Jiour=0; 
time_struct. tm_mday=l; 
time_struct . tm_mon=3 ; 
time_struct. tm_year=97; 
time_struct. tm_isdst=0; 
SetTimeTM(&time_struct,NULL) ; 

time_struct. tm_sec=0; 
time_struct . tm_min=0 ; 
time_struct. tm_hour=4; 
time_struct . tm_mday=l ; 
time_struct . tm_mon=3 ; 
time_struct . tm_year=97; 
timer=mktime (&time_struct) ; 
start_time- timer ; 

/*  System  frequency  &  baud  rate  */ 

printf ( "\nCurrent  baud  rate  ■  %ld" ,baud_rate=SerGetBaud(0, 0) ) ; 
printf ( " \nCurrent  system  frequency  =  %ld  Hz" , sys_freq=SimGetFSys ( ) ) ; 
if  (sys_freq 1=16000000) 

{ 

SimSetFSys  (16000000)  ,- 

SerSetBaud(baud_rate,0)  ; 

printf {" \nSystem  frequency  changed  to  %ld  Hz" , sys_f req=SimGetFSys ( ) ) 

} 
if  (baud_rate! =384  61) 

{ 

printf  ("  \nChanging  baud  rate  to  38400.   "),- 

printf (" \nChange  terminal  rate  and  press  RETURN"); 

SerSetBaud(38400, 0) ; 
baud_change : 

m- (char)menu_input (0,0,0) ; 
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if  (m!=10)  goto  baud_change; 
} 

/*  PS  ram  allocation  */ 
PSMemFreeAll () ; 
sqbase=PSMemAllocAll  (&tma)  ; 
if  (tma==0) 

{ 

printf ( " \nMemory  allocation  problem!"); 

) 
tmal28=tma-128;    /*  tma-128  *f 
rtbase-sqbase+4096; 
flsum=rtbase+4096; 

/*  x"2   lookup  table  */ 
#asm 

movem.l  dO-d7/aO-a6, - {a7)   ; Push  registers  on  stack 

;  x-2   table 

move.l  _sqbase,aO  ;use  aO  for  base  of  x"2  lookup  table 

move.l  #$0FFF,d0  ; start  with  highest  value 

sloop  move.l  dO,dl  ,- transfer  value  to  dl  for  squaring 

mulu  dO,dl  ; square  it 

move.l  dl, ($0000, aO,dO. 1*4)    ;store  it  using  root  as  index 

dbf  dO, sloop:  ;continue  for  entire  table 

;  A->D  setup 

move.b  #O,$F8CO00  ; Set  channel  zero  on  A->D  command  register 

move.b  #5,$F8CO0A  ;Turn  on  latch,  allow  TH  sampling  (A->D  SFR) 

movem. 1  (a7 ) +, d0-d7/a0-a6     ; Pull  registers  from  stack 
ttendasm 

/*  TPU  setup  */ 
tpu_setup ( ) ; 
} 

*  Instrument  test  routine 

void  inst_test (void) 
{ 

char  ack; 
char  m; 
int  i; 

int  transducer ; 
unsigned  int  xvalue; 
unsigned  int  yvalue; 
unsigned  long  midway; 
unsigned  short  *base; 

test_menu: 

printf (CLEAR) ; 
printf (UNDON) ; 

printf ( " Instrument  tests\n\n" ) ; 
printf (UNDOFF) ; 
printfC'l  -  View  Fl  transducer  (1.00  MHz)  k   SPI's\n"); 
printf ("2  -   View  F2  transducer  (2.20  MHz)  k   SPI's\n"); 
printf ("3  -  View  F3  transducer  (5.00  MHz)  &  SPI-s\n"); 
printf ("S  -  View  SPI ' s  only\n"); 
printf ( "ENTER  -  Return  to  main  menu\n\n"); 
printf ( "  %c" 

,13); 
printf ( "Selection  =>"); 
test_selection : 

m- (char)menu_input (0,0,0) ; 
if  ( (m==(-l) ) | | (m==10) )  goto  test_bail; 
if  (!  (  (m=='l')  j  |  (m==-2'l  |  |  (m-'3'  )  |  |  <m=='S'  ))) 
goto  test_selection; 
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printf ("%c   ENTER  to  accept,  any  other  key  to  reject .%c" ,m, 13) ; 
ack= ( char ) menu_input (0,0,0) j 
if  (ack!=10) 

{ 

printf (" 

,13); 

printf ( "Selection  ->"); 

goto  test_selection; 

} 


itch(m) 

case 

■1'  : 

transducer^ 

break; 

case 

•2-   : 

transducer^ 

break; 

case 

'3'  : 

transducer^ 

break; 

case 

'S'  : 

transducer=0 

break; 

} 

if  ( transducer ! =0 ) 
{ 

printf ( "\n\nEnabling  TPU  buffers."); 
PClear(E,2) ,- 
} 

printf { " \n\nlnitializing  memory. . . " ) ; 
nfpc=l; 
npp=12  0; 
prof_period~2  5; 
npa=22; 
total_jifp=l; 
mem_setup ( ) ; 
mat_headers ( ) ; 
clear_mem( ) ; 
printf (CLEAR) ; 
main_sample( ) ; 

/*  Set  base  address  of  data  storage  */ 
if  (transducer==l ) 

base=f lrms; 
else  if  (transducer==2) 

base=f 2rms; 
else  if  (transducer- =3) 

base=f3rms; 

/*  Screen  display  */ 

while(kbhit()==0> 
f 
if  (transducer ! -0) 

{ 

printf (CLEAR) ; 

printf ("\xlB[24;lf"}; 

printf  ("F%d  transducer   ".transducer)  ,- 

} 
printf ("SPI's  =>") ; 
for  (i=0; i<nosp;i++) 

{ 

printf ( "%04u  " , * (spmen+i) ) ; 

} 

if  (transducer! -0) 
{ 

for  (i=0;i<60;i++) 
{ 
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midway= ( (unsigned  long) * (base+i*2 } ) *2 

yvalue= (unsigned  short) (23-midway}; 

xvalue=i; 

print f ( "\xlB[%d;%df*",yvalue,xvalue) ; 

} 

/*  DelayMilliSecs  (250)  ,-  */ 


printf ( "\n" ) ; 
main_sample ( } j 
) 


/*  Disable  TPU  buffers 
PSet(E,2) ; 
goto  test_menu; 

test_bail : 
printf (CLEAR) ; 


Main  sampling  routine 


void  main_sample (void) 

( 
#asm 


movem.l  dO-d7/aO-a6, - (a7) 


; Push  registers  to  stack 


Initialization 


; Initialize  number  of  profiles  to  average  and  collect 

move.l  _npa,dO       ;get  value  of  npa 

addq.l  #$01, dO        ,-add  1  to  initialize  _nupc 

move.l  d0,_nupc       ;Init  nupc  so  nothing  governs  initialization 

move.l  _total_nfp,_nfpc  ;_total_nfp  won't  change,  nfpc  decreases 


move. 1  #0,_storf 1 


; initialize  to  zero  so  storage  occurs 


;Make  certain  nothing  exception  runs  first 


lea  noexc,a0 
move.l  aO, $0924 
move.l  aO, $0928 
move.l  aO, $092C 


Start  of  nothing  exception 
Set  TPU  9 ' s  exception  address 
Set  TPU  10' s  exception  address 
Set  TPU  11 's  exception  address 


; Initialize  bases  for  profile  storage 


move.l  _f lrms,_tf lbas 
move.l  _f 2rms,_tf 2bas 
move.l  _f 3rms,_tf 3bas 
move.l  _spmen, _tspmen 


,-copy  starting  base  to  variable  base  for  Fl 

; variable  base  for  F2 

,-variable  base  for  F3 

,- variable  base  for  SPI '  s 


; Enable  TPU  interrupts 
move.w  #$2400, sr 
move.w  #$OEO0, $FFFE0A 
move.w  #$2000, sr 


Mask  all  interrupts 
Enable  TPU  interrupts 
Unmask  all  interrupts 


,■  Main  sampling/averaging  loop 

sync     move.l  _npa,d0        ;Use  dO  to  compare  with  nupc 

cmp.l  _nupc,d0        ;Wait  here  until  nothing  initializes  nupc 
bne  sync:  ,-This  guarantees  Fl  will  be  sampled  first 

numtoav  move.l  _nupc,d0       ;Use  dO  to  check  number  of  profiles  collected 

bne  numtoav:         ;loop  until  npa  profiles  collected 


lea  noexc.aQ 
move.l  a0,$0924 
move.l  aO, $0928 


;get  address  of  nothing  exception 
; change  TPU  9 ' s  entry 
; change  TPU  10' s  entry 
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move.l  aO,S092C       .-change  TPU  ll's  entry 

; Averaging  routine  begins  here 

; Output  header  bytes  to  parallel  i/o  port 
move.l  _pioport,a5 
move.b  #$FF, (a5) 
move.b  __prof_index, (a5) 
add.b  #S01,_prof_index 
cmp.b  #$FE,_prof_index 
bne  nidxres : 
move.b  #$01,_j?rof_index 

; Initialize  RMS  routine  for  Fl 
nidxres  move.l  _flsqr,aO     ;Set  base  of  Fl ' s  sum(x~2) 

move.l  _tflbas,al    ; Set  base  of  Fl ' s  RMS  storage 
move.l  _rtbase,a2    ;Set  base  of  rms  table 

move.l  _wordprof,dO   ;get  number  of  bytes  to  add  to  previous  base 
add.l  dO,_tflbas      ; compute  new  base  for  comparison 

move.l  _npp,d6       ;initialize  number  of  points  per  profile 

;RMS  routine  -  Fl 
avloopfl  move.l  (aO),dO      ;get  byte  from  sum(x~2) 

move.l  #0,(a0)+       ,-initialize  sum(x"2)  storage 

move.l  #$0800, dl     ,-initialize  TEST 

move.l  #$0400, d4     ,-initialize  HALF 

move.l  #0,d3  initialize  BASE 

move.l  #0,d2         ,-initialize  ANSWER 
sredofl  lea  ($0000, a2,dl .w*4) ,a5    ; location  of  TABLE[TEST] 

cmp.l  (a5),d0 

bge  felsefl: 

move.w  dl, d2 

bra  dunfl: 
felsefl   move.w  dl,d3  ;BASE=TEST  ->  else  portion 

dunfl     move.w  d4,d5  ;WORK=HALF 

or.w  d3,d5 

move.w  d5,dl 

Isr.w  #l,d4 

bne  sredofl: 

bcs  sredofl : 


branch  to  else  portion  if  Y>=TABLE[TEST] 

ANSWER=TEST 

skip  else  portion 


WORK=HALF  or  BASE 
TEST=HALF  or  BASE 
shift  HALF  right  1  bit 
branch  back  10  times 
branch  back  1  final  time 


; Output  to  parallel  i/o 

move.w  d2,_bytecomb  ,-send  answer  to  buffer 

lea  _bytecomb,a4 

move.l  „pioport,a5 

cmpi.b  #$FF, (a4)      ;check  for  FF  in  msb 

bne  no_f  f la : 

move.b  #$FF, (a5)  ; send  double  FF ' s  for  each  single  FF 
no_ffla  move.b  (a4)+, (a5)     ;send  msb  to  parallel  i/o  port 

cmpi.b  #$FF, (a4)      ,-check  for  FF  in  lsb 

bne  no_fflb: 

move.b  #$FF, (a5)  ;send  double  FF's  for  each  single  FF 
no_fflb  move.b  (a4) , (a5)      ;send  lsb  to  parallel  i/o  port 

,-Check  _storflag  -  skip  to  nostorl  if  flag  nonzero 

cmpi . 1  #0,_storf lag 

bne  nostorl: 

move.w  d2 , (al)+      ;store  rms  value 

nostorl  subq  #$01,  d6         ,-decrement  profile  point  count 

bne  avloopfl:         ;jump  back  if  not  zero 

,- Initialize  RMS  routine  for  F2 

move.l  _f2sqr,a0      ;  Set  base  of  F2  '  s  sum(xA2) 

move.l  _tf2bas,al     ;Set  base  of  F2 ' s  RMS  storage 

move.l  __wordprof ,  dO   ;get  number  of  bytes  to  add  to  previous  base 
add.l  d0,_tf2bas     ;compute  new  base  for  comparison 
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move . 1  _npp , d6 


; initialize  number  of  points  per  profile 


,-RMS  routine 


F2 


avloopf2  move.l  (a0),d0 


,-get  byte  from  sum(xA2) 
move.l  #0,(a0)+      ;initialize  sum(x"2) 
move.l  #$0800, dl     jinitialize  TEST 
move.l  #$0400, d4      ;initialize  HALF 


sredof2   lea  ($0000, 


felsef2   move.w  dl , d3 
dunf  2    move.w  d4,d5 


move.l  #0,d3 

move.l  #0,d2 

,dl.w*4) ,a5 

cmp.l  (a5),d0 

bge  felsef 2 : 

move.w  dl,d2 

bra  dunf 2 : 

;BASE=TEST 
; WORK = HALF 

or.w  d3 , d5 

move.w  d5,dl 

lsr.w  #l,d4 

bne  sredof2 : 

bcs  sredof2 : 


.•initialize  BASE 
; initialize  ANSWER 
location  of  TABLE [TEST] 


branch  to  else  portion  if  Y>=TABLE[TEST] 
ANSWER=TEST 
skip  else  portion 
else  portion 

WORK=HALF  or  BASE 
TEST=HALF  or  BASE 
shift  HALF  right  1  bit 
branch  back  10  times 
branch  back  1  final  time 


,-Output  to  parallel  i/o 

move.w  d2,_bytecomb  ; send  answer  to  buffer 

lea  _bytecomb,a4 

move.l  _pioport,a5 

cmpi.b  #$FF, (a4)  ;check  for  FF  in  msb 

bne  no_ff2a: 

move.b  #$FF,  (a5)  ,-send  double  FF '  s  for  each  single  FF 

no_ff2a  move.b  (a4)-t-,(a5)     ,-send  msb  to  parallel  i/o  port 

cmpi.b  #$FF, (a4)  ; check  for  FF  in  lsb 

bne  no_ff2b: 

move.b  #$FF,(a5)  ;send  double  FF*s  for  each  single  FF 

no_ff2b  move.b  (a4),(a5)      ,-send  lsb  to  parallel  i/o  port 

,-Check  „storflag  -  skip  to  nostor2  if  flag  nonzero 

cmpi . 1  #0,_storf lag 

bne  nostor2: 

move.w  d2, (al)+      ; store  rms  value 

nostor2   subq  #$01, d6         .-decrement  profile  point  count 

bne  avloopf2:         ;jump  back  if  not  zero 

; Initialize  RMS  routine  for  F3 

move.l  _f3sqr,a0      ;Set  base  of  F3  '  s  sum(x'v2) 

move.l  _tf3bas,al    ,-Set  base  of  F3  '  s  RMS  storage 


move . 1  _wordprof , dO 
add.l  d0,_tf3bas 


move . 1  _npp , d6 


;get  number  of  bytes  to  add  to  previous  base 
; compute  new  base  for  comparison 


,- initialize  number  of  points  per  profile 


;RMS  routine 


F3 


avloopf3   move.l    (a0),d0 


move, 
move, 
move. 


;get  byte  f 
#0, (a0)+ 
#$0800, dl 
#$0400, d4 
#0,d3 
#0,d2 


sredof3   lea  ($0000, a2,dl ,w*4) , a5    ; 
cmp.l  (a5) ,d0 


sum(x"2) 
,-initialize  sum(xA2)  storage 
; initialize  TEST 
,-initialize  HALF 
,-initialize  BASE 
,-initialize  ANSWER 
cation  of  TABLE [TEST] 


felsef 3  move.w  dl,d3 
dunf 3    move.w  d4,dB 


bge  felsef 3 : 

move.w  dl, d2 

bra  dunf 3 : 

;BASE=TEST 
;WORK=HALF 

or .  w  d3  ,  d5 

move.w  d5, dl 

lsr.w  #l,d4 

bne  sredof3: 

bcs  sredof 3 : 


;branch  to  else  portion  if  Y>-TABLE [TEST] 
;ANSWER=TEST 
,-skip  else  portion 
else  portion 

;WORK=HALF  or  BASE 
;TEST=HALF  or  BASE 
,- shift  HALF  right  1  bit 
;branch  back  10  times 
;branch  back  1  final  time 
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.-Output  to  parallel  i/o 

move.w  d2,_bytecomb   ; send  answer  to  buffer 

lea  _bytecorab,a4 

move.l  _pioport,a5 

cmpi.b  #$FF, (a4)      ; check  for  FF  in  msb 

bne  no_ff3a: 

move.b  #$FF,  (a5)  ,-send  double  FF '  s  for  each  single  FF 
no_ff3a  move.b  (a4)+,(a5)     ; send  msb  to  parallel  i/o  port 

cmpi.b  #$FF, (a4)      ;check  for  FF  in  lsb 

bne  no_ff3b: 

move.b  #$FF, (a5)  ;send  double  FF ' s  for  each  single  FF 
no_ff3b  move.b  (a4)  ,  (a5)      ,-send  lsb  to  parallel  i/o  port 

,-Check  _storflag  -  skip  to  nostor3  if  flag  nonzero 

cmpi . 1  #0,_s tor flag 

bne  nostor3: 

move.w  62,  <al)  +      ,-store  rms  value 

nostor3   subq  #$01, d6         ;decrement  profile  point  count 

bne  avloopf3:        ;jump  back  if  not  zero 

;bra  skipsp: 

; Compute  mean  of  each  SPI  channel 

move.l  _spsum,aO     ,-set  base  of  SPI  sum(x)  storage 
move.l  _tspmen,al     ;set  base  of  SPI  mean  storage 

move.l  _spichans,dO   ;get  number  of  bytes  to  add  to  previous  base 
add.l  dO,_tspmen     .-compute  new  base  for  comparison 

move.l  _nosp,d6       ; initialize  spi  channel  count 

avloopsp  move.l  (a0),d0       ,-get  byte  from  sum(x) 

move.l  #$00,{a0)+    ^-initialize  sum(x)  storage 
divu.l  _npa,dO       ; divide  by  npa 

;0utput  to  parallel  i/o 

move.w  dO,_bytecomb   ; Send  answer  to  buffer 

lea  _bytecomb,a4 
move.l  _pioport,a5 

cmpi.b  #$FF, (a4)      ; check  for  FF  in  msb 

bne  no„ff4a: 

move.b  #$FF, (a5)  ; send  double  FF ' s  for  each  single  FF 
no_ff4a  move.b  (a4)+,(a5)     ;send  msb  to  parallel  i/o  port 

cmpi.b  #$FF, (a4)      ;check  for  FF  in  lsb 

bne  no_ff4b: 

move.b  #$FF, (a5)  ; send  double  FF's  for  each  single  FF 
no_ff4b  move.b  (a4),(a5)      ;send  lsb  to  parallel  i/o  port 

;Check  _storflag  -  skip  to  nostor4  if  flag  nonzero 

cmpi. 1  #0,_s tor flag 

bne  nostor4 : 

move.w  dO,  (al)+       ,-store  result 

nostor4   subq  #$01, d6  ;decrement  channel  count 

bne  avloopsp:         ;jump  back  if  not  zero 

^compare  nfpc  to  dif_nfp,  set  _storflag  (nonzero)  if  equal 

skipsp   move.l  _dif_nfp,d0 

cmp. 1  _nfpc,dO 

bne  stilstor: 

move .1  # 1 , _s tor flag 

;Check  if  sampling  is  complete,  return  if  not 
stilstor  subq.l  #l,_nfpc      ,-decrement  number  of  final  profiles 

bne  sync:  .-continue  until  all  final  profiles  collected 

;  End  of  main  sampling /averaging  loop 
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move.w  #$2400, sr 
move.w  #$0000,$FFFEOA 
move.w  #$2000, sr 


;Mask  all  interrupts 
,- Disable  TPU  interrupts 
.-Unmask  all  interrupts 


;  Return  to  C 

movera.l  (a7) +,d0-d7/a0-a6   ,-Pull  registers  from  stack 
#endasm 


Matlab  header  setup  routine 


id 

mat_headers (void) 

{ 

unsigned  long  *header; 

unsigned  char  "nameptr; 

header=f lhed; 

*header++=1040; 

/• 

type  */ 

*header++=npp; 

/* 

mrows  * / 

♦header ++=nfp ; 

/* 

ncols  */ 

*header++=0; 

/* 

imagf  * / 

*header++-6; 

/* 

namlen  */ 

nameptr= (unsigned 

char 

*) (header) 

*nameptr++- ' f ' 

*nameptr++=' 1 ' 

*nameptr++- ' r ' 

*nameptr++='nT 

*naraeptr++- ' s ' 

*nameptr++=0; 

header=f2hed; 

*header++-1040; 

/* 

type  *  1 

*header++=npp; 

/* 

mrows  * / 

*header++=nfp; 

/* 

ncols  */ 

*header++-0; 

/* 

imagf  */ 

*header++-6; 

/* 

namlen  */ 

nameptr= (unsigned 

char 

*) (header) 

*nameptr++= ' f ' 

*nameptr++= ' 2 ' 

*nameptr++= ' r ' 

*nameptr++= 'm' 

*nameptr++- ' s ' 

*nameptr++-0; 

header=f 3hed; 

*header++=1040; 

/* 

type  */ 

*header +  +  =npp ; 

/* 

mrows  * / 

* header ++=nfp  ; 

/* 

ncols  */ 

*header++=0; 

/* 

imagf  */ 

*header++=6; 

/* 

namlen  */ 

nameptr= (unsig 

led 

char 

*) (header) 

*nameptr++= ' f 
*nameptr++= ' 3 
*nameptr++= ' r 
*nameptr++= 'm 
*nameptr++= ' s 
*nameptr++=0; 


header=sphed; 
*header++=1040; 
* header ++-nosp; 
* header ++-nfp ; 
*header+  +  -0,- 
"header++=6; 


type  */ 
mrows  * 
ncols  * 
imagf  * 
namlen 


/*  variable  name 


variable  name 


/*  variable  name 


nameptr-  (unsigned  char  *}  (header)  ,- 


variable  name 


*nameptr++ 
*nameptr++ 
*nameptr++ 
*nameptr++ 
*nameptr++ 
*nameptr++ 
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heade r =mt lhed ; 

*header++=1040; 

*header++=36; 

•header ++=mta_numscans ; 

*header++=0; 

*header*-+=6; 

nameptr= (unsigned  char 

*nameptr++- 'm' ; 

*nameptr++= '  t '  ; 

*nameptr++= 'a'  ; 

*nameptr++= ' r '  ; 

*nameptr++= '  1 '  ; 

*nameptr++=0; 


/*  type  ' 
/*  mrows 
/*  ncols  n 
/*  imagf 
/*  namlei 
) (header) ; 


/*  variable  name  =  mtarl  */ 


header -mt2hed; 

*header++=1040;  /*  type  */ 

*header++=20;  /*  mrows  */ 

*header++=mta_numscans;  /*  ncols  */ 

*header++=0  ?  / *  imagf  * / 

*header++-6;  /*  namlen  * 
nameptr= {unsigned  char  *) (header);   /* 


variable  name  -  mtar2 


*nameptr++="m' 

*nameptr++= ' t ' 

*nameptr++= 'a' 

*nameptr++='r ' 

*nameptr+  +  - ' 2  ' 

*nameptr++=0; 

header =mt3hed; ! 

*header++=1040; 

/*  type 

*header++=20; 

/*  mrows 

*header++=mta_numscans; 

/*  ncols 

*header++=0; 

/*  imagf 

*header++-6; 

/*  namle 

nameptr- (unsigned 

char 

* )  (header)  ,- 

*nameptr++= 'm' 

*nameptr++= ' t* 

*nameptr++='a' 

"nameptr++= 'r ' 

*nameptr++= ' 3  ' 

*nameptr++=0; 

header- tmhed  ; 

*header++=1040; 

/* 

type  */ 

*header++=6; 

/* 

mrows  * / 

*header++-2; 

/* 

ncols  */ 

*header++=0; 

/* 

imagf  */ 

*header++=6; 

/* 

namlen  */ 

nameptr= (unsigned 

char 

*) (header) ; 

*nameptr++= ' t ' 

*nameptr++= ' i ' 

*nameptr++= 'm' 

*nameptr++= 'e ' 

*nameptr++= ' s ' 

*nameptr++=0; 

header =dmhed, • 

*header++=1040; 

/* 

type  */ 

•header ++=pad/ 2; 

/* 

mrows  */ 

*header++=l; 

/* 

ncols  */ 

*header++=0; 

/* 

imagf  * / 

*header++=6; 

/* 

namlen  * / 

nameptr= (unsigned 

char 

* ) (header) ; 

*nameptr++= 'd' 

/*  variable  name  =  mtar3 


/*  variable  name 


variable  name  =  dummy 
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*nameptr++= 'u- ; 
*nameptr++='m' ; 
*nameptr++= '  m '  ; 
*nameptr++=  'y'  ; 
*nameptr++=0; 
} 


Memory  setup  routine 


void  mem_setup (void) 
{ 

float  mta_per_prof ; 
unsigned  long  tmu; 

mta_per_prof=( (float) (pernuprof_period) ) / ( (float) (mta_delay) ) ,- 
nfp={tmal28-24*npp-8*nosp-3302  6) / (6*npp+2*nosp+ 
( (int) (mta_per_prof*152. 0+1.0) ) ) ; 
if  (nfpoO) 

{ 

nfp=nfpc; 
} 
if  (total_nfp==0) 
{ 

total_nfp=nfp; 
) 
mta_numscans=(int) <mta_per_prof * ( (float) (nfp) ) ) ; 
shore_mta_numscans= ( int) (mta_per_prof * ( ( float) (total_nfp) ) ) ; 
if  (shore_mta_numscans  >  MTA„MAXSCANS ) 
( 

shore_mta_nurascans  -  MTA_MAXSCANS; 


tmu-nfp* (6*npp+2*nosp) +mta_numscans*152+26 
f 1 sqr =  f 1 sum+npp ; 
f  2  sum- f 1 s  gr + npp ; 
f 2 sqr- f 2 sum+npp ; 
f 3sum=f2sqr+npp; 
f  3  sqr  =  f  3  sum+npp ; 
spsum=f3sqr+npp; 
spsgr-spsum+nosp; 
f lhed=spsqr+nosp ; 
flrms= (unsigned  short 
f2hed= (unsigned  long 
f2rms= (unsigned  short  *)( (unsigned  char 
f3hed- (unsigned  long  *) (f 2rms+npp*nfp) 
f 3rms= (unsigned  short 
sphed= (unsigned  long  *} (f3rms+npp*nfp) 
spmen=  (unsigned  short  M(  (unsigned  char 
mtlhed- (unsigned  long 
mtlbas- (unsigned  short 
mt2hed= (unsigned  long 
mt2bas- (unsigned  short 
mt3hed= (unsigned  long 


*) ( (unsigned  char 
*) (f lrms+npp*nfp) 


*)( (unsigned  char  *) (f3hed) +26) 


) (sphed)+26) 

26) 


) (spmen+nosp*nfp) ; 
*) ( (unsigned  char  * ) (mtlhed) 
) (mtlbas+mta_numscans*36) ; 
*) ((unsigned  char  *) (mt2hed) 
) (mt2bas+mta_numscans*20) ; 


26) 
26) 

6); 

6); 


mt3bas= (unsigned  short  *) ({unsigned  char  *) (mt3hed) 

tmhed= (unsigned  long  *) (mt3bas+mta_numscans*20) ; 

tmbas- (unsigned  short  *> ((unsigned  char  *)(tmhed)+2 

dmhed= (unsigned  long  *) (tmbas+12) ; 

dmbas= (unsigned  short  *) ( (unsigned  char  *) (dmhed) 

pad= (tmu/ 128)* 12 8 +12 8- tmu; 

of fstart- (unsigned  char  *) (fined); 

of fend- (unsigned  char  *) (dmbas) +pad-l ; 

/*  Set  values  for  number  of  bytes  in  word  profiles  and  all  SPI  channels  V 
wordprof -npp*2 ; 
spichans-nosp*2 ; 
) 


*  Menu  input  routine 


int  menu_input ( int  a,int  d,int  n) 
{ 

bool  x; 
char  b; 
int  m; 
int  v; 

1  ong  t  ime__ou  t = 0 ; 
void  display_n(int, int, int) ; 
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kbflushO  ; 

if  (a!=0)  display_n(a,d,n) ; 


clean  buffer  */ 

display  default  if  numeric 


check_buf : 

while (kbhit()= 


/*  wait  for  key 


time_out+=l; 

if  (time_out>4000000) 


10  key  pressed  in  time,  */ 
return  -1  */ 


goto  menu_bail; 


time_out=0; 

b- { char ) getchar { ) ; 

if  tb>='a')  b-=32; 


if 

(»" 

1) 
=  0) 

l(b>='0' 
((b>-'A' 

l(b=='     ' 
(b==S)  |  | 
(b==10) 
goto  check 

)&Mb< 
)&&(b<- 
)&&(a= 

_buf ; 

{ 

m= ( int)b; 

goto  menu_bail; 

if 

Cb= 

-») 

if    (n!=0) 

I*    reset  counter  */ 
/*  pull  character  from  buffer  */ 
/*  convert  lower  case  to  upper  */ 
/*  check  validity  */ 
=  '9'))|| 

=,Z')ii(a==0)) | | 
=0)  1  |  | 


/*  if  character  invalid, 
if  non-numeric  data,  */ 
/*  return  character  */ 


check  for  backspace  */ 
/*  avoid  divide  by  zero 


get  new  one  */ 


n/=10; 

display_n(a, d,n) ; 
J 

goto  check_buf; 

} 
if  (b==10)  /'  check  for  CR 

{ 

m=n; 

goto  menu_bail; 

} 
v=n*10+( (int) (b-'O* ) ) ; 
if  (  (a==4)&&(v>9999) )  goto  check_buf; 
if  ( (a==3)&&(v>999) )  goto  check_buf ; 
if  ( (a==2)&&(v>99) )  goto  check_buf; 
if  (  (a=  =  l)&&(v>9) }  goto  check_buf; 
n-v; 

display_n(a,d, n) ; 
goto  check_buf; 


menu_bail : 

return tm) ; 
} 


Nothing  exception 


/*  nothing  is  the  exception  which  executes  on  an  interrupt  from  either  TPU 
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TPU  10,  or  TPU  11.  This  exception  takes  the  place  of  flexc,  f2exc,  and 
f3exc  both  on  initialization  and  during  the  time  when  rms  and  mean  data 
is  being  calculated.  */ 

void  nothing(void) 

{ 
#asm 

;noexc  is  the  entry  point  for  the  nothing  exception 
noexc    movem.l  d0-d7/a0-a6, - (a7)      ; Push  registers  to  stack 


btst.b  #$03,SFFFE20 
beq  nfll: 

bclr.b  #$03,$FFFE20 
bra  nobail: 


Test  if  interrupt ' s  from  TPU  11 
Skip  next  portion  if  not  from  TPU  11 
Acknowledge  and  turn  off  interrupt 
Exit  exception 


nfll     btst.b  #$02,$FFFE20        ,-Test  if  interrupt '  s  from  TPU  10 

beq  nflO:  ; Skip  next  portion  if  not  from  TPU  10 

bclr.b  #$02,$FFFE20        ;Acknowledge  and  turn  off  interrupt 
bra  nobail:  ;Exit  exception 

; Interrupt  must  be  from  TPU  9,  so  restore  sampling  exceptions 
nflO     bclr.b  #$01,$FFFE20        .-Acknowledge  and  turn  off  interrupt 


Look  at  nupc  valu 

If  not  zero,  skip  next  portion 

nupc-prof_period 


move . 1  _nupc , dO 
bne  decnupc: 
move.l  _prof_period, dO 
decnupc   subq . 1  #  $0 1 , dO  ; nupc=nupc - 1 

move.l  d0,_nupc  ;put  back  in  nupc 

cmp.l  _npa,d0  ; Compare  nupc  to  npa 

bne  nobail:  ;Exit  exception  if  nupc!=npa 

fit  nupc=npa,  restore  fl,  12,    &  f3  exceptions. 

lea  f lentry, aO 

lea  f 2entry, al 

lea  f 3entry,a2 

move.l  aO,$092C 

move.l  al,$0928 

move.l  a2, $0924 


Get  address  of  11 's  sampling  exception 
Get  address  of  10 's  sampling  exception 
Get  address  of  9's  sampling  exception 
Restore  TPU  11 's  exception 
Restore  TPU  10 's  exception 
Restore  TPU  9's  exception 


nobail   movem.l  {a7) +, d0-d7/a0-a6   ;Pull  registers  from  stack 
rte 

(tendasm 


Offload  routine 


void  of f load (void) 
{ 

bool  speed; 
char  command [20]; 
int  m,n; 

short  read_file,- 
unsigned  char  *header_byte; 
unsigned  long  tmu,- 
unsigned  short  "disk_header; 

printf (CLEAR) ; 

printf (UNDON) ; 

printf ( "Of f load  data  files."); 

printf (UNDOFF) ; 

printf (" \n\nRETURN  -  offload  at  current  baud  rate"); 
printf ("\nl  -  High  speed  offload  (115200) \n\n" ) ; 
printf ( "%c  ") ; 

printf  (deselection  =>•' ,13); 
get_speed: 

m=menu_input (0,0,0) ; 

if  (m==(-l))  goto  of f load_bail; 

if  { ! ( (m==(10) ) | | (m==(49) ) ) )  goto  get_speed; 

if  (m== (10) )  speed=0; 

else  speed=l; 
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if  (speed"l) 


printf  ( *' \n\nChanging  baud  rate  to  115200  -  change  terminal  baud  rate"); 
printf!"  and  press  SPACE"); 
SiraSetFSys(14720000) ; 
SerSetBaud(115200, 0) ; 
verify_change : 

m=menu_input (0,0,0)  ; 

if  (m=={-l)) 

SimSetFSys (16000000)  ; 

SerSetBaud(38400, 0) ; 

goto  of f load_bail; 

} 
if  (m!~32)  goto  verify_change; 
} 

printf (" \n\nNext  data  file  «  %d" , data_f ile_pointer) ; 
printf (" \nEnter  0  to  return  to  main  menu."); 
printf (" \nOff load  which  data  file?\n"); 
get_of f load: 

if  (data_f ile_po inter ==1) 

n*l; 

else 

n-data_f ile_po inter- 1; 
m=menu_input ( 3  ,  0 ,  n )  ; 

if  ((bss(-I)) | |  {ra«s0)J  goto  offload_bail; 
if  (m>max_f iles)  goto  get_offload; 
read_f ile- { short) m; 


of f_loop: 

printf {" \n\nReading  file..."); 

if  ( (de-DFRead( (ptr) sqbase, read_f ile) ) ! =0) 

printf (" \nDrive  error  ->  %d\n",de); 
else 

printf ( " \nSuccess . \n" ) r 

DriveOff  0  ,- 

/*  Determine  data  file  size  */ 
disk_header= (unsigned  short  *) sqbase; 
printf ( "\nnfp  =  %ld  ", 

nfp- (unsigned  long) ( *disk_header++) ) ; 
printf ("npp  -  %ld  ", 

npp= (unsigned  long) ( *disk_header++) ) ; 
printf ("nosp  -  %ld  ", 

nosp= (unsigned  long) (*disk_header++) ) ; 
printf ( "terap_mta_numscans  =  %ld  ", 

temp_mta_numscans= (unsigned  long) (*disk_header++) ) ; 
of f start- (unsigned  char  *) (sqbase) + (24*npp+8*nosp) +32768; 
tmu-nfp* (6*npp+2*nosp) + (152*temp_mta_numscans) +2  6*9+24; 
offlen=( (tmu/128)*128+128) ; 
of f end=of f start+of f len-1; 

printf ( " \nAddress  of  start  -  %p" , of f start) ; 
printf (" \nLength  of  file  =  %lu  bytes" , of flen) ; 
printf { " \n\nBegin  transfer . " ) ; 

XmodemSendMem(off start, of flen, 5) ; 

read_f ile+=l ; 

if  (read_f ile>max_f iles)  goto  of f load_bail; 

delay„no_clock( )  ,- 

printf { "\n\nOff load  next  file?  (data  file  #%03d)  -  enter  Y  or  N. 

, read_f ile) ; 
get_next : 

m=menu_input (0,0,0); 

if  (m==  (-1)  )  goto  f  ix„baud,- 

if  (! ( (m== (78) ) | | (m==(89))))  goto  get_next; 

if  (m=={89))  goto  off_loop; 
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f ix_baud: 

if  (speed--!) 


SimSetFSys(14720000)  ,- 

SerSetBaud{115200, 0) ; 

printf (CLEAR) ; 

printf (" \n\nChanging  baud  to  38400  -  change  terminal  baud  rate" 

printf {"  and  press  SPACE"); 

SimSetFSys (16000000)  j 

SerSetBaud(38400, 0)  ; 

veri  f y_change_back : 

m=menu_input (0,0,0); 

if  (m=={-in  goto  of fload_bail; 

if  (m!=(32) )  goto  verify_change_back; 


} 
of f load_bail : 
printf (CLEAR) ; 


Reset  recovery  routine 


void  ret_from„reset (void) 
{ 

printf (" \nReset  occurred  -  continuing  pre-reset  execution\n" ) ; 
} 


Time  and  date  set  routine 


void  set_time (void) 
{ 

int  m,n; 

struct  tm  time_struct, *vtime_struct; 
time_t  timer; 

timer=RtcToCtm( ) ; 
vtime_struct=localtime (&timer) ; 

printf (CLEAR) ; 

printf (UNDON) ; 

printf ("Set  time  and  date"); 

printf (UNDOFF) ; 

printf ("\n\nMonth  (l-12)\n"); 
get_month: 

n= (vtime_struct->tm_mon) +1; 

m=menu_input (2 , 0,n) ; 

if  (m== (-1) )  goto  time_bail; 

m--l; 

if  ( (m<0) | | (m>ll) )  goto  get_month; 

time_struct . tm_mon=m; 

printf ("\n\nDay  (l-31)\n"); 
get_day: 

n=vtime_struct->tm_mday; 

m=menu_input (2, 0,n) ; 

if  (m==(-l))  goto  time_bail; 

if  ( (m<l) | | (m>31) )  goto  get_day; 

time_struct . tm_mday=m; 

printf ("\n\nVear  (0-99) \n") ; 
get_year : 

n-vtime_struct->tm_year; 

m-menu_input ( 2 , 0 , n ) ; 

if  (m==(-l))  goto  time_bail; 

if  (m<4)  goto  get_year; 
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time_struct. tm_year=m; 

printf ( " \n\nHour  (0-23)\n">; 
get_hour : 

n=vtime_struct->tm_hour ; 

m=menu_input (2 , 0,n) ; 

if  (m==(-l})  goto  time_bail; 

if  (m>23)  goto  get_hourr 

time_struct . tm_hour=m; 

printf ("\n\nMinute  (0-59) \n") ; 
get_rainute: 

n=vtime_struct->tm_min; 

ra -me nu_ input { 2 , 0 , n) ; 

if  (m==(-l))  goto  time_bail; 

if  (m>59)  goto  get_minute; 

time_struct . tm_min=m; 

printf ("\n\nSeconds  {0-59) \n" ) ; 
get_second: 

n=vtime_struct->tm„sec; 

m=menu_input (2 , 0,n) ; 

if  (m==(-l))  goto  time_bail; 

if  (m>59)  goto  get_second; 

time_struct .  tm_sec=in; 

SetTimeTM(&tirae_struct,NULL) ; 
time_bail : 

printf (CLEAR) ; 
} 


*  Main  data  collection  routine 

void  start_collection (void) 
{ 

float  totaLsecs; 
int  m,n; 

int  total_time_hours ; 
int  total_time_jnins; 
int  total_time_secs ; 
int  total_time_hunds; 
unsigned  short  *disk_header; 
unsigned  long  testtime; 
int  counter; 

struct  tm  time_struct, *vtime_struct; 
time_t  timer; 

printf (CLEAR) j 

if  (samp__mode==0)  f irst_time=l; 

else 

{ 

num_to_fill=num_col lections ; 
timer=RtcToCtm( ) ; 
if  (start_time< (timer+30) ) 
{ 

vtime_struct=localtime (&timer) ; 
t  ime_s  true  t . tm_year- vt ime_s  true t ->  tm_y ear ; 
time_struct . tm_mon=vtime_struct->tm_mon; 
time_struct . tm_mday=vtime_struct->tm_mday,- 
time_struct . tm_hour=0; 
time_struct .  tm_min=0  ,- 
time_struct . tm_sec=0; 
testtime=0; 
counter=0; 

while ( testtime< (timer+3  0) ) 
{ 

tes ttime=mktime{ Sit ime_s true t)  +  (of  fset+hours_between* counter)  * 
36; 
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star t_time= test time; 
counter+=l; 


} 

collection_loop: 

i f  { sarap_raode- - 1 ) 

( 

printf ( "Next  collection  at  "); 
vtime_struct-localtime  (&start_time)  ,- 
printf  ( "%02d:%02d:%02d.OO" , vtime_struct->tm_hour, vtime_struct->tm_min, 

vtime_struct->tm_sec) ; 
printf ( " \n\n*  to  interrupt ..."); 
while  (timer< (start_time-30) ) 
{ 

timer=RtcToCtm( ) ; 
sleep  (GetTickRate  ( )  )  ,- 
if  {kbhitO  1*0) 
{ 

if  ( (char)getchar ()=='*') 
{ 

nura_to_f  i 1 1 =0 ; 
printf (CLEAR) ; 
goto  sample_bail; 
} 
} 
} 
an_power (1) ; 
delay_no_clock( ) ; 
start_mta ( ) ; 
} 

printf (CLEAR) ; 

printf  (UNDON)  ,- 

printf ( "Start  data  collection"); 

printf (UNDOFF) ; 

if  (data_f ile_pointer>max_f iles) 
{ 

printf (CLEAR) ; 
printf (FLSHON) ; 

printf {" \n\nNo  free  data  files!"); 
printf (FLSHOFF) ; 

printf (" \nOff load  data  and  reset  file  pointer."); 
num_to_fill-0; 
goto  sample_bail; 
} 

npp=perm_npp ; 

npa-perm_npa  ,- 

prof_period=perra_prof_period; 

nfpc=mem_nfp; 

total_nfp=shore_nfp; 

mem_setup( ) ; 

total_secs-nfp*prof_period/100. 0; 
total_time__hours=total_secs/3600; 

total_time_mins=(total_secs-3  600.Q*total_time_hours) /60.0; 
total_time_secs= (total_secs-3  600 .O*total_time_hours 

-60. O*total_time_mins) ; 
total_time_hunds= (total_secs- (int) total_secs) *100 . 0; 
sprint f (acquisition_time, "%02d: %02d: %02d.%02d" , total_time_hours 

, total_time_mins, total_time_secs, total_time_hunds) ; 


printf("\n*  Every  %4.2f  seconds,  one  rms  profile,  computed  from  " 

,prof_period/100.0)  ,- 
printf ("%ld  profiles,  will  be",npa); 
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printf ("\n   collected  from  the  ACP's."); 

printf ("\n*  Each  ACP  profile  will  contain  "); 

printf("%ld  points ." ,perra_npp) ; 

printf ("\n*  %ld  single  point  instruments  will  be  sampled. ", nosp) ; 

printf  C\n*  Every  %4.2f  seconds,  one  mean  value  from  each  single 

,prof_period/100.0) ; 
printf ( "point  instrument,  " ) ; 

printf("\n  computed  from  %ld  values,  will  be  collected. ", npa) j 
printf (" \n\nTime  to  fill  this  data  file  ->  %s" , acquisition_time} ; 
printf (" SnTime  to  completion  of  sampling  ->  %d  hours" ,perm_sdt) ; 

printf (" \n\nlnitializing  memory. .  .  " )  ; 
mat_headers ( ) ; 
clear_mem( ) ; 

if  ( (f irst_time--l) && (samp_mode--0) ) 
{ 

printf (" \nEnter  0  files  to  return  to  main  menu"); 
printf (" \nFill  how  many  files?  (l-%d)\n" 

,max_f iles-data_f ile_pointer+l) ; 
get_num_f iles : 
n-lj 

m=menu_input (3 , 0,n} ; 
if  (  (m«(-l))  |  |  fm==0)  ) 
{ 

printf (CLEAR) ; 
goto  sample_bail; 
} 
if  (m>  (max__f  iles-data_f  ile_pointer)  )  goto  get_num_f  iles,- 
num_t  o_  f  i  1 1  -m  ,- 

an_power { 1 ) ; 

delay_no_clock( } ; 

start_mta ( ) ; 

timer=RtcToCtm( ) ; 

s tar t_t ime- 1  imer  +  15; 

printf (CLEAR) ; 

printf ( "Warming  up. . . \n" ) ; 


if  ( (samp_mode==0) && (f irst_time==0) ) 
{ 

start_mta ( ) ; 
timer=RtcToCtm ( ) ; 
start_time=timer+2; 
} 

if  ( ( first_time==l) | | (samp_mode=-l) ) 
{ 
printf (" \nTurning  on  instrument  power,  setting  full  charge..."); 

PClear{E,2) ; 

e2  =  0; 

printf  ("\nTPU  buffers  enabled."),- 

f irst_time=0; 

} 

/*  Write  starting  time  to  data  file  */ 
t imer - s tar t_time; 
vtime_struct-localtime (&timer) ; 

( *tmbas++) ■ (unsigned  short) vt irae„struct-> tm_year; 
( *tmbas++) - (unsigned  short) (vtime_struct->tm_mon) +1; 
( *tmbas++) - (unsigned  short) vtime_struct->tm_mday; 
( *tmbas++) - (unsigned  short) vtime_struct-> tm_hour; 
{ *tmbas++) = (unsigned  short )vtime_s true t->tm_min, ■ 
(*tmbas++) - (unsigned  short )vtime_s true t->tm_sec; 

printf (" \n\nData  file  %d  will  begin  filling  at  " , data_f ile_pointer) ; 
vtime_s  true  t= local  time  (&start_time)  ,- 

printf ("%02d:%O2d:%O2d.O0", vtime_struct->tm_hour, vtime_struct->tm_minJ 
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vtime_struct->tm_sec) ; 

delay_no_clock ( ) ; 

/*  Grab  parallel  i/o  port  address  */ 

PAR_ENABLE ; 

pioport-PPORT;       /*  Assigns  a  pointer  to  the  port 

/*  Initial  time  header  to  parallel  port  */ 

*pioport=OxFF; 

*pioport=0x00; 

*pioport= (unsigned  char) vtime_s true t->tm_year; 

*pioport= (unsigned  char) (vtirae_struct->tm_mon) +1; 

*pioport- {unsigned  char) vtime_s true t->tm_mday, ■ 

*pioport- (unsigned  char) vtime_struct->tm_hour; 

*pioport- (unsigned  char) vtime_s true t->tm_min; 

*pioport= (unsigned  char) vtime_s true t->tm_sec; 
*pioport= (unsigned  char ) prof _period; 
*pioport- (unsigned  char)npp; 
*pioport- (unsigned  char)nosp; 
*pioport- (unsigned  char) (nfp/256); 
*pioport- (unsigned  char) (nfp%256); 
*pioport= (unsigned  char) (shore_mta_numscans/256) ; 
*pioport= (unsigned  char) (shore_mta_numscans%256) ; 

prof_index-l;        /*  Initial  shore  data  header  index  */ 

wait_for„start: 

timer=RtcToCtm() ; 

if  (start_time! -timer)  goto  wait_£or_start; 

/*  Run  main  sampling  routine  */ 

main_sample ( ) ; 

/*  Write  ending  time  to  file  */ 
timer=RtcToCtm() ; 
vtime_struct=localtime(&timer) ; 

( *tmbas++) - (unsigned  short) vtirae_struct->tm_year ; 
( *tmbas++) - (unsigned  short) (vtime_struct->tm_mon) +1; 
( *tmbas++) - (unsigned  short) vtime_struct->tm_mday ; r 
(*tmbas++) ■ (unsigned  short) vtime_struct- >tm_hour ; 
(*tmbas++) = (unsigned  short) vtime_struct->tm_min; 
(*tmbas++) = (unsigned  short) vtirae_struct-> tm_sec; 

/*  Final  time  header  to  parallel  port  */ 

*pioport=OxFF; 
*pioport=OxFE; 

*pioport= (unsigned  char) vtime_struct->tm__year; 

*pioport- (unsigned  char)  (vtime_struct->tm_mon)  +1  ,- 

*pioport- (unsigned  char) vtime_struct->tm_mday; 

*pioport- (unsigned  char) vtime_struct->tm_hour; 

*pioport= (unsigned  char) vtime„struct->tm_min; 

*pioport= (unsigned  char) vtime_struct->tm_sec; 

/*  Disable  parallel  port 

PAR_DISABLE; 

/*  Get  data  from  mta  if  TT7  data  file  only!  */ 
if  (dif_nfp==l) 

collect_mta ( ) ; 
else 

talk_thru() ; 

/*         talk_thru() ;        •/ 
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/*  Save  to  disk  */ 

printf ( "\nSaving  data  file. .  .  " )  ; 
disk_header= (unsigned  short  *) sqbase; 
*disk_header++- (unsigned  short) nfp; 
*disk_header++= (unsigned  short )npp, ■ 
*disk_header++= (unsigned  short) nosp; 
*disk_header- (unsigned  short )mta_numscans; 

printf ( "\nnfp  =  %ld  npp  =  %ld  nosp  ■  %ld  mta_numscans  =  %ld\n" , 
nfp, npp, nosp, mta_numscans) ; 

if  ( (de=DFWrite( (ptr) sqbase, (short) data„file_pointer) ) !-0) 

printf (" SnDrive  error  ->  %d\n",de); 
else 

printf  ( "  \nSuccess  .  \n" )  ,- 

DriveOff {) ; 
data_£ile_pointer+=l; 

num_ t o_ f i 1 1 - - 1 ; 

if  ( (num_to_f ill>0) && (samp_roode=-l) ) 
{ 

timer-RtcToCtmf ) ; 
vtime_struct=localtime (stimer)  j 
time_struct .  tm_year=vtime_struct->tm_year,- 
time_struct . tm_mon=vtime_struct->tm_mon; 
time^struct . tm_mday=vtime_struct->tm_mday; 
tirae_struct . tm_hour-0; 
time_struct . tm„min=0; 
time_struct. tm_sec=0; 
testtime=0; 
counter-0; 
while! test time< (timer+30) ) 

( 

testtime-mktime(&time_struct) + (of fset+hours_between*counter) * 
3  6; 

s tar t_time- test time; 

counter +"1 j 

> 

printf (CLEAR) ; 

an_power (2) ; 

PSet(E,2) ; 

printf ("Power  off,  TPU  buffers  disabled. \n" ) ; 

} 

if  (num_to_f ill>0)  goto  collection_loop; 

printf (CLEAR) ; 
sample_bail: 

an_power ( 2 )  ; 

printf ( "Analog  power  turned  off  -  charge  set  to  floating."); 

printf  ( "\nTPU  buffers  disabled."),- 

PSet(E,2) ; 
} 

*  MTA  setup 

void  start_mta (void) 
{ 

char  mta_id; 
char  comp_string[5] ; 
int  i ,  j  ,- 
long  timeout; 

struct  tm  time_struct, *vtime_struct; 
time_t  timer, tplus4; 

/*  Initialize  TPU  communications  */ 

printf (CLEAR) ; 
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printf ("Talking  to  MTA's"); 
printf ( " \n\nPausing. . . " ) ; 

PClear(E,2);  /*  Enable  TPU  buffers  */ 

i=TSerOpen(inchan, INP, 100, 9600, 'N' , 8, 1) ;     /*  input  channel  */ 
printf ( "\nOpen  channel  %d  -  status  ■  %d,  INP  =  %d\n" , inchan, i, INP) ; 
i=TSerOpen(outchan,OUTP, 100, 9600, 'N' , 8, 1) ;   /*  output  channel  */ 
printf ("Open  channel  %d  -  status  *  %d  OUTP  =  %d\n\n" , outchan, i.OUTP) ; 

/*  Program  HTA  -  Rl  command  */ 

delay_no_clock ( ) ; 

for  (j=0;j<3; j++) 
{ 

mta_id- (char)  j+4  9  ; 
printf ("MTA  #%c\n" ,mta_id) ; 

sprint f(mta_prog, "P%02dB%O3dT%04dG%O3dS%051dD%04d%c" ,mta_pulse, 

mta_blank,mta_thresh[ j+1] ,mta_gain, shore_mta_numscans,mta_delay, 13)  ; 

prog_request: 

printf ( "Sending  program  request  (%cRl) . . . \n" ,mta_id) ; 

TSerPutByte (outchan, 0) ; 

TSerPutByte (outchan, mta_id) ; 

TSerPutByte (outchan, ' R' ) ; 

TSerPutByte (outchan, ' 1' ) ; 

TSerPutByte (outchan, 13 ) ; 

printf  ( "Awaiting  acknowledge  (%cAl)  .  .  .  \n"  ,mta_id)  ,- 
grabstring(mta„id) ; 

if  (stringerror ! -0) 

goto  start_mta_bail; 
sprint f (comp_string, "%cAl" ,mta_id) ; 
if  (strcmpf  instring, comp_string) ! =0) 
{ 

sprint f  (comp_string,  "%cK"  ,mta_id)  ,- 
if  (strcmp( instring, comp_string) ==0) 
{ 

printf ( "Command  not  understood,  resending. . . \n" ) ; 
timeout=0; 
while (timeout<2 00000) 

timeout+=l; 
goto  prog_request; 
} 
goto  start_rata_baii,- 
} 

printf ( "Sending  program  string  -  %c%s\n" ,mta_id,mta_prog) ; 
TSerPutByte ( outchan, 0 ) ; 
TSerPutByte  (outchan,  mta_id)  ,- 
for  (i  =  0;  i<28,-i  +  +) 

TSerPutByte  ( outchan,  mta_prog  [i  ]  )  ,- 

printf ( "Awaiting  acknowledge  (%cA) . . . \n" ,mta_id) ; 
grabstring{mta_id) ; 
if  (stringerror ! -0) 

goto  start_jnta_bail; 
sprintf {comp_string, "%cA" ,mta_id) ; 
if  (strcmp{ instring, comp_string) !=0) 
{ 

sprintf (comp_string, "%cN" ,mta_id) ; 
if  (strcmp( instring, comp_string) ==0) 
{ 

printf ( "Command  not  understood,  resending. .. \n" ) ; 
timeout=0; 
while (timeout<2 00000) 

timeout+=l; 
goto  prog_request; 
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/*  Set  MTA  time 
sendtime 


>tm_year , 
> tm_hour , 


goto  start_mta_bail; 
} 

R3  command  */ 

tplus4=RtcToCtm  0+4; 
vtime_struct=localtime (&tplus4) ; 

sprintf fmta_time, "S%02d%02d%02d%02d%02d%02d%c" , vtime_struct- 

vtime_struct->tm_mon, vtime_s true t->tm_mday,vt ime_s true t- 

vtime_s true t- >tm_min, vtime_s true t->tm_sec, 13)  ; 

printf ( "Sending  time-sync  request  (%cR3 ) . . . \n" ,mta_id) ; 

TSerPutBy te ( outchan ,  0 )  ; 

TSerPutByte(outchan,mta_id) ; 

TSerPutByte (outchan, 'R'  )  ; 

TSerPutByte ( outchan , '  3  '  )  ; 

TSerPutByte (outchan, 13) ; 

printf ( "Awaiting  acknowledge  (%cA3) . . . \n" ,mta_id) ; 
grabstring (mta„id) ; 
if  (stringerror ! =0) 

goto  start_mta_bail; 
sprintf  (comp__string,  "%cA3  "  ,mta_id)  ; 
if  (stremp (instring, comp_string) ! -0) 
{ 

sprintf (comp_string, "%cN" ,mta_id) ; 
if  (stremp (instring, comp_string) --0) 
{ 

printf [ "Command  not  understood,  resending. . . \n" ) ; 
timeout=0; 
while (timeout<2 00000) 

timeout+-l; 
goto  sendtime; 
} 
goto  start_mta„bail; 
} 


printf ( "Sending  present  time  +5  seconds  -  %c%s\n" ,mta_id,mta_time) ; 
TSerPutBytefoutchan, 0) ; 
TSerPutByte(outchan,mta_id) ; 
for  (i^0;i<14;i++) 

TSerPutByte(outchan,mta_time[i] } ; 

printf ( "Awaiting  acknowledge  (%cA) . . . \n" ,mta_id) ; 
grabstring (mta_id) ; 
if  (stringerror! -0) 

goto  start_mta_bail; 
sprintf (comp_string, "%cA" ,mta_id) ; 
if  (stremp (instring, comp_string} ! =0) 
{ 

sprintf  (comp_string,  "%cN"  ,mta_id)  ,- 
if  (stremp (instring, comp_string) ==0) 
{ 

printf ( "Command  not  understood,  resending. .. \n" ) ; 
timeout=0; 
while (timeout<2 00000) 

timeout+=l; 
goto  sendtime ; 
} 
goto  start_mta_bail; 
} 

/*  Wait  for  tplus4  time  to  arrive,  then  signal  MTA  with  CR  */ 
printf ( "Awaiting  sync  time. . . \n" ) ; 
while{tplus4>RtcToCtm{) ) ; 
TSerPutByte(outchan, 13) ;       /*  Synchronizing  CR  */ 
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printffSync  pulse  sent.   MTA  #%c  should  be  started.  \n"  ,mta_id)  ; 


start_mta_bail: 
printf ( "\n" ) ; 


TSerClose(inchan) ; 
TSerClose {outchan) ; 


/*  Disable  TPU  buffers 


Storm  mode 


void  storm_mode (void) 
{ 

struct  tm  time_struct, *vtime_struct; 
time_t  timer; 

samp_mode=l; 

f irst_time=l; 

timer=RtcToCtm()  ,- 

vtime_struct-localtimG (&timer) ; 

time_struct . tm_year=vtime_struct->tm_year; 

time_struct. tm_mon-vtime_struct->tm_mon; 

time_struct . tm_mday=vtime_struct->tm_mday; 

time_struct. tm_hour=0; 

time_struct . tm_min-0 ; 

time_struct . tm_sec=0; 

start_time-mktime (&time_struct ) + ( {RtcToCtm( ) -mktime (&time_struct) ) 

/ (3  600*hours_between) ) * (3600*hours_between) 

+3600*hours_between; 
start_collection{ ) ; 
} 


Talk-through  portion  to  communicate  directly  to  MTA's  from  shore 


void  talk_thru (void) 
1 

char  ch=0; 
int  i ,  j  ,- 

PConfOutp{E,2) ; 
PClear(E,2} ; 


Enable  TPU  buffers  */ 


PConf0utp{F,7) ; 
PSet (F,7) ; 


£7-1; 


MTA  profiler  on 


printf (CLEAR) ; 
printf ("MTA  Talkthru\n\n" ) ; 
printf ("TPU  buffers  enabled,  MTA  power  on.\n"l; 
printf ( "Hit  !  to  exit,  use  -  for  null.Vn"); 

i=TSerOpen(inchan, INP, 100, 9600, 'N', 8,1} ; 
printf ("Open  channel  8  -  status  =  %d,  INP 
i=TSerOpen (outchan, OUTP, 100,  9600,  'N' ,8, 1) ; 
printf ( "Open  channel  9  -  status  =  %d  OUTP 


/*  input  channel  */ 
%d\n" , i, INP) ; 
/*  output  channel  */ 
%d\n\n",i,OUTP) ; 


whilefch! 


!  ') 


while(kbhit{)==0) 
( 
if  (TSerByteAvail (inchan) ) 
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ch=TSerGetByte(inchan) ; 


SerPutByte(ch) ; 


ch=getchar ( ) ; 

if  (ch=='-')  ch=0,- 

if  (ch==10)  ch=13; 

TSerPutByte {outchan, ch} ; 


TSerClose ( inchan) ; 
TSerClose (outchan) ; 

PSet(E,2);        /*  Disable  TPU  buffers  */ 
e2  =  l; 

printfC'TPU  buffers  disabled,  MTA  power  off.\n"); 
} 


TPU  and  interrupt  setup  routine 


void  tpu_setup(void) 

{ 
#asm 

movem . 1 


; Push  registers  to  stack 


; Enable  interrupts  and  set  exception  addresses  to  nothing  routine 


move.w  #$2400, sr 
lea  noexc,a0 
move.l  a0,$0924 
move.l  aO, $0928 
move.l  a0,$092C 


Mask  all  interrupts  (+  set  supervisor  mode) 
Start  of  nothing  exception 
Set  TPU  9's  exception  address 
Set  TPU  10 's  exception  address 
Set  TPU  11' s  exception  address 


;Turn  off  channels 
move.w  #0,$FFFE1C 
move.w  #0,$FFFE1E 


,-Clear  chs.  8-15  priorities  (CPRO) 
;Clear  chs.  0-7  priorities  (CPRl) 


.-Interrupt  set  up 
move.w  #$0000,$FFFEOA 


/Function  select  registers  (CFSR0,  CFSR1,  and  CFSR2) 


move.w  $FFFE0C,d4 
andi.w  #$F0OO,d4 
ori.w  #$0777, d4 
move.w  d4, $FFFE0C 


Get  current  functions  from  12-15 
Clear  functions  from  12-14 
Set  chs.  12-14  to  SPWM 
Put  back  in  CFSR0 


move.w  $FFFE0E,d4 
andi.w  #$0O0F,d4 
ori.w  #$7770, d4 
move.w  d4,$FFFE0E 


Get  current  function  from  8-11 
Clear  functions  from  9-11 
Set  chs.  9-11  to  SPWM 
Put  back  in  CFSR1 


;Host  sequence  registers  (HSQR0  &  HSQR1) 


move.w  $FFFE14,d4 
andi.w  #$C003,d4 
ori.w  #$0008, d4 
move.w  d4, $FFFE14 


Get  current  values  for  8-15 
Clear  bits  2-13  (chs.  9-14) 
Ch.  9  mode  2  chs.  10-14  mode  0 
Put  back  in  HSQR0 


; Parameter  RAM  values 


F3  software  trigger- 


move 
move 
move 
move 
move 


w  #$0092,$FFFF90 
w  #$0005,$FFFF94 
w  #$2710,$FFFF96 
W  #$A590,$FFFF98 
w  #$0000, SFFFF9A 


Channel  control 
5  microsecond  hightime 
100  Hz  frequency 
Links  and  ref_addrl 
No  delay 


;Channel  10  -  F2  software  trigger 
move.w  #$0O92,$FFFFA0    ,-Channel  control 
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move.w  #$0005, $FFFFA4 
move.w  #$2710,$FFFFA6 
move.w  #$0090,$FFFFA8 
move.w  #S1D4C,$FFFFAA 


;5  microsecond  hightime 
;100  Hz  frequency 
;Ref_addrl  to  ch.  9's  lastris 
;7.5  ms  delay 


;Channel  11  -  Fl  software  trigger 


move 
move 


w  #$0092,$FFFFBO 
w  #$0005,$FFFFB4 
W  #$2710, $FFFFB6 
W  #$0090,$FFFFB8 
w  #$1388,$FFFFBA 


Channel  control 

5  microsecond  hightime 

100  Hz  frequency 

Ref_addrl  to  ch.  9's  lastrise 

5  ms  delay 


; Channe 1  1 2  -  sync  channe 1 


move.w  #$OO92,$FFFFC0 
move.w  #$0014, $FFFFC4 
move.w  #$2710,$FFFFC6 
move.w  #$0090,$FFFFC8 
,-move.w  #$0F82,  SFFFFCA 
move.w  #$05BE, SFFFFCA 


Channel  control 
20  microsecond  hightime 
100  Hz  frequency 
Ref^addrl  to  ch.  9's  lastrise 
,■3.97  ms  delay 
r 1 .47  ms  delay 


; Channel  13  -  hardware  trigger  channel 


move 
move 
move 
move 
move 


W  #$0092,$FFFFDO 
W  #$0014,$FFFFD4 
w  #$09C4,$FFFFD6 
w  #$0090,$FFFFD8 
W  #$09A6, $FFFFDA 


Channel  control 

20  microsecond  hightime 

400  Hz  frequency 

Ref_addrl  to  ch.  9's  lastrise 

2.47  ms  delay 


; Channel  14 


move 
move 
move 
move 


W  #$0O92,$FFFFE0 
w  #$0007,$FFFFE4 
w  #$0O0A, $FFFFE6 
W  #$0090,$FFFFE8 
W  #$0000,$FFFFEA 


>D  trigger  channel 


Channel  control 

7  microsecond  hightime 

100  kHz  frequency 

Ref„addrl  to  ch.  9's  lastrise 

No  delay 


;Host  service  request  registers  (HSRR0) 


move.w  $FFFE18,d4 
andi.w  #$FFF3,d4 
ori.w  #$0008, d4 
move.w  d4,$FFFE18 
move.W  #$0OOC,$FFFElC 


,-Get  chs.  8-15 

; Clear  ch.  9 

,-  Initialize  ch.  9 

;Put  back  in  HSRR0 

;Set  ch.  9  priority  high 


statl  move.w  $FFFE18,d4     ;Pull  status  register 

andi.w  #$O0OC,d4      ;Look  at  just  ch.  9 
bne  statl:  ;Keep  looping  "til  0 


move.w  $FFFE18,d4 
andi.w  #$C0OF,d4 
ori.w  #$2AA0,d4 
move.w  d4,$FFFE18 
move.w  #$2AAE, $FFFElC 


Get  chs.  8-15 

Clear  chs.  10-14 

Initialize  chs.  10-14 

Put  back  in  HSRR0 

Ch.  9  high  priority,  10-14  middle 


stat2  move.w  $FFFE18,d4     ;Pull  status  register 

andi.w  #$3FF0,d4     ;Look  at  just  chs.  10-14 
bne  stat2:  ;Keep  looping  'til  0 

move.w  #$2000, sr       ;Unmask  interrupts 

movem.l  (a7)+,d0-d7/a0-a6     ;Pull  registers  from  stack 


void  voltage (void) 

{ 
#asm 


movem.l  d0-d7/a0-a6, - (a7)      ;Push  registers  to  stack 

bclr  #$03, $FFFFFC15      ;Set  port  2  on  multiplexer  board 
bset  #$04, $FFFFFC15 
bclr  #$05, $FFFFFC15 


move.l  #$50, dl; 


;Set  index  for  delay  purposes 
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vuwpa    move.b  $F8C00A,dO       ;Read  SFR  to  check  busy  flag 

bpl  vuwpa:  ;Stay  here  until  busy  flag  goes  low 

vuwpb    move.b  $F8C0OA,dO       ,-SFR  -  busy  goes  high  when  conversions  done 


bmi  vuwpb: 
move.b  #2,$F8C0OO 
subq.l  #$01, dl 
bne  vuwpa 


Stay  here  until  conversion  complete 
Set  channel  2  for  next  sample 
Decrement  index 
Loop  through  80  times 


; Channel  2  and  port  2  are  now  set  for  the  next  sample 
,-which  is  V (unregulated) 


vuloopa  move, 
vuloopb  move. 


$F8CO0A,dO       ;Read  SFR  to  check  busy  flag 

bpl  vuloopa:  ;Stay  here  until  busy  flag  goes  low 

$F8CO0A,d0       ;SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 

Get  MSB  from  A->D 

Get  LSB  from  A->D 

Set  channel  2  for  next  sample 

Bring  entire  word  into  dl 

store  result  in  variable  vunreg 


bmi  vuloopb: 

move 

b 

SF8C002, (a4) 

+ 

move 

b 

SF8CO00, (a4) 

+ 

move 

b 

#2,$F8C000 

move 

w 

-<a4) ,dl 

move 

w 

dl , _vunreg 

set  #$03, $FFFFFC15 


;Set  port  3  on  multiplexer  board 


move.l  #$50, dl  ,-Set  index  for  delay  purposes 

vrwpa    move.b  $F8C00A,dO       ,-Read  SFR  to  check  busy  flag 

bpl  vrwpa:  ;Stay  here  until  busy  flag  goes  low 

vrwpb    move.b  $F8C00A,dO       ;SFR  -  busy  goes  high  when  conversions  done 


bmi  vrwpb : 
move.b  #2, $F8CO00 
subq.l  #$01, dl 
bne  vrwpa 


Stay  here  until  conversion  complete 
Set  channel  2  for  next  sample 
Decrement  index 
Loop  through  80  times 


; Channel  2  and  port  3  are 
;which  is  V(regulated) 


set  for  the  next  sample 


vrloopa  move, 
vrloopb  move. 


$F8C00A,d0        ;Read  SFR  to  check  busy  flag 

bpl  vrloopa:  ;Stay  here  until  busy  flag  goes  low 

$F8C00A,d0       ;SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 

Get  MSB  from  A->D 

Get  LSB  from  A->D 

Set  channel  2  for  next  sample 

Bring  entire  word  into  dl 

store  result  in  variable  vreg 


bmi  vrloopb: 

move 

b 

SF8C002, 

a4) 

move 

b 

SF8COO0, 

a  4  ) 

move 

b 

#2,$F8COO0 

move 

w 

-(a4) ,dl 

move 

w 

dl , _vreg 

bclr  #$03,$FFFFFC15 
bclr  #$04,$FFFFFC15 
bclr  #$05,$FFFFFC15 


;Set  port  0  on  multiplexer  board 


move.l  #$50, dl;  ;Set  index  for  delay  purposes 

lklwpa   move.b  $F8CO0A,d0       ,-Read  SFR  to  check  busy  flag 

bpl  lklwpa:  ,-Stay  here  until  busy  flag  goes  low 

lklwpb   move.b  $F8CO0A,d0       ,-SFR  -  busy  goes  high  when  conversions  done 

bmi  lklwpb:  ,-Stay  here  until  conversion  complete 

move.b  #2,$FBCOO0  ;Set  channel  2  for  next  sample 

subq . 1  # $01 , dl  ; Decrement  index 

bne  lklwpa  ,-Loop  through  80  times 


; Channel  2  and  port  0  are  now 
;which  is  leak  detector  #1 


let  for  the  next  sample 


Iklloopa  move.b  $F8COOA,dO       ,-Read  SFR 
bpl  Iklloopa: 

lklloopb  move.b  $F8CO0A,d0        ;SFR 
bmi  lklloopb: 
move.b  $F8C002, (a4) 


move 
move 
move 
move 


b  $F8C0O0, (a 
b  #2,$F8C00O 
w  -{a4),dl 
w  dl,_volt_l 


check  busy  flag 

,-Stay  here  until  busy  flag  goes  low 
busy  goes  high  when  conversions  done 

,-Stay  here  until  conversion  complete 
;Get  MSB  from  A->D 
)+     ;Get  LSB  from  A->D 

;Set  channel  2  for  next  sample 
; Bring  entire  word  into  dl 
3kl    ,- store  result  in  variable  volt_leakl 
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bset  #$03,$FFFFFC15 
bclr  #$04,$FFFFFC15 
bclr  #$05,$FFFFFC15 


;Set  port  1  on  multiplexer  board 


move.l  #$50, dl;  ; Set  index  for  delay  purposes 

lk2wpa   move.b  $F8C00A,d0        ;Read  SFR  to  check  busy  flag 

bpl  lk2wpa:  ."Stay  here  until  busy  flag  goes  low 

lk2wpb   move.b  $F8CO0A,d0        ;SFR  -  busy  goes  high  when  conversions  done 

bmi  lk2wpb:  ;Stay  here  until  conversion  complete 

move.b  #2,$F8CO00  ,-Set  channel  2  for  next  sample 

subq.l  #$01, dl  ; Decrement  index 

bne  lk2wpa  ;Loop  through  80  times 


; Channe 1 
;  which  i: 


2  and  port  1  are  now  set  for  the  next  sample 
leak  detector  #2 


lk21oopa  move.b  $F8C0OA,d0       ;Read  SFR  to  check  busy  flag 

bpl  lk21oopa:  ;Stay  here  until  busy  flag  goes  low 

lk21oopb  move.b  $F8C0OA,dO       ; SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 

Get  MSB  from  A->D 

Get  LSB  from  A->D 

Set  channel  2  for  next  sample 

Bring  entire  word  into  dl 

store  result  in  variable  volt_leak2 


bmi  lk21oopb: 

move 

b 

$F8C002, 

a4)  + 

move 

b 

$F8C00O, 

a4)  + 

move 

b 

#2,$F8C00O 

move 

w 

-(a4),dl 

move 

w 

dl,_volt. 

_leak2 

bclr  #$03,$FFFFFC15 
bset  #$04, $FFFFFC15 
bset  #$05, $FFFFFC15 


; Set  port  6  on  multiplexer  board 


move.l  #$50, dl;         ;Set  index  for  delay  purposes 
tiwpa    move.b  $F8C00A,dO        ;Read  SFR  to  check  busy  flag 

bpl  tiwpa:  ,-Stay  here  until  busy  flag  goes  low 

tiwpb    move.b  $F8C00A,dO        ;SFR  -  busy  goes  high  when  conversions  done 


bmi  tiwpb: 
move.b  #3,$F8CO00 
subq.l  #$01, dl 
bne  tiwpa 


Stay  here  until  conversion  complete 
Set  channel  3  for  next  sample 
Decrement  index 
Loop  through  80  times 


/Channel  3  and  port  6  are  now  set  for  the  next  sample 
,-which  is  Internal  temperature 


tiloopa  move.b  $F8C0OA,d0       ;Read  SFR  to  check  busy  flag 

bpl  tiloopa:  ;Stay  here  until  busy  flag  goes  low 

tiloopb  move.b  $F8C00A,d0        ;SFR  -  busy  goes  high  when  conversions  done 

Stay  here  until  conversion  complete 

Get  MSB  from  A->D 

Get  LSB  from  A->D 

Set  channel  3  for  next  sample 

Bring  entire  word  into  dl 

store  result  in  variable  tempin 


bmi  ti 

oopb: 

move 

b 

$F8C002, 

(a4)  + 

move 

b 

$F8C00O, 

(a4M 

move 

b 

#3,$F8C00O 

move 

w 

-(a4) ,dl 

move 

w 

dl , _temp 

in 

FFFC1 

5 

;Set 

port 

5  on  multiplexer  board 


bclr  #$04, $FFFFFC15 

bset  #$05,$FFFFFC15 


move.l  #$50, dl; 


vbwpa 
vbwpb 


move.b  $F8C00A,d0 
bpl  vbwpa: 
move.b  $F8COOA,dO 
bmi  vbwpb : 
move.b  #2, $F8CO00 
subq.l  #$01, dl 
bne  vbwpa 


; Set  index  for  delay  purposes 


,-Read  SFR  to  check  busy  flag 

,-Stay  here  until  busy  flag  goes  low 

;SFR  -  busy  goes  high  when  conversions  done 

,-Stay  here  until  conversion  complete 

,-Set  channel  2  for  next  sample 

; Decrement  index 

;Loop  through  80  times 


,-Channel  2  and  part  5  are  now  set  for  the  next  sample 
; which  is  Bus  Voltage 

vbloopa  move.b  $F8CO0A,d0        ;Read  SFR  to  check  busy  flag 

bpl  vbloopa:  ,-Stay  here  until  busy  flag  goes  low 

vbloopb  move.b  $F8CO0A,d0        ;SFR  -  busy  goes  high  when  conversions  done 
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bmi  vbloopb: 
move.b  $F8C002, (a4)  + 
move.b  $F8C0OO, (a4}+ 
move.b  #3,$FBC0O0 
move. w  - (a4) , dl 

move .w  dl, _vbus 


;Stay  here  until  conversion  complete 
Get  MSB  from  A->D 
Get  LSB  from  A->D 
Set  channel  3  for  next  sample 
Bring  entire  word  into  dl 

; store  result  in  variable  vbus 


movem.l  (a7) +, dO-d7/aO-a6 


;Pull  registers  from  stack 
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